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41. — Introduction. 


In classical mechanics a particle is considered to be a point particle, that 
is, without extension in space. For a free particle it is a well-known result 
that the angular momentum Z is conserved. 

In Dirae’s theory of the electron (in which the electron is considered to be 
a point particle) the corresponding quantity is not conserved. To preserve the 
conservation theorem a generalized angular momentum M is ascribed to the 
particle with M = L-+S, L being the quantity corresponding to the classical 
angular momentum and S a quantity chosen so that M is conserved. L is 
called the orbital angular momentum and S the spin angular momentum of 
he particle. 

By considering, instead of the motion of a particle, the motion of a fluid 


(*) Now at the Department of Applied Mathematics, University College of North 
Vales, Bangor. 
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droplet in terms of a conserved energy-momentum tensor which is not sym-. 
metric WEYSSENHOFF (!) showed that, for the droplet, it is the sum of external 
and internal angular momentum densities that is conserved. 

The former is called the moment of the linear momentum density, that is, 
it corresponds to the angular momentum of the point particle. The latter is 
called the density of spin and its existence is connected with the asymmetry 
of the energy-tensor. 

Weyssenhofi’s theory has been extended by Boum and VIGIER (?). They 
have given a physical interpretation to the theory and have shown that Weys- 
senhoff’s postulate, that the time components of the spin density tensor vanish 
in the rest system, is not necessary but defines a particular type of motion. 
Boum and VIGIER assume that certain average motions of the particle can be 
described by a. symmetric energy momentum tensor density which vanishes 
outside a space-like three dimensional limited region. Two points are defined 
within the particle—the centre of mass and the centre of matter density—and 
the spin angular momentum is shown to be the total angular momentum of 
the finite particle with respect to the centre of matter density. 

To formulate a generalized Hamiltonian dynamics for particles with spin 
we observe that for the free relativistic point particle, if the world momentum p, 
is not equal to a constant (mass) times the world velocity dx /d7 (where 7 
is the proper time) angular momentum will not be conserved. 

For particles with spin, defining the momentum by 


da 
1 ‘a Zz a 
( ) Pa dt Tr Ya 


leads to non-conservation of angular momentum for the free particle and the 
possibility of introducing spin. If relation (1) holds however the motion of 
the particles cannot be described in terms of position co-ordinates and momenta 
alone. Thus to introduce spin we must introduce extra dynamical variables 
for the particles. This will be done by considering particles as finite and not 
as point particles. 

It is well known in classical mechanics that the general motion of a fluid 
element can be considered as the resultant of a rigid body motion together 
with a deformation of the element. A rigid body motion can be considered as 
the translation of any specified point P of the body together with a rotation 
about some axis through P. 

By analogy we shall consider some point P within the particle and discuss 
the motion in terms of the motion of P and the motion relative to P. 

To do this we will consider a frame of reference whose origin is P and which 

(à 

(2 


Jo. 


WEYSSENHOFF and A. RAABE: Acta Phys. Polon., 9, 7 (1947). 


ral 
) D. Boum and J. P. Vicrer: Phys. Rev., 109, 1882 (1952). 
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moves with P (*). To start with at least P will not be specified. Let (ONE 
(u=9,1, 2,3; «—1,2,3, 4) be the components, relative to a Lorentz frame Usi 
of a tetrad of unit vectors U“ which are orthogonal and with U° a vector tan- 
gential to the world line of P in 4. The U*, will be considered as functions 
of the length È of the world line of P. These will be considered as the dynam- 
ical variables in a generalized form of Hamiltonian dynamics (8). 


2. — Generalized hamiltonian dynamics. 


Generalized Hamiltonian equations for a classical system have been ob- 
tained by MARTIN (*). MARTIN discussed a system defined by N dynamical 
variables 91, Pas... Py- 

The Poisson bracket of two functions « and v of the g, is defined by (5) 


1 ou dv 
= —— x ——— 
(2. ) (u, (i Cp RS CPs b) 


where the N? functions «,;(9,) Satisfy 


[ Œns + sr = 0, 


J 
(2-2) ter tne ao Care, 
| 25 ès + ps eqs + Los 09, . 


These are necessary and sufficient conditions for the P.b. to satisfy all the usual 
conditions including the Jacobi identity. 
The time derivative of any function w(y,t) will be given by 


du du 
(2.3) alp RE 


where H, the Hamiltonian, is chosen to represent the development of the system 


in time. 


(3) A four dimensional set of axes moving with the particle has been used to 
give a Lagrangian formulation of Bohm and Vigier’s theory see F. HaLBwaACHS, P. HIL- 
LION and J. P. Viarer: Nuovo Cimento, 10, 817 (1958); F. HaLBwacus and J. P. VIGIER: 
Comp. Rend., 248, 1124 (1959). Tetrads associated with the particle have also been 
considered by T. TAKaBayast: Nuovo Cimento, 13, 532 (1959). 

(4) J. L. MARTIN: Proc. Roy. Soc., A 251, 536 (1959). 

(5) The summation convention for repeated suffixes is used. 
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In particular 
d oH 


2, Sip auld Hide 
(2.4) di Pr (Pr ) {rs 0Ps 


Eq. (2.4) are the generalized Hamiltonian equations. 


3. — The equations of motion for the U”, 
The U", are orthogonal unit vectors. Therefore 
(3.1) CEU SOs 


Consider U a 4x4 matrix where U,,=U",. Then eq. (3.1) can be written 
as the matrix equation 
UE 


where — denotes the transposed matrix. 
Obviously UU =1 which gives the equations 
(3.2) UM UE, = dp - 


nates Orie 


Consider now the tetrads U'“,, U“ at two successive positions on the 
world line of P. These will be related by an infinitesimal transformation 


ORE, aU oO UE 


(with do,,+ d0,,=0) or 


SU", = 30,0", . 


Dividing by dé and taking the limit as SÉ +0 we get 


dues 
Be Sig LE v 
( ) a dé 0e U 62) 
where 

dos 

fn = dé ’ 
so that 
(3.4) 2. + Qi = (|); 


Eq. (3.3) will be regarded as the equations of motion for the ee, 
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If x, be the co-ordinates of P in Z then 


dés=. da_ dr, 


or 
dt, de, — 
dé dé | 
Thus 
ro Vs 
os dé © 


In our generalized Hamiltonian dynamics the variables will be the x, and 
the U“ —twenty variables in all of which only ten are independent on account 
of eq. (3.1)—with the arc length £ (of the world line of P) the independent 


parameter. 


4. — The generalized hamiltonian dynamics for the finite particle. 


In non relativistic mechanics Hamilton’s equations are 


dx, oH dir OH (r = 1,2, 3), 


"diy Ope, di Oui, 
and for the free point particle 


Ha PrPr è 


2m 


The P.b. of two functions « and v of the «,, p, is defined by 


This scheme can be obtained from Martin’s equations (with N= 6) by 


writing 

Pr = PsAtr (R=1,2,...,6; A—0,1), 
with 

Pr =, for A=; 

Pr = Dr for A1 
= 


eas 


94 J. B. HUGHES 


and choosing 
Ang = En@ Ors (S=3@-+3s; @®=0,1), 
where 


1 =— im = 1, Eo = Eu = 0. 
Similarly for the relativistic point particle, for which 


dea Bea (4 = 1,2,3,4) 
dt CPx dr Gé 


(where 7 is the proper time). 


iL 
4.1 He D> Pa : 
(4.1) Sy Pala (°) ; 
du. dv du ov 
(10) = — — , 
CL, OP x OD, OU 
write 
Pr Pa Ato Cae estos A=0,3) 
and choose 
Ces — E4098 . 


This suggests that, for the twenty dynamical variable case, we write 


(4.2a) Pr = Vine (A 0, 1, 2; 34412 2028 
with 
(4.2) ere: for Aj, 

pee Ola tors ete 


and choose 


(4.3) 


Ons E49 Ong . 


The generalized Hamiltonian for the free particle will be defined by 


Mm 
(4.4) H = DI > Pipe. 


(°) More generally, if 22 = (1/m?)p,px(= c?) we could have 


= mck (2), 


with (A) an arbitrary function satisfying F’(A),_,= 1. 
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where m is a constant and the summation is over all the Pr except p,=%,. 
x 
Thus 


(4.5) H = 4mU*,U", 


‘which is an obvious generalization of (4.1). 
With «,, given by (4.2) eq. (2.1) becomes 


a\ [ai A 
= i du. Ov du Ov du Ou 
(4.6) (u, v) = Eau ==" QU” QU“ EA Ew ATT any? 
Che OU ee wou! ats CU. aU 


(u,v = 0,1,2,3; « =1, 2, 3, 4), 


and (4.2b) in (2.4) gives 


6 da, 
(4.7) TÈ = (ay, H) = mel", 
durs 
(4.8) de AU HN me ay = 


Comparing (4.7) and (4.8) with (3.3) and (3.5) we see that 


(nes, Ut Us 
(4.9) | and 
ue eee 


Thus the P.b. is completely defined 


2 È A 
uw ov Ou Cv ou Cv 
a | Q 


7, 00, OU, 00, "OU, OU, 


© 
À 
0. 


(4.10) m(u,v) = 


This, together with 


where 
a ail 4 TJU 
‘Him UT, 


specifies the dynamics scheme for the free finite particle. 

Note that for the point particle mU°,=p, so that the generalized P.b. 
(4.10) is obtained from that for the point particle by adding the term 
Q,,(dudu!,)(00/Qu",) which will be non-zero only if w and v are functions of 
the «new» variables U*,, U?,, U*,. Thus point dynamics is the limiting case 
of our dynamics. 


Uni 
au 
+ 
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5. - The momentum and spin of the finite particle. 


Having introduced extra dynamical variables and the dynamics scheme 
they satisfy, we can now define a generalized momentum for the finite 
particle of the type given by eq. (1.1) and therefore introduce spin angular 
momentum. 

It will be assumed that the motion of the finite particle can be described 
in terms of its momentum p,, its external angular momentum L,,, and its spin, 
angular momentum S,,. These will be defined by the equations 


(5.1) Py = 4, a ears 


(4, = multipliers; we assume that a) will have the dimensions of mass), 


(5.2) Lig = D, Pg LP,» 
(5 .3) 93 =s,,U",U", (s,, = multipliers) 
with 

Sis + Su = 0 


For the free particle these satisfy the conservation laws 


a Apa 
(5.4) ar = (Ps, H)=0, 
dM, d 
(5 5) ne Gide (Lap Supe ON, 
M g=Lipt S,, is called the total angular momentum (7). 


Moreover, if we assume that the same P.b. relations hold between the Des 
Pas M,ss as hold between the corresponding quantities for the classical 
point particle; and that the a, and the s,, for the free particle are inde- 
pendent of position along the world line of P then these latter will be de- 
termined completely in terms of the 9, This means that the motion of the 
particle can be described in terms of the motion of the point P and the rota- 
tion of a four-dimensional set of axes whose origin is P. In particular it will 


(*) Similarly p, could be considered as the total linear momentum, being the sum 
of the external linear momentums açU°, and the internal momentum oUF; 
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be shown that for this set of axes we could take the set of vectors whose di- 
rection cosines satisfy the Frenet formulae in four dimensions. This set we 
shall call the Frenet frame. 


6. — Determination of the a,,, s,,, 


We have the following equations to determine the Was Say 


(6.1) (Px H) == 0 
(6.2) (M,,,H) = 
(6.3a) (©) Pp) = da 
(6.3b) (Gr Cp) = (Pas Pg) = 0 
(6.4) (M3, M uv) = M, Og» =r M,,0,, "A Mo dei Mon oF a 
i) (p,, H) = 0 
OP, OH 
—— —— — ( 
Qu U 00" 
that is 
Q,a,U", = 0 


Multiplying by U°, and summing over x we get 


(6.5) 26%, = 0 


no po 
ii) (M,,,H) =0. 
(@,P,—%pD, + 8,,U",0",, H) = 0, 
UT Dg el rai pat Spl seal gl pe UO) 0, 
0,58 9 4 Der D, vu UT = Gi UR Uo ae U#,U”°,) : 
Multiplying by U%,U°, and summing « and f we get 
(6.6) Quy 89 — Lo Sry = 1009 — 00 oy » 


ili) a) (0,0) = (Pas Pp) = 0 
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These follow automatically from the definition of the P.b. For example, 


since 
OU AO) 
a oP T y a B cel 
(6.7) m(L ds U 8) ai Qe: SIT MAT TC Duo Sap ; 
CU ss Us 
then 


M(P ys Ps) = ma,a, (UT; DS) Fe BO, 2 iy Ong = 0 ? 


because of the asymmetry of 2,,,. 


This merely establishes the relation 
(6.8) Go itn, 


where m is the constant appearing in the definition of H, and thus has the 
dimensions of mass. 


iv) Moy) = Mug, + Mp ban — Mas Spu— M pu 8 


op By a on” 


Because of eq. (6.3) the above is automatically satisfied for M,,= L,,: There- 
fore it reduces to the condition 


84, 8,04, Up VU) = SU UT, dug + U%s Ud = UU, 8g, — VU Sap) - 


u 


Using (6.7) and equating the coefficients of the 6’s appearing we get in 
each case the condition 


(6.9) ae — Miro. 


Eq. (6.5), (6.6), (6.8), (6.9) are sufficient to determine te a, s 
the 2, and m, 


pv 


in terms of 


uv 


This will be demonstrated more easily if we consider 2,,,, $,,, a,—which 
formally appear as 4-tensors—as two six-vectors and a scalar plus a 3-vector, 
respectively. 


Let 
(a) a, = [dc a] 
(6.10) | (b) 92, =[w; v] 
| SRE sal 
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so that 
Or == (00% 
(6.11) 
Q ix = Ex Vi 
| Si = li 
(6.12) 
| Six = Enid: 


where Latin suffixes have the values (1, 2, 3), and ¢,,;, is the completely anti- 
symmetric tensor of rank 3. 
(6.5) is now equivalent to the equation 


(6.13) AKU = 0,0" 


Eq. (6.5) is satisfied identically when y=9=0. When 0=0, y—% it 
becomes 


(6.14) qxw+Ixv=a 


and when 0=%, y=j (#1) 


(6.15) Ixw+qxv=0. 


Eq. (6.9) is satisfied identically when t=o. When t=0, o=k it be- 
comes 


(6.16) (l-w + q-v)l — (l-q)v = ml 


Andivhentr 1 0) (>= Kk) 


(6.17) (lew + q-v)q— (l-qyo = mq. 


We thus have the five equations: 
aXv = GW , 
qxXw+ixa =a, 
Ixw+qxv=0, 


Al+ pov —= ml, 


Aq + uw = mq, 
where 


A=lw+q'v, 


22 = l:q, 
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Qi = Wi Q;; = Ein 4 
|’ so that 


as (6.22) @? = Qi Boi , v = Ho OI 29,9 


uv? 


% 


DE == mU® + a;U', 


and 
Mii, À Sap = 8y,U",U"s ) 
PE. 
(0-23) a, È PEN à 
Sop == 0 
2) aon 


Then p,p,=m*+a?. To evaluate a?=a,a; we get from (6.19) 


4 
at pas | m2 m2 ©? 
nes (6.254) a? = sa (DX Ww)-(vXw) = E [v2@2 — (v-w)2] = m? = 
“A ; | 
LM 
Ar; 


RCE EE Pa 
xh ° a #7 ee FF eee AR Bee? y 
= 7. è Se tu, 

We, oa ete GENERALIZED HAMILTONIAN DYNAMICS ETC. - 1 


» dgr = 


52 


Thus 
(6.25b) Pas = Mm? (i Hr 2) =m? - - 
| v? -v? 
Also 
48 apap = 35,8 ni 38,55% * 
but 
1m : 
. (6.26) 5585 = 3 © 2,9, — = 7 
whence 


Eg. (6.18) is a relation to be satisfied by the 2, 
(6.27) Ein Doi Qn = 0 . 


_ It can be proved (see Appendix) that the ©,, satisfy also the relation 


(6.28) 92,0 yy Qu = Ra 


ap pv 


where 7? =o? + v2=39,,9,,. The s,, satisfy a similar relation 


mi 
(6.29) Sau woop = — Ta Bap - 
From (6.29) follows that the S,, satisfy 
_ (6.30) 8 hey fot pe = — KS - 


ba 


7. - The equation of motion for the point P. 


Substituting the values for the a,, given by (6.23), in (5.1) we have 
m 
Pa = mu, à à pr Don us a. ’ 
m = 
= mU*, * pr (Lou Dur U eer 0 yp U? 5) , 


wo? m = 
= m(1 +2) U°, Wise onl : 


CL A) 
Pa 
Do 


= 
PL 
“ 
PS 
a 


i * 
= 
or 
LA 


+ = 


o x 
a, J 
‘ = Ca 
‘7 è 
ae tae 1e 


- 
/ ; 
‘bee 
din 


MOTI 
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where dot =d/dé. Since the U*, do not depend explicitly on & it follows 
from (4.11) that 


COE ke È Li ii 


(03 


Hence 


Aon ee Hers 
(7.1) Bae ar cae 


Eq. (7.1) is a second order differential equation for U°,: Since US =, where 
æ, are the co-ordinates of P, it can be regarded as the equation of motion 
or P. 

For the free particle p, is constant and eq. (7.1) is easily solved giving 


fi toi 
(7.2) U°,= A, cos AE + B, sin SEPE ARE 
where A,, B, are arbitrary constant vectors. 
This shows that the motion of the point P consists of an oscillatory motion 
Superimposed on a constant rectilinear motion. The constant rectilinear part 
is proportional to the total momentum and is 


MPa “{ M\ Pa 
Pa Ds, a M’ 
where M=Vp,p,=€ times the proper mass of the particle. The oscillatory 
part has a frequency À where 4° = p,p,/38,,8,, which will be very large if 
the spin is small compared to the momentum. Eq. (7.2) is similar to Dirac’s 
solution for the motion of the free electron (8). 
The result (7.2) is a particular case of the application of equation (6.28) 


to the equations of motion (3.3) for the U”. Multiplying (6.28) by U°; and 
summing over f we see that each of the U*, satisfies the equation 


DQ 82 yy $248 U = — VQ 0°; 
which can be written as the differential equation 


(7.3) (D?+ 42D)U",=0, 


where D=d/dé. 


(*) P. A. M. Dirac: The Principles of Quantum Mechanics (Oxford, 1947), § 69. 
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Hence each of the U“ has the form 
(7.4) UF = AH cos AE + B! sin ZÉ CH 
PAC EG SULA == CN, 


where A”, BY, C“ are arbitrary constant vectors. 


8. — Comparison with the Weyssenhoff and the Bohm-Vigier theories. 


In this theory as in the Weyssenhoff and Bohm-Vigier theories the motion 
of the free particle is described in terms of a velocity, a momentum, and a 
Spin angular momentum, for which a) the momentum is not parallel to the 
velocity, and b) the spin angular momentum satisfies an equation which in 
our notation is 


(8.1) SUD "Up, 


WEYSSENHOFF assumed as well that the inner product of the spin angular 
momentum and the velocity is zero. BoHM and VIGIER showed that for a 
particle, for which a) and b) above are satisfied, an additional assumption 
is necessary to define the motion completely. Such an assumption could be 
the one assumed by WEYSSENHOFF. Thus Weyssenhoff’s theory is a particular 
case of the Bohm-Vigier theory representing one possible motion for the par- 
ticle. But in the latter theory other motions are possible depending on the 
assumption made to close the system. 

It follows from (5.3) and (6.24) that in our theory 


m2 Gale 
(8.2) Sa LISA asa ye Qin DE, J*, US: == (ll), 


Thus Weyssenhoff's condition is satisfied. It is not, however, an a priori as- 
sumption but results from our choice of the quantities CPE since we took 
So: —0. The s,, we determined by equations (6.1) to (6.4). Of these eq. (6.1) 
and (6.2) occur also in the Weyssenhoff and Bohm-Vigier theories and (6.3) 
gave no conditions on the s,,. Thus (8.2) results essentially because of the 
assumption regarding the P.b. (eq. (6.4)) of the components of the total angular 
momentum tensor. 
Consider now the pseudovector s, defined by 


== 7 dea yf 0 
(8.3) Su Ci 2 Euvof Sap Le fi; So U vy? 
where Si, is the dual tensor of S,,. Obviously s,0°,=0. Also 
a ae ee AG di 
(8.4) 884 = 39 89 489 
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that is, s, has the same magnitude as S,;. It is also constant (as in the other 
theories). 
From (6.24) 
cae le Ti 0 
Sn Fazi pve 8 ag Ut 0e sU v 
The quantity &,,,gU',U';U°, is obviously orthogonal to U*,, U?,, U°,, and 
hence we can write 


Ti k 
(8.5) nus Ut, U%,U, = erin U", 
This gives 
oy ya hgh les Ge 
(5.6) SE on D SijEoiik Ÿ n 3 
m k Tk 
lg A Bi 7 
(3 1) = Spi En Qul TA 202 VEU pu 
Then 
m E 
Sin = dp Eijk Qi Qua U IDE) 


m = ap 
= — [gw 9250 U iy Bai 835 Dis QU id . 


2% «2 
‘The first term is zero on account of (6.27) and the second since 
(8.8) Es Qui = £0 Dy Dire = 4 VrErimUm = 0. 
Thus s, is constant. Following Bohm-Vigier we call it the spin pseudovector. 


For future reference we note that the inner roduct of 5. with its dual 
È Pp op 
S is zero 
af ? 


2 

C I a fl 1 à 

DA) HER Even Q,;2 m Epyus U pe pU “ DE, ’ 
vi 0 

aa LIO U ra y (9) , 


which is zero on account of (8.8). Hence 


(8.9) ESSA 
and in particular 
(8.10) SES EU 


(9) Seyi pU Uno = € Ge see (8.5) > 


jlmk y 


U°, (since 7, !, m, #0). 


Ejlm 
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9. — The Frenet tetrad. 


If we write 


then 

(91) Tem Ur 
and 

(9.2) r,U°,=0. 


With the vectors U°,, 7, we can form with S,, the following tensor and pseudo- 
tensor quantities 


(9.3a) ISAS PRE beet be 
(9.3b) Sr Dane 


Of the quantities (9.3a) the first is identically zero, and the other is the spin 
pseudovector s,, 


à m 
$s i EE k 
Sa = Sap U BE 203 v,U ae 


Of the quantities (9.36) the second is identically zero 


gl : 
Pages | MEN j ke a 
Sp + dp? 27.4 ExpuyU mo U's 0 
1 
0 
re Lo QU œ ? 
ye 
1 3 
= Oe, Us, = 0, 
vp? be 


and the other is a quantity which we shall call q, 


1 
(9.4) qy = Sya lp = pe D PAU , 


m i 
i x Tie de ri 
= De Exit; U AE" ae v2 i U a * 


Uni 
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The four quantities U°,, q,, 7,, 8, are mutually orthogonal 


9.5 
Sa Ur Ur UV 3 == 0% 
Vu". = Sopra = 9%; 
(9.6) 
Ibe Se pre SNE AU 
n m Seats m 
(9.7) Ps a Ue a MEI n A.V, = 0. 


Let V*,, V1,, V?,, V3, be unit vectors in the directions of U°,, —q,, 1,38 
respectively. Then 


a 


7 70 
J ioe =U œ 
vire oi Ua 
(ot ali \/O# 9 
9.8 ; 
(0.8) DOD 
>. Tt 
V3 ST o U & 
œ Vv ’ 


and. 


i, TU Ti Ti STI 
Ve. => Ou (E x OU — US = 4/02 J me 


Vi, = a DOUÉ = oe (osent, Lo QUE Vor V9, + 4/22, , 
using (6.254) 
Dé = 2 OU, = vat (EU, + ©, 0°) = — 4/8 Vi, , 
Vi, = 0 


(since V*, is parallel to s,). 
Introducing quantities o and o defined by 


(9.9) 


i) 
~ 


1 
AE — = 4) 
o 
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the equations for the derivatives of the V“, are 


Vie Saves 
a 0 XI 

7 ] ] 2 

SIA +-< Vi 
(9.10) 0 lo) 

SE Lee 

V = = 5 J 5) 9 

Vie= 05 


which are the Frenet equations in 4-space (1°) for the direction cosines of an 
orthogonal tetrad of vectors moving along a curve which has a constant tri- 
normal, that is, zero tilt. V1,, V?,, V3 are the direction-cosines of the normal 
(or principal normal), binormal, and trinormal to the curve, respectively. 1/0 is 
the curvature; 1/0 the tilt of the curve. We shall refer to the orthogonal frame 
specified by the direction-cosines V“, as the Frenet-frame. 

Since they are orthogonal unit vectors the V”, satisfy 


(9.11) NET mea LEE 


op * 
Moreover, since they are linear combinations of the U“, we can write 


(9.12) V4 =0,,U",. 


py 


Regarding V“, (for w varying, « fixed) as a column vector V,, (9.12) can be 
written as the matrix equations 


(9.13) V,=0U 


x a? 


where, because of the orthogonality relations (3.1), (9.11), O must be an 
orthogonal matrix 


(9.14) 00 = 00 =1. 


(This is easily verified using eq. (9.8), (6.18), (6.19).) 


(19) See A. R. ForsyTH: Geometry of Four Dimensions (Cambridge, 1930), § 164 
whose terminology we use. 
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10. — Invariance of the dynamics scheme. 


Our basic equations are 


i) the definition of the P.b. 


and ii) the definition of the Hamiltonian 
Sami" Uke 
It is easy to show that for transformations of the type (9.12) these equa- 


tions retain their form. 
Eq. (9.12) can be written 


(10.1) US = QE ’ 
so that 
OU", OV", 
av, cs O uv ave O60 wp 
Hence 
ou 2 ou OU's 0 ou 
ave OU”; dee Lia 
or 
, Ou ai 
(10.2) BU. = Ono Ta 


(10.3) (u, v) = + n se Sa “ PO a ee 
where 

(10.4) DOr Om 

or in matrix notation 

(10.5) Q' = 020 
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and from (10.1) 


m 


H aS m Ue Ut BRAS ye yt | 
(10.6) DUR AE ep ET ETATS EE 
PRO Ve. 


Thus the dynamics scheme is invariant to an orthogonal transformation of 
the moving frame or reference. Consequently the motion could be discussed 
in terms of the Frenet frame (11). It will be seen from (9.10) that by doing 
so our equations will be simplified since then only four of the quantities QU 


will be non-zero. 


11. — Concluding remarks. 


In this paper a generalized Hamiltonian dynamics for classical finite par- 
ticles with spin has been introduced. The main assumptions made are: 


i) that the motion of the particle can be described in terms of the 
position co-ordinates of some point P within the body, and the direction-cosines 
of an orthogonal tetrad of vectors forming a moving frame of reference whose 
origin is P; 


ii) that the motion of the particle can be described in terms of its mo- 
mentum, its external angular momentum, and its spin angular momentum, 


Accepting these assumptions it was found that the Frenet frame could be 
used as the moving frame of reference. 

In a subsequent paper the motion will be discussed relative to the Frenet 
frame. This results in great simplification since only four of the quantities Li, 
will be non-zero. The motion will also be discussed in terms of a frame obtained 
from the Frenet frame by a rotation leaving V! and V* unchanged. The re- 
sulting scheme has a quantum analogue yielding a wave equation in which 
the mass term includes an operator involving internal co-ordinates of the 
particle only, thus introducing the possibility of mass quantization. 


(1) Use of the Frenet equations in obtaining spinor equation of motion for classical 
point particles was made by GURSEY. F. Giirspy: Nuovo Cimento, 5, 784 (1957). 
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APPENDIX 


Proof of Relation (6.28). 


First of all we note that because of (6.24) and (6.11) 


DEI 
ORO == QD Qe Qo5 fi a (to d;0;) ; 
‘0 
p2 Me. 
= —[axw},, (kK A, 9), 
o 
= [(vXw) +o], = — 0%. 
Hence 
GAZ) 2 9.21.21; ala Q 5,82 9203 Set Se wQ;; O 
From (6.19) 
871810500 = — MS>5 n 
but since s,, = 0 
Se = — MS; ; ° 
Substituting for s,;, from (6.23) 
a a a 
# OTOTOS — = mo, = . 
A "2 p2 
This result together with (A-1) gives 
LOR KON AB, —- ONTO + QUOTO = — v20,, —= oo Oe , 
that is, 
(A-2) DQ unhdy; = — (@? + v)Q,. 
Now consider 
(A-3) Fe eens Fr Die + 90,,9,9n ’ 


Æ 2 7 
= — wD) + E 16€ 5m 0 Um ro , 
=— 002; = 0,0;0; — 0°72 


LIS) 
= — (04 09)Q. 


Therefore from (A-2) and (A-3) it follows that 


VASE) Qing gun (0° + 0) Qyp n 20, si per (32,20) Lap È 
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A similar relation holds for the s,,: 
as à a a Aan Le > 
(A-5) Sous pvS»B nuit (ES pSuv) Sop 


Since s,,=90 it is obviously true when Alter iGO OlN ia 0s Otherwise 
the left-hand side is 


SiuSuoSvi = Sid ki 15 3 
> His Oe com O..° 
by (6.23) = RE Sines m9 1j 9 
. do | 
by (6.9) = (— G8i3) ; 
m°? 
= — Soi S 55 = — (4 Suv Say) Sis , 


Hence (A-5) is established. 
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Further Experiments with PAPA. 


A. GAMBA, L. GAMBERINI, G. PALMIERI and R. SANNA 


Istituto di Fisica dell Universita - Genova 


(ricevuto il 27 Marzo 1961) 


Contents: Introduction. — 1. The checkerboard experiment. — 2. The 
music experiment. 


Introduction. 


The PAPA machine, built in our Institute, has been described elsewhere, 
and its ability to recognize simple geometrical patterns has been demonstrated (*). 
In the experiments reported here we wanted to test the ability of PAPA in 
recognizing more sophisticated « cultural» patterns. By this term we mean 
patterns whose significance is given by « hidden rules » which are not obviously 
apparent to a superficial observer. For example, in experiment No. 1 (Sec- 
tion 1 below) we presented a checkerboard in which holes connected by rook 
moves and knight moves respectively were given to the machine for recog- 
nition; in experiment No. 2 (Section 2 below) we coded a few notes from 
Bach’s chorals and from «fake» music made by one of us (L. GAMBERINT) 
and we asked the machine to learn these two classes of patterns. 

The results turned out to be exceedingly encouraging and we think that 
this proves the «intelligence » of our machine. Given the appropriate devel- 
opment by increasing the number of association units, PAPA can be a de- 
cision making machine that, for instance, by looking at an electrical wiring 
diagram can decide whether it is a «good» one or not, given the pertinent 
instruction by examples. Or it can become a weather forecasting machine when 
shown examples of « good » and « bad » forecasting from a collection of physical 
data. In other words, PAPA will behave «as if» it had guessed the rules. 


(*) A. GAMBA: Proc. IRE, 49, 349 (1961); G. PALMIERI and R. Sanna: Methodos, 
to be published. 
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FURTHER EXPERIMENTS WITH PAPA 


1. — The checkerboard experiment. 


113 


In cards representing a checkerboard, such as in Fig. 1, we punched 12 holes. 
In the first class the 12 holes were connected continuously by rook moves and 


in the second class by knight 
moves. We punched a total of 
300 cards (150 of each type) 
all different. One hundred cards 
of each type were shown to 
PAPA during the learning pe- 
riod and the remaining 50 of 
each type were kept apart to 
be shown later to the machine 
for recognition. Again we stress 
the fact that these were all 
different from each other and 
from those shown during the 
training period. 

As described elsewhere (1) 
our PAPA has 79 A-units. In 
order to see how it would be- 
have with a greater number of 
A-units, one could change the 
random masks in each A-unit 
and repeat the experiment, 
combining the 
those obtained before. 
simple approach is to rotate 
the PAPA (or, equivalently, 
the patterns being shown) and 
consider it as another PAPA 
with another (rotated) 79 A- 


answers with 
A more 


units. In our experiments we rotated the patterns by 
could repeat the learning curves as the number of 


core | oo | arma 
ae EHH Aas 
pte IE (E 
dir | EE 
49 12 26 =e 43 
FREE | eme | cure | porc a Se 88) 
a A | sa fH Hee || ‘at : 7 \ 
He | eee || eee | OS | ca 
Mieli | «"s | [sans | eee | n 
16 -49 4 38 27 
Fig. 1. — Rook moves (top) and knight moves 
(bottom) recognized by Papa. The recognition 


factors (negative when wrong) refer to 12x79 


equivalent A-units. 


LANDES TRS NAT Rea 
Fig. 2. — Learning curves for rook moves as the 


79 to 12> 
The abscissa gives the number of A-units as a mul- 
tiple of 79. The ordinate gives the number of pat- 
terns recognized correctly (D), uncorrectly (x), and 
given even (©) out of the 50 examples presented. 


number of A-units increases from 1 > 79. 


30° each time and thus 
A-units increased from 
79 to 12x79. The results are 


and 3 


Fig. 2 


summarized in 
for the rook and knight moves, 


SPACE 


463 


Fig. 3. — Learning curves for 


knight moves. Same notations as 


chav, AT eve 


10° 11 12 
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respectively. The recognition factor (given as logarithms to base 10) for the 
equivalent 1279 A-units ranged from a maximum of 70 down to — 49 


with an average of 22.70. 


2. — The music experiment. 


A few notes from various Bach chorals were coded into cards, so that each 


card corresponded to a simple musical phrase. The coding is more easily 


Bae 


| 
| 
| 


È °° 1% 


1. 


Fig. 4 — Examples 


showing how 
coded for experiment no. : 


music was 


explained by giving two exam- 
ples (see Fig. 4). The sign at 
the bottom right stands for the 
pause at the end of the tune. 
By such a device the total area 
of each phrase was then always 
the same. 225 tunes such as 
these—all in the same key— 
were selected by 
(GAMBERINI). 
Another 225 tunes for the 
other class were then prepared 
as examples of «fake » music. 
They were made so that they 
sound «bad» when played on 
the piano keyboard and com- 
pared with the «good» Bach 
music. They were, however, 
chosen so that it was not im- 
mediately apparent from the 
coded cards which music be- 


one of us 


longed to which author (see Fig. 5). We wanted to see whether PAPA was 
able to form the concept of « Bach music » vs. «non-Bach music », whatever 


this might mean. 
50 examples of each class 
were then selected at random 


Fig. 5. — Examples of Bach music 
(bottom) and non-Bach music (top) 
recognized by PAPA. The recogni- 
tion factors (negative when wrong) 
refer to 12 x 79 equivalent A-units. 


me LT CT mn, at 
tim (i ili | 11 Im 
21 9 PE 13 

Pa lug | x a 
ae | ms i DE ; a i pei 
11 
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and put aside to be shown to PAPA after instructing it with the remaining 
175 examples. 


The results are given in Fig. 6 and 7. 


Fig. 6. — Learning curves for 
Bach music. Same notations as 
in Pig. 2. 


Fig. 7. — Learning curves for non- 
Bach music. Same notations as 
ma ira Dr 


SS We ih 12 


After this experiment was completed we discovered that, on the average, 
the Bach music had a shorter pause at the end than non-Bach music. This 
has certainly introduced some bias, whose amount should not, however, sub- 
stantially alter the result of our experiment. The recognition factors here 
ranged from a maximum of 71 down to — 57 with an average of 14.42. 
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Inelastic Collisions and Threshold Effects. 


L. FONDA 


Institute for Advanced Study - Princeton, N.4J. 
Istituto di Fisica dell'Università - Palermo (*). 


(ricevuto il 5 Maggio 1961) 


CONTENTS. — 1. Introduetion. — 2. The many-channel problem. — 3. Time- 
independent scattering formalism. — 4. Low-energy behaviour of cross 
sections: a. Endoergic case. b. Exoergic case. c. Elastic scattering. 
d. Three-body final channel. — 5. Threshold effeets. — 6. Coulomb effects. 


1. — Introduction. 


In a number of recent theoretical papers the many-channel problem has 
been extensively discussed with particular emphasis on the anomalous energy 


dependence which is observed in the scattering and reaction cross-sections at 
the threshold of a new inelastic process (1). It is the purpose of this paper 
to review such phenomena by introducing a new method which is based 


*) Present address. 
A 
G. 
G. 
A. I. Baz: Zurn. Eksp. Teor. Fiz. 33, 923 (1957); translation: Soviet Phys. 
6, 709 (1958). 

Re 

R. G. Newton: Phys. Rev., 114, 1611 (1959). 

. FONDA and R. G. Newton: Ann Phys., 7, 133 (1959). 

. Fonpa: Nuovo Cimento, 13, 956 (1959). 

. G. NEWTON and L. Fonpa: Ann. Phys., 9, 416 (1960). 

. M. Detves: Nuel. Phys., 9, 391 (1958-59). 


P. WIGNER: Phys. Rev., 73, 1002 (1948). 
Breit: Phys. Rev., 107, 1612 (1957). 
Breit: Handbuch der Physik, Band x11/1 (Berlin, 1959), p. 1. 


G. Newton: Ann. Phys., 4, 29 (1958). 
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on a thorough investigation of the properties of the complete Green’s function 
which describes the multi-channel system. The threshold behaviour of the 
Green’s function which gives the propagation of probability in the newly opened 
channel, will prove to be the origin of the anomalous effects in the complete 
Green’s function and consequently in the various elements of the transition 
amplitude. The method is able to describe both non-Couiomb and Coulomb 
effects. In the former case, anomalies of the kind of a cusp or of a rounded 
step show up in the old cross-sections at the threshold for the new inelastic 
process, due to the sudden removal of flux from the incident beam at the onset 
of the new inelastic cross-section starting with an infinite slope. On the contrary, 
in the latter case no cusp or rounded step appear; instead in the case of Coulomb 
attraction a finite jump is expected in the various old cross-sections due to 
the non-zero value of the new inelastic cross-section at its own threshold (7). 
No anomaly is observed, however, at the threshold for a reaction leading to 
three (or more) particles, apart from an unobservable discontinuity in the 
second derivative of the old cross-sections. If consideration is given to re- 
actions whose final three-body channels have at least one particle in common, 
a threshold effect then appears as a function of the momentum of the common 
particle with the total energy kept led) 

To introduce the reader to the many-channel problem, we will discuss in 
Section 2 a simple example, already considered in ref. (#7), which is an 
extension of the method of the distorted waves developed in the theory of 
atomic collisions (23). In this way we will be able to define our formalism and 
to understand various features of the many-channel problem. 

The general time-independent scattering theory which covers the case of 


(11) R. K. Aparr: Phys. Rev., 111, 632 (1958). 

(12) A. I. Baz and L. B. Okun: Zurn. Eksp. Teor. Fiz. 35, 757 (1958); transla- 
tion: Soviet Phys. JETP., 8, 526 (1959). 

(13). J. D. JACKSON and H. W. Wy LD, jr.: Nuovo Cimento, 13, 85 (1959). 

(14) L. Fonpa and R. G. NEWTON: Nuovo Cimento, 14, 1027 (1959). 

(15) L. I. LAPIDUS and CHou KuaNne-cuao: Zurn. Eksp. Teor. Fiz., 38, 201 (1960); 
translation: Soviet Phys. JETP, 11, 147 (1960). 

(16) L. I. LAPIDUS and Cou Kuane-cnao: Zurn. Eksp. Teor. Fiz., 39, 112 (1960); 
translation: Soviet Phys. JETP, 12, 82 (1961). 

(17) L. Fonpa and R. G. NEWTON: Phys. Rev., 119, 1394 (1960). 

(18) L. I. LAPIDUS and CHOU Kuane-cuao: Zurn. Eksp. Teor. Fiz., 39, 364 (1960); 
translation: Soviet Phys. JETP, 12, 258 (1961). 

(19) P. Bupini and L. FONDA: Phys. Rev. Lett. 6, 419 (1961). 


(29) G. BreIr: Phys. Rev., 69, 472 (1946). 

(21) M. Cini and S. FUBINI: Nuovo Cimento, 2, 75 (1955). 

(22) H. FESHBACH: Ann. Phys., 5, 357 (1958). 
(23) N. F. Morr and H. SW. Massey: Theory of Atomic Collisions, Il Ed. 


(Oxford, 1949). 


lod 
© 
= 


118 L. FONDA 


Coulomb forces and many-channels, is considered in detail in Section 3. In 
Section 4 the low energy properties of scattering and reaction cross-sections 
are reviewed. The threshold effects are discussed in Sections 5 and 6. A dis- 
cussion on the importance of such phenomena for the determination of rela- 
tive parities and spins of the reaction products, of scattering phase shifts, of 
cross-sections for processes which are often not feasible experimentally, and 
in the search for new particles, is given in Section 5. 


2. — The many-channel problem. 


We will consider in this section the collision of a particle with a target capable 
of a certain number of excited states. The target particle will be supposed 
to be composed by A particles, while the incident particle will be taken, for 
simplicity, to have no internal structure. The interactions between the À +1 
particles will be described by a set of potentials. In this example only the 
possibility for scattering (elastic or inelastic) of the incident particle on the 
target will be discussed. Also effects due to the identity of particles will not 
be taken into account. All the A-+1 particles will be, for that matter, as dif- 
ferent from one another. Such a description would be, therefore, an approxi- 
mation to the actual situation of, for instance, collision of a nucleon with a 
nucleus since both the possibility of exchange scattering, of production of 
deuterons or «-particles, and the requirement of an antisymmetrical total wave 
function will not be considered. 

Let #' be the total hamiltonian describing the system target incident 
particle. We write for it 


eu HET ET +H AV,» 


where H, describes the target nucleus in its center of mass frame of reference, 
T, is the kinetic energy of the incoming nucleon, 7° is the kinetic energy of 
the center of mass of the nucleus. V,, is the sum over all the elementary po- 
tentials between the incident particle and the A particles of the target nucleus. 

Introducing co-ordinates ro, Fi, T2,..., r, measured, for convenience, from 
the center of mass of the system A, we separate the motion of the center of 
mass of the A+1 system and obtain the Schrédinger equation 


(2.2) HIT Rae rs.) = EV(r, Fip Tale 


he 
(2.3) «x = "Ag Ht Von, 


2 Lo 4 
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where 2, = mom,/(m+m,) is the reduced mass for the incident particle. 
Let «, be the complete set of eigenfunctions of H ,: 


(2.4) a. 


H, will have in general a discrete plus a continuous spectrum. The discrete 
spectrum represents the various excited states of the target nucleus which, 
for our collision problem, will characterize the «discrete or two-body chan- 
nels ». The continuous spectrum represents the various possibilities for the 
nucleus to break-up, i.e. the realization for our collision problem of « continuous 
(three or more body) channels». In the process of collision of the incident 
nucleon with the target nucleus in a certain state u,;, the potential V,, will 
induce transitions between the various energy levels #, of the nucleus and 
finally the incident nucleon will leave the nucleus in a state w;. We will say 
that the reaction channel i > channel f has occurred. In so far as channel f 
is a two-body channel, the kinetic energies of the two final products will be 
completely specified once the total energy £ is given (the initial channel % is 
a two-body channel for obvious experimental reasons). On the contrary, if 
channel f is a three or more body channel, we will have a continuous range 
of energies for the three or more products present in it for a given total 
energy E. The continuous channels will be neglected in the further develop- 
ments of our example. This will be a good approximation if the energy £ is 
far from the threshold for break-up of the nucleus. 

In the absence of the exclusion principle, we can write for the total wave 
function Ÿ the following expansion: 


(2.5) Wir tages tg) ae D Ulis cs 74) Palo) - 


i 


Introducing such an expansion for Y into eq. (2.2), we get the following set 
of coupled differential equations: 


(2.6) | (E — Hii)yi = > Vis Ps 
j 
where 
Di 
2.7a LEE, 
= 2U: 
(2.7b) Vis == (Wi, Vo Wi) ’ 


where u, is the reduced mass in channel è and in our simple example is in- 
dependent of 7. po is the momentum operator of the incident particle. V;; 
depends only on ry but can in general be non-local and energy-dependent. 


oa 
ce 
nul 
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Eq. (2.6) can be written in matrix notation 
(2.8) (E — H,)p=Vy, 


where the potential V is a square matrix whose elements are given by (2.7b), 
H, is diagonal with elements given by (2.7a), and y is a column matrix with 
elements y;. Hermiticity of V,, implies hermiticity for the matrix V. If the 
potential V,, is invariant under that particular time reversal transformation 
which reverses only the motion of the incident particle considered as iso- 
lated from the system A, then from invariance under the usually defined time 
reversal it follows that the matrix V is symmetric. 

It is sometimes useful to introduce projection operators P, for the matrix 


representation given above: 
(2.9) IE. Pia 04; P; ’ PA i= eh HE 


If we consider the solutions of the homogeneous equation associated to eq. (2.6): 


h2k3 


> 


ir 


(2.10) H,@ = | 3h E, oo = Eq, 


where k; is the wave number in channel 7, we can express the projection oper- 
ators P,; for example in the form 


(2.11) (Pise = On Os > P02) (Per | ’ 


where the symbol > stands for the integrations and summations over all the 
continuous and discrete quantum numbers necessary to specify gf) completely. 
Eq. (2.10) can be written in matrix notation 


(2512) (E—H)g =0, 


where gj’ is a column matrix whose elements in the co-ordinate representation 


are given by 

(2.13) Per (Se, r= d,(2x) ? x, exp [tk;-r], 

where &, indicates the momentum, the spin and its 2-component in channel è 
Bi = (ki, 8, v) 


and y! is the appropriate real spin function in channel è. 
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All the collision properties are simply expressed in terms of V and y; we 
have in fact for à reaction leading to the channel f 


(2.14) (97 Us, Ve) = (g®, Vy) 4 


Eqs. (2.8), (2.12) and (2.14) show that one can really consider y as the 
wave function and H,+-V as the hamiltonian for the many-channel system 
in substitution of the complicated Y and #. In the following we will take 
over these arguments to the description of the collision of any two systems. 


We will base our investigation of the many-channel problem on the properties 


of the general set of eq. (8) that we will assume to have been introduced in 
one way or another for any kind of reactions, even for reactions leading to 
three or more particles (*), in such a way that the collision properties be de- 
fined by the right hand side of eq. (2.14). For simplicity the potential V will 
be taken to be energy-independent. 


3. - Time-independent scattering formalism. 


In this Section I will consider in some detail by following standard meth- 
ods (295) the general time-independent scattering formalism which applies 
to the case of modified Coulomb forces and many channels. 

We separate the interaction hamiltonian matrix V into two parts: 


(3.1) V=Ve+V', 


where V, is the diagonal matrix which describes the Coulomb interaction in 
the various channels (we assume that Coulomb-like tails exist only in the 
diagonal elements of J’): 


(3.2) Vor) = 44; 42;€7/7 - 


The intégral equation for the state vector y reads as follows (**) 


: (3.3) yp? = g + G.V'y®, 


(*) A generalization of H,, V and y is of course needed in these cases. 

(24) C. MoLLER: Kgl. Danske Videnskab. Selskab, Mat.-Fys. Medd., 23, No. 1 (1945); 
B. A. LippMANN and J. ScHWINGER: Phys. Rev., 79, 469 (1950); M. GELL-MANN and 
M. L. GOLDBERGER: Phys. Rev., 91, 398 (1953). | ‘ 

(25) M. H. Hurt, jr. and G. Broir: Handbuch der Physik, Band x11/1, (Berlin, 1959), 
p. 408. | | | 

(**) Equation (3.3) does not lead to the peculiar amplitude renormalization (see 
ref. (26)) one has to perform when the free wave po and the free Green’s function Gy 
) 


û n ‘o (4) ‘ 
are taken in place of the modified pg, and (,. 


9 - Supplemento al Nuovo Cimento. 
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where the Coulomb Green’s functions G, is given by 
(3.4) G, = (H+ te — Hy— V.) 


with the usual e>0 device. Like G,, both y° and g in (3.3) satisfy an 
outgoing wave boundary condition. Sometime we will place a superscript (+) 
to these wave functions. The superscript (7) indicates that the incoming wave 
packet is in the i-th channel. In the co-ordinate representation the Coulomb 
wave function g® is given by 


a ; k,- A k; 
(3.5) DOUÉ, r) = 0,,(2x) ty D (21 + 1) exp fio] P.(* = 3 < 7 
1 il | vy 


where 
o = arg (+1 + ini) 


(3.6) 
Mi = Ailk; = Zi; 22, €7/N0; 


and v; is the relative velocity of the two particles in channel è. F',(@) is the 
real Coulomb wave function regular at the origin in the notation of HULL 
and BREIT (2°): 


exp [— (x/2)n1 


Min, v44(200) 


with 
be 1 2m (tae) 
3.8 M; 2 210) 11 1 4 — 0 + RCE 
( ) in.143(200)/(2î0) RE] CF (1 + 1)(2 + 3) & 
Asymptotically we have 
(3.9) Lim F,(o) = sin (0 — n log 20 — Si + 21) 5 
ee a 


It is useful to recall also the definition of the irregular Coulomb wave 
function F,(0) —1G;(0): 


(3.10) HP (0) = F:(0) — iG;(0) = exp 


n n exp [to] (— 1)'*1W in, i44(— 210), 
which asymptotically behaves like (*) 


(3.11) Lim HP (o) = exp [i(o — 7 log 2@)] (—i)+1 exp [io] 


+0 


(28) S. OkuBo and D. FELDMAN: Phys. Rev., 117, 292 (1960). 
() Fo), Gg) and H{ (0) go over to 2Î:(2), — 0n,(2) and oh) (0) for ZZ,e®—0 
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and contains therefore only outgoing waves for large o. The functions M 
and W used above are the regular and irregular confluent hypergeometric 
functions in the notation of WHITTAKER and WATSON (2). 

As can be easily seen from the asymptotic formulas (3.9) and (3.11), the 
wronskian of F, with HA is 


dA (0) 
- H® MET; Sl NET, Esp FS oh 
do 1 (0) F(0) do Ù 


(3.12) = 
Eq. (3.12) together with (3.11), enable us to write for the channel f, in the 
co-ordinate representation, the Coulomb Green’s function G,; satisfying an 
outgoing wave boundary condition: 


/ 


(8.18) ler) = #2 yet +P | 
L 


For VE Whe) Heidi 
27h? 


rr kr - ee 
with the usual meaning for r_ and r_. 
In terms of the confluent hypergeometric functions M and W, we have 


72 1 / My 2! re Teti nu 1++(2 ik, UE Aly 
(3.14) naro = 3 > (21)! r_ l (2 | (Qik, re BUTS 


21h? 
(= ik,)'I(1 "ls 41) Wei (— 2ik,;r-)]. 


The asymptotic expression for <r|G.,|r’> is obtained using (3.11): 


Von 1 
(3.15) him ic = za = exp [i(k,r — n, log 2kyr)]: 
Dx (por (E)lr') 


He 
where ky=k,r/r. g is the Coulomb wave function which satisfies an in- 


coming wave boundary condition with the outgoing wave packet in channel f: 


v k,-r\ F(k;r) 
(8.16) 99 (Er) = 6h27) 32 2 È ) exp [— set] Pe (SL a n 


p? is related to pg? through time reversal and space reflection: 


(3.17) le eles Sty v;) = (ir PUT les, Sty —Y}) » 


(27) E. T. Wuirraker and G. N. Warson: A Course of Modern Analysis (New 
York, 1948). 


en 
= 
+ 
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where the operators P and 7 are defined by (see ref. (28)) 


ee) r|\P=Crl|, T=) Pili Gok , 


ÿ 


where Æ is the complex conjugation operator and ol is the Pauli spin matrix 
for the n-th particle in channel j if it has spin 4, or io = 1 if it has spin zero. 
The sum on j extends over all the channels (*). 

The asymptotic behaviour of eq. (3.3) determines the collision properties 
of the problem. Besides (3.15), we need the asymptotic form of g® for which 
we refer for example to ref. (>). The result is 


(3.19) Lim atx? = exp [i(k,-r + n; log [kur — k;-r])]0, ds, + 
1 rA D, z Di 
- ED [t(kyr — 4, log 2k;r)] (2) OE 
f 


where, defining the T-matrix by 


2gt\? = 
(3.20) O,; = — (=) (Gus My ky ki Ts, 
we have 


(COLT {=o hood eS 


3fSi  VpVi FF (ge 


; Fun) = aT oes Ôô + (4 AN Kg?) À 


Ssssi VpVi 


T,; describes the scattering in channel + in the absence of everything but 
Coulomb forces. T,; is related, through (3.20), to the amplitude 


No Hi EXP [25061 _ ‘ U 
(3.22) DES ke) = 2h, sin? 48 exp [— in, log (sin? 40)] = 


= 55g, 2 (21 + Wlexp [oi] — 1) Ps(eos 9), 0 = kyk,. 


The square magnitude of @,, is the differential cross-section from channel à 
to channel f: 


d fi , î 
(3.23) Mao = |0;;(6,, 6.) |? 


(28) E. P. WIGNER: Group Theory, (New York and London, 1959), Ch. 26; G. C. Wick: 
Ann. Rev. Nucl. Sci., 8, 1 (1958). 

(*) For three or more body channels we need as many (io("?) quantities as the 
number of particles in the considered channel j, see ref. (17). 
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The matrix 7 given by (3.21) is simply related to the S-matrix defined by 
(3.24) Spe = (YP, pO) = (95°, Spa ) 

which shows manifestly the unitarity of the S-operator. wO is the scattering 


solution of the complete Schrédinger equation satisfying an incoming wave 
boundary condition with the outgoing wave packet in channel f; it is related 


to pw by an equation of the type (3.17). We have 
(3.254) Sn = Soe — 2ri BE — BGP, V'y®), 
where 

(3.25b) Son fai (oar. pee) 0719.5, ORO 


is the S-matrix for pure Coulomb forces. 
The unitarity of the S,-matrix furnishes the equation for 0,;: 


D7 


3.26 = 
(3.26) à 


[0.(k;, k,) — OÙ(k;, ky agent, k)OStER AE 


Use of (3.26) and of the unitarity of the full S-operator yields the following 
equation for the complete amplitude @0,,: 


"I I 
(3.27) > 2 Ot) 7 DE, E)| = DA LOTO RAC ry tas yo 
up Ua 


NSnVn 


where the sum over # runs only over the channels which are open at the con- 
sidered energy. If in (3.27) we put i=f and &,=&,, we get the optical 
theorem: 


~ 


(3.28) / È Im O67, E,) =“ > dQ, | QE: PIE = DES) Ù 


NSn Vn 
It is useful sometimes to use the S-matrix expanded in spherical harmonics: 


h20(E — E' i 
229)" Sy sa E o Di Vee) IMM E vr) 
Mi fp hi lp)" yh; 


Ca], M; M — vi; vi) Sfippitoa(P) » 


where C,;(J, M; m,v) and Y? are the Clebsch-Gordan coefficients and the 
spherical harmonics in the notation and with the phase convention of Blatt 
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and Weisskopf (2°). The unitarity of the S-operator implies that the matrix 


J DI e DARA 
5, re unitary. 


Use of (3.21), (3.22) and (3.25) yields the expansion of O, in terms of 
the S-matrix, in which the pure Coulomb scattering is separated out: 


(3.30) O,;( (E;, 4) a(k;, k; 000, ale 


Sx Do YIU(k) VEC) Ch}, My M — 95, 94)" 


k; i JMlpzl, 


>: Bee Od 
Cuald; UE M 5 Vis Vi) (exp [2407,] 0720p Ogg Sr) e. 


Using (3.30) and (3.21) we get the S-matrix in terms of the potential V': 


(3.01) Sue, tus, = CXD [27041010140 spa a (2 9) gUu=t exp [io7,]- 


> far arr By, (kV aye? (0) Wye, Pilsen 0") 


G1j8j 


where we have supposed V’ to be in general non-local (*). V'7(r' |r") and y’(r) 
are defined by 


92 9 È Rf U'* rl M TT 
(3:32) jure D'ns(T) Ver’ CRE = Ve | LE Honor 


(3.83) UE, r) = (nt D CJ, M; M— me, 0) VER) 


JMUsjli 
Pius (Kils 8s, 7) 
IGT) «E 
where 
(3.34) Yor > Cis(J ;M, y) Y IP Wer 
is a total angular momentum eigenfunction. 
yp’ (r) satisfies the differential equation: 
È 5 Ge dalla 2; e? 
(3.99) k : dr n he Ao =| pis, (Rs 168, r) = 
=> È qe Oe ie ores ane |) Prali, r'); 
0 


(29) J. M. BLATT and V. F. Weisskopr: Theoretical Nuclear Physics (New York, 1952), 
Appendix A. 
(*) In the local limit V'¥(r' |r") > (de 1") /r2) V(r’). 
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or, in turn, the integral equation: 


(3.36) mod r) = 4x exp [ici] Fi (kr) 0 500,0 


à 


x Ar! dr'r'rP,(kyr)BO(kyro) DEVI, mal? 0) Winn (lilies 11), 


kt,3x 1 h 2 


where r_ and r_ refer to > and r’. The asymptotic form of (3.36) is related, 
through (3.31), to the S-matrix: 


(13:37) Lim Whe, (KiliSi, 1) = 270441 Pyne exp [— i(k,r — n; log 2k;r)] — 


k;u;\? 
es (ea) exp [i(kjr — n; log 2k;r)]S; FR aes 

Use of time-reversal invariance and of (3.17) furnishes the principle of 
detailed balance: 


(3.38) Orleans) hee 


t 


where by —& we mean 
—g= (—k, 8, —»). 


Comparison of (3.30) and (3.38) shows that, under the hypothesis of time- 
reversal invariance, the matrix Sh, Bates is symmetric (*) (see also ref. (7) for 
three-body channels). 

We now proceed to the study of the properties of the cross-sections for 
scattering and reactions just above the threshold of the initial or final channel. 


4. — Low-energy behaviour of cross-sections. 


We will consider here the behaviour of the various elements of the tran- 
sition amplitude when either the final channel f is considered at its own thresh- 
old (endoergie case) or the initial channel i is considered in the limit k;=0 
(exoergic case). We will consider first the cases if, after that the proper- 
ties for the scattering case will be deduced as a consequence. We will consider, 
quite in general, the possibility that Coulomb forces are present in the various 


(*) We remind of the rules: 


Le ET de Toi ’ Cals M; Mm, y) =(—) 8-4 (1, (J, Mim, v). 
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channels; the case of no Coulomb interactions will be obtained by simply 
placing e?Z,Z,=0 in the final formulae. 

The following considerations will hold provided the potential V’ vanishes 
beyond a finite value of r. In most cases of interest this condition can be re- 
laxed though. Let us first consider the 


a) Endoergic case. For an endoergic reaction channel à — channel f 
(two-body), an expression for the transition amplitude 0,; in terms of the 
potential V' is obtained using (3.16) and the first version of (3.21). If we take 
the incident beam in the direction of the z-axis, we get 


CADE DELIO 


1 Mi Li de = 
CAD Palena Gai Z (2h + Vite’ exp [os]: 


bo 
"Cal, Vij da Dir Vr) Cig, (Fy Dis 07%)? 


(4.1) 
i Zi 
far dr'r'r Fy (kv) > 3, Vogt AAA Se, 10) 
il;s; n : ‘ 2 


In the limit as k, goes to zero, one can determine from (3.7) the low energy 
behaviour of the function F,. Since M,,,,13(2i0)/(2io)"! goes to a constant, 
the behaviour of OY” at k,— 0 is determined by the function (*) 


Ky exp [— tog] exp [—n a2], +1 + im. 


The limit gives 


kuti Ay = 0 4 

(4.2) Lim OÙ ~ V2a exp [— qa] Alti ”/1,! A, > 05 
(Pina 

Vox|A;,|##4(— i) 4/1! A7 08 


(endoergie case, f Zi). 

We see that in the case of Coulomb repulsion (A,> 0) in the final channel f, 
OF? starts from k,=0 like the exponential exp [— 24,/k,], while for Coulomb 
attraction (A,< 0) for all Z,, O? is finite there. From the low energy prop- 
erties of (3.8), we see that the next term in the expansion of OY for A,<0 
is proportional to k?. We have therefore that the derivative of the reaction 
cross-section o,; with respect to %, is finite at threshold (%,= 0). The case 
of no Coulomb forces yields the well known energy dependence. We now 
proceed to the discussion of the 


* pi | = = ig 4 Di 
(*) The term exp[—40,] has been introduced to avoid unobservable oscillations 
at zero energy. 
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b) Exoergic case. In order to obtain the behaviour of the transition 
amplitude ©,, when the initial channel 7 is considered at its own threshold 
(k; = 0), it is convenient to take into consideration the second version of (3.21). 
Following the same procedure as before, in the limit as k; goes to zero, we 
get (*). 


[hi ARA 

(4.3) Lim Of ~ | yak; exp | mos |AIM le Ar 0, 
ki —0 3 

3x k:+| À | hki(— j)u/l,! AS A0" 


(exoergic case, f# 1). 

The characteristic exponential exp{—xA,/k;] is still present for Coulomb 
repulsion. In the case of Coulomb attraction, O? goes to infinity in the 
limit k;>0. The energy dependence for A,= 0 is well known and yields the 
1/v law for the absorption of slow neutrons by nuclei”. 


c) Elastic scattering. If Coulomb forces are present in channel 7, then the 
combination of both (4.2) and (4.3) gives us the energy dependence of the non- 
Coulomb part of the scattered wave. Consideration of the pure Coulomb scat- 
tering amplitude 0,, shows that at very low energy only pure Coulomb scat- 
tering survives, since Lim Ou k;°. 

In the case of no Coulomb forces, writing 0,,(E,,&)= DOMME, €) YY"(k,) 
we get the well-known formula tale 


(4.4) Lim Quo LUE, A,=0 (scattering, f= 4), 
=> îi DSi 


which shows that the scattering cross-section is finite at zero energy when 
Coulomb forces are not present. 

Use of eq. (4.2) and (4.3), in conjunction with (3.30), gives us the threshold 
behaviour of the various elements of the S-matrix: 


(4.5) exp [240] dx da Oh ni plie a K 1,3 CHRIS Kits, 9 
where 
[ ket A, = 0: 
(4.6) fice | Var exp [— nl Ayrin|t,! MAY 205 
Va |A,|ttk(— i)! A; <0,; 


with a similar expression for Ky. Gr, i, IS à quantity which is finite 
(*) Eq. (4.3) can be clearly obtained from (4.2) using the principle of detailed 
balance (3.38). 
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and in general non-zero as either k, or k; approaches zero. The next term in 
the low energy expansions of .#?, for A,<0 (for A,<0) is proportional to 
k° (to ki). 

Tf there are no Coulomb forces in channel a and if we consider energies near 
the threshold of that channel, then it can be easily seen by the use of (4.5) 
and (3.30) that 

[tee a eee 

A | ai\ Say Si Org VAI ’ 

(4.7) È | 


al (ou (ROS = Oy a Va One Ore VA J 


i In Fig. 1 the energy behaviour 
a of the cross-sections (differential or 
A5=0 A;>0 Ap<0 integrated over angles) is plotted for 

34 the endoergic case (f 7?) vs. the mo- 

di DA mentum #, of the initial channel 7. 

ki_ We notice that in this case o,,; starts 

se NE hr with an infinite slope in the case of 

Fig. 1. — Energy behaviour of 0;; vs. the no Coulomb forces in the final chan- 


momentum %,;, in the endoergic case. 


nels} sli wesgdefimeWa == 2 EE 
= Y|0% |} for the reaction cross-section integrated over angles (*), remem- 


da 
bering that dk,/dk; = w,k;/u;k,, for the partial cross-section of! we get 


(Lp) 233 
(4.8) Lim 224. ~ (2,4 1) attra, | = 2? ele 
kr—0 dk, RE Ese 
(endoergic case, ft, A;= 0). 


So that we recognize that the /,=0 part of 6,; exhibits the infinite deriv- 
ative with respect to k; at k,=0. Notice that the derivative of o,, (differ- 
ential or integrated over angles, indifferently) with respect to k,, is finite at 
ka==<0: 

If there are Coulomb forces in channel f, then the reaction cross-section 6; 
(differential or integrated over angles, endoergic case) starts with zero slope 
in the repulsive case, with finite slope and non-zero value in the attractive case. 


(°) The same is true for the differential cross section; in that case we define 


ay) = A (0) (1) ¢ f (0) (1) 1 na Il 
10 pes oS) = 2 Re LOVE VE 
An if | f de 4/40 1 ( a 
(2) 


oi? 2 Re lemon 1 Yur * (ke l OO Yr-r(k,) |2, ete 
fi os pe fi 4/40 2 rf + | Vi fi ( A) |2 , ete. 
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d) Three-body final channel. We want to consider here the low-energy 
properties of an endoergic reaction cross-section leading to a final three-body 
channel. Let us consider first the case with no Coulomb forces in channel f. 
Formulae (3.19), (3.20) and (3.21) are easily generalized to cover this case. 
In particular we have for the transition amplitude @,, 


(4.9) DEEE = = (E) (pus ttrgliss]kE(PÎ (Ex) PE); Vp) - 


O,,, as given here, yields the cross-section for the reaction leading to the 
three-body channel f in which the first particle (n. 1) obtains a momentum 
between hk, and hk,+dhk,, and the particles 2 and 3 receive a relative mo- 
mentum #k,, in the direction between ,, and Ord Ge: 


do; 


. (D D = 
(230) dep 


19x:(€1, E033 6.) h 


gP (E) and go (6) are given as follows: 


4.11) | Po Hi toa ri) = (27%) 4% exp [ik,-r,], 
a one T3) = (277) Ÿy2s exp [ik23° res] ; 


where r, is the distance separating particle 1 from the center of mass of the 
system 243; Pos is the relative distance between 2 and 3. y, is the reduced 
mass in the incoming two-body channel è, us is the reduced mass of the system 
213 in the final channel f. 

We assume, as before, that the interaction hamiltonian V vanishes beyond 
finite values of 7, and rx. In the limit as the total energy E approaches the 
threshold for the channel f from above, we have 


ae Lim (ly la) = Mg, 


where we have again separated O,, into the contributions from the various 
values of J, and 1, orbital angular momentum of particle 1 and of particle 2 
relative to 3, respectively. It is clear that (4.12) holds also in the limit as 
either k, or k, approaches zero, the total energy Æ being kept fixed in the 
limiting process. 

The following property of 6, is established : 


19, 
(4.13) Lim <= = finite, 


Kan >0 ( 
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os 


6, is zero for k,,=0, for k,=0 and also for H= E,, we. right at the thresh- 


old for channel f (*). Remembering that 


Keke, ki h2ki 


2 Us 2 lh; 2 uy 


ni 7 
E; Ei; , 


where 
Ma(Mo + Mz) 
dy = DE dc: 
Mi Mik Mok Ms 


we get for the integrated cross-section: 


max 
Rae 


Lim | d?k, Or =) (is 
E> Ey 
0 


eae x 


d 

Lim 

SES) pony dle. 
0 


| d'k,6;; = (0%, 


jem x 


2 


Na = wh 
Lim —, | d*k,0,; = finite. 
bE, dk; 

0 


The second derivative of the integrated cross-section with respect to k,, 
is in general different from zero, and finite. 

If two of the particles present in channel f are charged, labelling them with 
the number 2 and 3, we find for @,,(1,, l3), as a function of k., in the limit 
k., — 0, the behaviour (4.2). For Coulomb repulsion we have that also the 
second derivative of the integrated cross-section with respect to k, is zero in 
the limit E > E, (*). For Coulomb attraction, while the integrated cross- 
section and its first derivative are still zero, the second derivative turns out 
to be infinite. 


5. — Threshold effects. 


We are interested now in the study of the energy behaviour of the various 
elements of the transition amplitude at the threshold for the channel a, a new 


() Notice that for 6,; only the part relative to 43=0 contributes to (4.13), waves 
with 4,3 >0 give zero contribution in the limit as ky, goes to zero. 

('*) This is true also in the case in which the third particle has a charge of the same 
sign of the other two. 
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channel which is opening up at the considered energy. In this section we will 
consider the case in which Coulomb forces are not present in channel a. 
Coulomb effects will be discussed in Section 6. The discussion will be led for 
a~i and af, since the other cases have been considered in Section 4. 

Our method starts with the realization that the complete Green’s func- 
tion G, defined by the equation 


(5.1) ; (E-H,—V)6=1 


and the outgoing wave boundary condition, must contain all the information 
we want to get on the energy behaviour of the various cross-sections at the 
threshold for channel a. In fact G appears in the transition matrix as follows: 


(5.2) Ire Loro de (in [Vi+tV'G Fan) ; 


SFSi Oven 
since a~i, af and T.,, has nothing to do with the channel a, all the 
threshold effects must be in G. 

To work this out we have to isolate the part in G which refers to channel 4 
only. This is most easily done by the introduction of two projection operators, 
P,and P,, the P, refers to the new channel a and has been defined in 
Section 2, and P, is given by P,= Dip where the sum extends over all the 


IF a 


other channels. This separation has been carried out in all details on p. 498 
of ref. (3°), and we give here the result: 


(5.3) Gee IEE C PSP. (MG EU), 


where G, and Y, are Green’s functions which satisfy the equations: 


(5.4) (HE —H,P,— P,VP;)G@ = re 


with the outgoing wave boundary condition. AE) is the effective hamilton- 


ian for channel a: 


(5.6) HE) = oa ae pe SE > VLC) Joo Vern ’ 
bb 


tilt 
£a 


which is clearly energy-dependent and, since at the considered energy some 
channels D’ are open, non-hermitian. In the following the term with the sum, 
at the right hand side of (5.6), will be indicated shortly as Va Gi. 


(30) L. Fonpa and R. G. Newton: Ann. Phys., 10, 490 (1960). 
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The various elements of G, depend only on the momenta in channels d'. 
Therefore when we take the derivative of G, with respect to k, for k,=0 
(threshold of channel a) we get zero. Near the threshold we have then 


= (Le GV) Poa a eal 


aw 04, 
ALA DIE 


where we have taken the derivative with respect to the modulus of k, so 
that eq. (5.7) holds both above and below the threshold for the channel a. 
We remind that k, goes into i|k,| below the threshold (see ref. (158)). This 
can be most easily understood from eq. (5.10) below (*). 

It is by now clear that all the threshold effects are in %, and they can 
be best exploited by considering the integral equation for %,: 


(5.8) Ga = Goa + Goal Vo + Var O5 Vo) Fa » 

where G,, satisfies the equation 

(5.9) (BH), = 1 

and the outgoing wave boundary condition. In the co-ordinate representation: 


exp [ik,|r—r’'|] above , 


(5.10) <r |Gya[r’) = — At Laas 
aah? |r—r'| lexp(—|kel|r—r'[] below. 


We take now the derivative with respect to |k,| of (5.8) near the thresh- 
old: 


0G, LG ae : oe 
Ti] = oa] (+ eet VaGaVia) Le) + Goa Vas + Va Goo) a 


(5.11) a 
) a] kal 


Eq. (5.11), together with (5.10), establishes the important fact that 09,/0|k,| 
contains only outgoing waves. If we multiply (5.11) on the left by G,} we 
get the equation 


TO 0G 
Dal == A, E E — {ESSE 
(512) [ (grey = Ga aa 


{1 wer i VavGy Va) Gat Ê 


(*) The rule holds true also when Coulomb forces are present in channel a. 
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(5.12) can be solved formally: 


0G, Fr 00 i 7 > 
(5.13) El ==> Gol Tk] {1 (Vaz | VG, aa) = 
a 0 rc a T T à T 
= {1 sa Ga(Vaa + VarGyVia)} ° Î {l + QE =F AXEL Fa) Sai È 


7 © [tal 
Since 0Z,/0|k,| contains only outgoing waves, no solutions of the homo- 
geneous equation (E- A(E))u = 0 can be added to (5.13). In fact the non- 
hermitian hamiltonian (E) contains a skew-hermitian term — irVa: 
-3(H — H,P,— P,VP.)V,a which prevents the presence of only outgoing waves 
in wu. 

We are left with the discussion of the derivative of Gé. Writing Go in 
the form 


| 
(5.14) 
| 


we obtain near the threshold: 


QG ia li 
(5.15) ne ARRE 


(5.15) tells us that the threshold effects come out only from the S-wave of 
channel a. (") is a symbol to indicate that above the threshold the upper 


quantity has to be used, below threshold the lower. 
Substitution of (5.15) into (5.13) and of the result into (5.7), after using 
the identity at the threshold 


(5.16) map. farai <r! [Pa = ¥ |9P> <9" | for ka = 0, 


gives the final equation for the transition amplitude at threshold (°): 


0 ) : 
(5.17) ait 0,(E;, E.) = (i i) > Ores, 126) Ox (Eas EO gn . 


Sa Va 


Only the S-wave of channel a contributes to (5.17). The Kronecker’s 0, 


comes from (4.7). 
The derivative of @, with respect to any other channel momentum is 


infinite both from above and from below, since 


(5.18) fires eles Vale (EN ("| a 


ka 0 dk, fo’ | ka 


\ 
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From (5.17) one can easily obtain the derivative of the differential cross- 


section 

" 0 do;(6; ,6:) , (Im | | 

(5.19) 3 ke] | 4Q, E i A 2041 (ona ) Oa (Ene A) O (Fe E.) 6 Va Vi ? 
and of the total cross-section defined by (3.28) 

h':90( Ò tot ae Siro Re O, ; D 

(b:20) 5 Kal TO; (&3) _ k, Im 2 Pei (Er E) O ai (ee) E,)0 Va Vi 


Eq. (5.17) holds true also for O(£,, É,) with the substitution of O,,(&,, Éa) 
with OŸ(E,, £.) on the right- heed side. The following equation for the /,-par- 
tial cross-section integrated over angles is then obtained: 


0 Im Î 
(5.21) se | OF Ses 2 ee SOY 1 E.) ONE sh) DICESSE Es) Over. 5 
C | Va] af SaVa 
From (5.17), by the use of (3.30), we obtain the energy behaviour of the various 
elements of the S-matrix (for fai) at the threshold for channel a. 


Writing in terms of the .#-matrix defined by (4.5), we have 


—~ 
Cl 
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bo 
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Algo ‘in (5.19), (5.20), (5.21), and (5.22) the derivative with respect to |k, 

is finite, while the derivative with respect to any other channel momentum 
(for example k; or k,) is infinite both 

Gi or œtot from above and from below the thresh- 
isa+f old (*). The measurement of the cross- 


sections (differential, integrated over 
LAS AA si angles, total), at the threshold of chan- 
nel a, will exibit a characteristic « cusp» 


k; or «rounded step » as shown in Fig. 2 


Fig. 2. — Various forms of cross-section The downward at a ee S-like 
behaviour at the threshold for chan- Step (3rd and 4th cases in Fig. 2) are 
nel a. not possible for the scattering cross- 


section when below the considered 
threshold there is only one channel (channel 1) open. To see this in detail 
we follow an argument given by NEWTON (5). 


(*) It is clear from our method that the reason why the energy derivative 1s infinite 
also below the threshold for channel a is a purely quantum mechanical one. No expla- 
nation on classical grounds can be found for it as for the same phenomenon from above 
which can be interpreted on the base of the conservation of total flux see (also ref. (8)). 
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Let us consider for simplicity the case in which the threshold effect appears 
in the /=0 wave of channel 1 and let us forget about the spins. We have 
for the {=0 partial cross-section at threshold 


(5.23) N 


where S,, refers to J=/,=/;=0 and ky is k, evaluated at %,=0. Below 
the threshold for channel a, S,, is unitary and can be expressed as exp [2i6,] 
with 6,, real number. Above the threshold the scattering phase shift 6,, be- 
comes complex, but since |S,,|?<1 from unitarity, its imaginary part must 
be non-negative. We have therefore near the threshold 


À i 

5 9 ) 

(5.24) Oy, = Osso + n | CAM 

where we have neglected terms in k°, k3, ete.; «>0 and dò is ò,, evaluated 


at k,=0. By the use of (5.24) we easily evaluate eq. (5.23) above and below 
the threshold for channel a. In the limit as 4, >0, we get 


x 
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(5.25b) Toul = Sn Do COS Ore below. 
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Eq. (5.25a) gives us the desired result that above the threshold the derivative 


of a is always negative. 


The ratio of (5.25a) to (5.25b) gives 


a 
(0) | 
bd 


TE (0) 
COR On 
0/0 | ka Ou 


= tg du0- 


Because of (5.18) we have that if tg O1 is positive then the threshold anom- 
aly is a cusp (1st case in Fig. 2), if tg dx is negative then the anomaly is 
a rounded step (2nd case in Fig. 2) (*). 

Ea. (5.25) enables us to obtain the scattering phase shift ò,, at the thresh- 
old for channel a by the measurement of the slope of 9% above and below the 
threshold in the experimental curve. 

Various other examples in which the direct use of eq. (5.19) and (5.20) 
furnishes information about scattering phase shifts, have been considered in 
detail by Baz (4) and NEWTON (°). 


(*) An equation similar to (5.26) holds true for any scattering cross section, at the 
in the perturbation approximation which assumes 


threshold of a generic channel a, 
the off-diagonal potentials to be weak; see equation (5.9) of reference (°). 
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It is now clear how a phenomenon of this kind can give information about 
relative parities and spins. In fact, only the {, — 0 part of the wave in chan- 
nel a contributes to (5.19), (5.20) and (5.21). Therefore, if it iS pos- 
sible to tell experimentally in which partial wave, in the initial or in the final 
channel, the cusp (or step) appears, the parity or spin of that channel will 
result from conservation of total angular momentum and parity. For example: 
if one discovers that in the cross-section o,, the cusp appears for 1,=0, then 
P,,=+1, P,, being the relative parity between channel a and channel j. Ja 
instead the cusp appears for l,=1, then P,,—=—1. This kind of analysis can 
be profitably applied, for example, to strange particle processes. The example 
just given works beautifully for è = xp, f= AK, a= XK, (see ref. (11:12,14)) (*). 
The experimental determination of the final partial wave in which the cusp 
appears for the reaction mp — AK at the threshold for X production, is able 
to determine the P,, relative parity. And, for P,;—-+1, only a spin } for 
ZX is compatible with a spin $ of A (spin of the K-meson taken = 0)(*). 

Threshold effects can be also useful for the determination of the reaction 
cross-section 0,, which in many instances is out of reach of direct experimental 
measurement. For the case in which one of the particles in channel a has 
spin zero so that s, is given, taking the square of (5.21) above and below we 
get at threshold (4) 


a) sull 
lo) | Ne lo) DI 
97 : Ù _ FR = 
(5.27) ag On | | + lag on) | = to (ho ka oot) , 
Oka |a d| kal lp 
a 3 4 : : 
si goes to infinity linearly in the 
limit k, +0, a k, has been placed to counterbalance this effect. Analogously 
Œa=0) : 7er inear a È o > pre) 
i; goes to zero linearly and %, makes it non-zero in the limit. Clearly 
only l,waves which are coupled to the 1,=0 wave satisfy (5.27). From 
agir q ‘ (te) ¢ 
measurements of the slope of 6,’ above and below the threshold and of 
—1 _(La=0) 3 EE: de 5 : 
(ko), one can determine the of cross-section at zero energy. If its 
right-hand side is known, (5.27) can be used to predict the size of the cusp (or 
rounded step) since near the threshold we can write: 0,; = 0;;|, ,+k,(006;:/0ka 
neglecting powers of k, greater than the first. 
Eq. (5.17) does not hold if the final channel is the new channel a. In fact, 
one can easily see that for f=a an extra term appears at the right hand 


side due to the derivative of gf in the transition amplitude. This term is 


where it is understood that »,=y;. Since o 


G 


Per 


(*) See also the recent paper by J. SucHER, G. A. Snow and T. B. Do Phys. 
Rev., 122, 1645 (1961). 

(**) Experimental work is now in progress: F. EISLER, P. FRANZINI, J. M. GAIL- 
LARD, A. GARFINKEL, J. KEREN, R. PLANO, A. PRODELL and M. Scuwartz: Rev. Mod. 
Phys., 33, 436 (1961); S. E. Wozr, N. Scumrirz, L. J. Luoyp, W. LASKar, F. 8. Craw- 
FORD jr., J. Burton, J. A. ANDERSON and G. ALEXANDER: Rev. Mod. Phys., 33 
439 (1961). po 
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proportional to the 1,=1 wave in channel a, so that it gives no contribution 
when consideration is given to cross-sections integrated over angles. The anal- 
ogous of eq. (5.21) then yields(*) 


A 


(6) 


ay, (ea'eet”) = — = (Hub) (keto), 
where »,=y;. (5.28) allows the determination of 6," from the slope and inter- 
cept of the curve (ko °) near the threshold. 

The same procedure applied to the case f =? =a gives no more than the 
known zero energy property: 

(5.29) NEI? 

Oka 

where the equality holds when below the threshold for channel a there is no 
channel open. For the example considered above of production of strange par- 
ticles in pion-nucleon collisions, by using eq. (5.27) one can determine, at zero 
energy, the cross-section ott %x for the process XK > AK which is not fea- 
sible experimentally. Analogously, use of (5.28) at the AK and XK thresholds 
yields the non directly measurable cross-sections oh and ok, respectively. 

Threshold anomalies can also give information about the existence of new 
particles. For example, the discovery of a narrow Singularity in the energy 
dependence of the cross-sections for the processes T+p>r+pandr+pti 
>r°{n would support the existence of the reaction T+pop°+n with 
production of the p° meson (**) (*). On the same line of thought, proton- 
proton collisions, and elastic neutron scattering by nuclei near the threshold 
for the (n,2n) reaction, give information about a « pion+nucleon » bound 
System (*) and about the dineutron (8), respectively. 

The threshold anomalies considered in this section have been observed 
experimentally in the elastic scattering *H(p, p)*H at the threshold for the 
reaction *H(p, n)*He in the form of a downward cusp (#), and in the scattering 
Li(p, p)Li at the threshold for Li(p, n)Be in the form of a cusp (5). 


(SAVE G. Zinov, A. D. Konin, 8. M. KORENCHENKO and B. PontEcoRvO: Zurn. 
Eksp. Teor. Fiz., 86, 1948 (1959); translation: Soviet Phys. JETP, 9, 1386 (1959). 

(82) V. I. GOLDANSKI and Ya. A. SMORODINSKI: Zurn. Eksp. Teor. Fiz., 36, 1950 
(1959); translation: Soviet Phys. JETP, 9, 1387 (1959). 

(*) The equation for the cross-sections which appears in (22). can be obtained 
from our (5.27). 

(3) Yu. D. PROKOSHKIN, V. I. RYKALIN and I. M. VASILEVSKI: Report D-678, 
Dubna (February 1961). 

(34) A. I. Baz and Ya. A. SMORODINSKI: Compt. Rend. du Congrès International 
de Physique Nucléaire (Paris, 1958), p. 579. 

(5) M. E. Ennis and A. HEMMENDINGER: Phys. Rev., 95, 772 (1954). 

(86) P. R. MALMBERG: Phys. Rev., 101, 114 (1956). 
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We end this section with a brief discussion of the case in which three par- 
ticles are present in the new channel a (e.g., photoproduction of two pions on 
hydrogen, production of an extra nucleon in nucleon-nucleus collisions, etc.). 

The guantity which exerts an influence on the scattering and reaction 
cross-sections o,; via removal of flux, is now the reaction cross-section Ga: 
integrated over the energy range available to one of the particles in channel a: 


wa 
pax 


| d*k;0,;. This quantity (see (4.14)) and its first derivative with respect to &, 
0 
are zero at threshold so that no threshold effect of the kind of a cusp (or rounded 


step) is expected in o,,. Since the second derivative is different from zero 
and finite, a discontinuity is instead expected in the second derivative of 6,; 
with respect to k; at threshold. Such expectations are confirmed by a quanti- 
tative analysis, on the same line as before, that we omit. The quantity that 
plays now the decisive role is the free three-particle Green’s function for 
channel a. The effect predicted is, however, not useful (*). 

Special attention must be paid instead to the case in which there are (at 
least) two different reactions with three particles in the final state. If one 
of the final particles («) is common to both the considered reactions, then a 
threshold effect is expected in the spectrum of the particle « coming from the 
lighter three-body channel due to (4.13) (with f=a). The effect appears as 
a function of the momentum of the particle « while the total energy is fixed 
and lends itself useful for the determination of relative parities and spins of 
the reaction products, and of cross-sections for processes not feasible experi- 
mentally. Details can be found in ref. (17-19) (**). 


6. — Coulomb effects. 


From Fig. 1 we see that when a repulsive Coulomb force is present in 
channel a, then o,; starts at the threshold like exp[—xA,/k,], while in the 
case of Coulomb attraction in a, 0,; starts with a finite value at the threshold 
and in general with non-zero and finite slope. We expect, therefore, that in 
the case of Coulomb repulsion the new channel a makes itself felt in the other 
cross-sections very smoothly and no cusp (or rounded step) will consequently 
appear in them (**). Also in the case of Coulomb attraction no cusp is expected, 
but instead another type of anomaly is to be obxerved. As the energy of the 


() The size of the discontinuity in the second derivative of o,; With respect to ky 
has been given by L. M. DELVES: n-d elastic and inelastic scattering near the inelastic 
threshold (preprint, M.I.T., 1961). 


(”) The method discussed. in this section is also applicable in this case with some 
changes. 
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incident particle approaches the threshold for channel a from below, we will 
observe in G,; a series of ever more rapid oscillations due to the physical pos- 
sibility of excitation of the infinitely many Coulomb bound states of channel a. 
The limit of o,; below the threshold will not exist, since these Coulomb bound 
states, and therefore the corresponding resonances in 6,;, have the threshold 
energy as an accumulation point. What will be seen experimentally is the 
average of o,; below the threshold of channel a, and this quantity will eventually 
exhibit a step-like discontinuity at threshold, to counterbalance the sudden 
leakage of flux to the new channel a (7). 

The Coulomb effects are able therefore, in both the repulsive and in the 
attractive case, to wash out the cusp (or rounded step), as given in Fig. 2, 
in the old cross-sections at the threshold for the new channel a. In most cases 
of interest, however, such as in processes in which pion and/or strange par- 
ticles take place, the Coulomb effects described above cover an energetic region 
around the threshold of the order of 10-* MeV, which is out of reach to the 
present techniques due to the poor energy resolution in the incoming beam. 
In these cases the overall energy behaviour of the various o,,; at the threshold 
for channel a will look like a cusp (or a rounded step) even though o,,; does 
not exibit infinite derivatives at either side of the threshold in either case 
(repulsive or attractive). In these cases therefore arguments like the ones 
given before on the determination of parity, spins, ete., will still be ap- 
plicable. 

To put into a quantitative form the above considerations, we start again 
with taking into account eq. (5.3) for the complete Green’s function. Since we 
are interested in a reaction for which both the initial and the final channels 
are different from the newly opened channel a, the threshold anomalies will 
be given exclusively by the behaviour of ©. We define for a function f(k,) 


(6.1) Af(la) = Lim, f(a) — (ka), 


with k, taken below and near the threshold for channel a. Then, for the 
complete Green’s function we have 


(6.2) AG = (14+ GV) PAG Pl + VG) 


since G, does not «see» channel a. 
In order to evaluate AY, we write for Y, the integral equation 


(6.3) Ya Si Gea + Glas r Vi, G, 1409 9, > 


Taking again into consideration the fact that we want only outgoing waves 


rd 
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into AY,, we get 


(6.4) AG= [1+ (AGca— Gea)(Waat Vis Go el AGeal1 + Vit Vas Go Moa) Sal = 
STi GAVE VG Ve) AG) tg ALP te ene |. 
where with the operator ie we mean 
(6.5) Tee > Ve. Onl ae 
ee 
which is so defined that 


(6.6) ATO CRU RS CS Poy 

All the quantities on the right-hand side of (6.4) are defined in the limit as 
the threshold energy for channel a is reached from above. Substituting (6.4) 
into (6.2) and using (5.3), one gets the following equation for 0,,: 


è 7 A re 27 \° (f)(—) 1 T' (à) 
(6.7) AO Crs E.) e (=) (i My k,/k; )} (oi ’ 153 [AG as Te ai Pe 


Obviously, there is no effect in 0,, since this quantity is diagonal on the channel 
index. We see, therefore, that the anomaly in ®,; comes exclusively from Gq. 

By using eq. (3.14), we can evaluate explicitly the jump experienced by 
the Coulomb Green’s function G,,. The quantity which is responsible for such 
a jump across the threshold is 


(6.8) (ie VIA) 
while at the threshold Merci x(20kar_)/(2ikar_)*! is a continuous function of 


the energy. The quantity (6.8 3) has been considered in great detail in ref. (7). 
From eq. (16) of that paper we have 


(6.9) ATK.) P +1 + ing) Won, 143 (— 2er) = 


0, fori A208 


== 2a | A,| 2242 M itng|.244( 24 | ka | 7s) 


en Dix ]ma]] (2l + 1)! Ghee ae? 


for A, Oe 


For the Coulomb Green’s function we get 


te for). Ape ADR 
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(6.10) Ar|Galr'> = phi  OXP [— 2ta0| Mal NT > (2756) P. (ee | 
Fi(lkalr<) F(|kalr>) | | 
ae lea lp n: i¥ | k, | rs, ; one AP == (H) 


492 


INELASTIC COLLISIONS AND THRESHOLD EFFECTS 143 


As expected, we do not get any effect in G,,, and therefore in ®,;, in the 
repulsive case. For attractive Coulomb forces in channel a we obtain instead 
a jump in G,,, and in O,,, when we go across the threshold. In the following 
we will therefore continue to discuss the attractive case. 

Substitution of (6.10) into (6.7) yields a simple expression for A0;; only 
if we expand in spherical harmonics: 


CINA era VEO) 
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k; JMlils 
LaSalaSa 
Y M - Î , V4 TXT. WV 7 a7 n 
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where S7 is the submatrix whose elements refer only to channel a. exp[210;] 
is a diagonal matrix in the /-representation for channel a. 
If we renormalize the definition of the S-matrix by 


(6.12) Sin = exp 107] Si, exp[— to] , 


which does not alter any of the observable quantities, and taking into account 
the fact that 
(6.13) Lim, exp [2i(o,— oo)] = exp [tal], 
we obtain an expression for AO, = exp [— ict] AO,, exp [— ioj] in which only 
S-matrices appear and the exp [2¢(o; — ziml)l is substituted with 
exp [ia(l—2|y.|)]; again, exp[ial] is a matrix in the l-representation for 
channel a. A@,,; is more suitable for applications in that it does not contain 
the unobservable oscillations at k,—0 which are present in (6.11) due to 

= Di, ~@ 
exp [2t0%]. x 

In terms of the S-matrices we get from (3.28) 


2 = 
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From (6.11) and (3.30) we obtain the energy behaviour of the elements of 
the S-matrix (for f4a#1) at the threshold for channel a: 


È = 1 == 
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The oscillatory denominator in (6.11), (6.14) and (6.15) represents the res- 
onances at each of the infinitely many Coulomb bound states of the channel a. 
These oscillations become infinitely fast near k,=0-, point of accumulation 
for the Coulomb bound states. Therefore only the average of the various ob- 
servable quantities will be detected experimentally. The averaging procedures 
for the differential reaction cross-sections and for the total cross-section have 
been given in ref. (9). The result is that, while the various averaged differ- 
ential cross-sections exhibit a jump across the threshold for channel a, the 
averaged total cross-section turns out to be continuous there. 

A phenomenon of this kind is for example expected to occur in the elastic 
scattering of X-rays by atoms at the threshold for the photoelectrie effect (*). 

There the step-like discontin- 


a-photoeffect 4 o-X-ray-elastic A ototor index of uity in the scattering Cross- 
refraction or mass ae F è 
absorption coefficient section is of the same order 


of magnitude as the cross- 
VAL /\NW section itself and is equal 
RA and opposite to the initial 
step in the o-photoeffect. 

frequency 4 : 
rea Fig. 3 shows the details of 
Fig. 3. — Photoeffect. the phenomenon. The effect 


has never been observed up 
to now since the experiments have been set up either to measure the index of 


refraction that, when measured directly, does not show the anomaly, or the 
mass absorption coefficient which does not exhibit the anomaly either. 

At the threshold for a three-body channel a, from the considerations at the 
end of Section 4, no effect is expected to occur in the case of Coulomb repulsion. 
If two of the particles in channel a are attractively charged, then an anomaly 
should be observed in the second derivative of the old cross-sections at the 
threshold for channel a. Coulomb effects of the kind discussed in this section 
are expected to occur also in the case mentioned at the end of Section 5. 

The method used in this paper for the discussion of the threshold anomalies 
applies in a straightforward way when Coulomb forces are present in the newly 
opened channel a. This is the reason why, for uniformity of treatment, it has 
been preferred over other methods which are able to reach the conclusions 
of Section 5 by following a different line of thought. 

Beside the well known #-matrix approach employed by WIGNER (1) and 
BREIT (2), essentially two other methods can be found in the literature: the 
one by BAz abd OKUX (1?) based on the unitarity on the Suns, Matrix used 
above and below the threshold for the new channel a, and the one by NEW- 
TON (°) based on a K-matrix formalism. These two methods are not, however, 
easily extensible to the Coulomb case. 
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The PAPA machine (1) can learn by itself without teacher, 7.e. the machine 
is able to organize its experimental data into classes without human inter- 
vention. Here is how it operates. 

The machine separates at random the N objects (examples) shown to it 
into two sets of N/2 objects—called class a and class b respectively—and learns 
these two classes in the usual way, as if they were really two objectively dif- 
ferent classes taught to it by a human teacher. After this, the same N objects 
are shown again to the machine for recognition. PAPA will then recognize 
N’ objects as belonging to class a, and N—N' objects as belonging to class b. 
Of course the N’ objects assigned to class a will not in general be those N/2 
objects belonging to the original class a, since PAPA does make mistakes with 
a finite number of A units. The reason for this is as follows. 

In the original class a, due to random fluctuations, there might be a higher 
percentage than in class b of objects satisfying a certain property. For example, 
let the objects of class a be on the average of greater size than those of class b. 
This introduce a bias (?) in that, other things being equal, there is a greater 
chance that PAPA will make mistakes when a large object of the original 
class b or a small object of class a is presented for recognition. It is now clear 
that the average size of the N' objects recognized as belonging to «class a will 
be even greater than it was in the original class. Similarly the N—N' objects, 


(1) A. GamBa: Proc. IRE, 49, 349 (1961); G. PALMIERI and A. Sanna: Methodos. 
to be published; A. GAMBA, L. GAMBERINI, G. PALMIERI and R. SANNA: Suppl. Nuovo 
Cimento, 20, 112 (1561). 


(*) A bias from our point of view: a perfectly valid criterion from the point of 
view of PAPA. 
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recognized as belonging to class b will now have smaller average sizes. We 
can now define two new classes a’ and b’ made up by the N’ and N—N’ objects 
recognized as a and b respectively, and start the process over again. A stage 
will be reached when the objects will be separated according to their sizes, 
or to the property that originally fluctuated the most. PAPA operates there- 
fore as a truly «intelligence amplifier ». 

The process can be applied separately on each of the two classes (large 
and small objects) in order to obtain a further classification into sub-classes, ete. 
It is clear that the complete ordering of objects obtained in this way is not 
necessarily the one a human would do, although it satisfies some of the more 
obvious requirements of a man-made classification. For example, the first 
dichotomie separation into two classes is based on «the most widely fluctuating 
property », which is also in most cases the most evident or simple property 
for the sense-organs of the machine. Therefore the machine starts from its 
« most simple and evident » properties and proceeds to the more sophisticated 
ones using a «criterion of simplicity ». 

Experiments on this type of learning are now in progress and will be 
reported later. 
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CoNTENTS. — 1. Introduction. — 2. The dynamics scheme relative to 
the Frenet frame. — 8. Some geometrical results. — 4. The equivalent two 
mass-point representation. — 5. Introduction of internal co-ordinates 
and their conjugate momenta. — 6. Concluding remarks. 


1. — Introduction. 


In a previous paper (1) (henceforward referred to as I) the motion of a 
finite relativistic free particle with spin was discussed in terms of the co-ordi- 
nates æ, («=1, 2, 3, 4) of some point P within the body and the direction 
cosines U". (u=0, 1, 2, 3) of a tetrad of orthogonal vectors forming a frame 
of reference moving with the particle. A generalized Hamiltonian dynamics 
was proposed for the particle by defining Poisson brackets (P.b.), the Hamil- 
tonian, and the equations of motion. It was shown that the dynamics is 
invariant to an orthogonal transformation of the U"“.; and that in particular 
for the U", could be taken quantities V“ which satisfy the Frenet equations 
in four dimensions. The frame of reference, whose vectors have for direction 
cosines the quantities V" will be called the Frenet frame. 

It was postulated that the free particle had a conserved momentum p, 
and total angular momentum M,;, the latter being the sum of the external 


() Now at the Department of Applied Mathematics, University College of North 
Wales, Bangor. 


(*) J. B. Hucnes: Suppl. Nuovo Cimento, 20, 89 (1961). 
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angular momentum (or moment of momentum) L,, and the spin angular 
momentum S,,. The momentum was defined as a linear function of the U“; 
the spin angular momentum, a quadratic function of the VU“. By requiring 
that certain P.b. relations should be satisfied (*) the coefficients appearing in 
the definition of p, and S,; were uniquely defined in terms of quantities Li, 
which appeared in the equations of motion for the U".. It was assumed that, 
for a free particle, these coefficients are independent of position along the world 
linetof P. 


2. — The dynamics scheme relative to the Frenet frame. 


The dynamical variables are x, and V“, where V°,, V*,, V°,, V* are the 
direction-cosines of the tangent, normal, binormal, and trinormal, respectively, 
of the world line of P. Let & be the arc length measured along the world line 
of P, and for any quantity u = u(x,, V") let w=du/dé. 


DARÀ 


Then 
(2.1) Vs pae 1 
(2.2) ZIA AS 


where the quantities 2, are given by the matrix 


1 
0 = 0 0 
© 
| 
DE. 0 = 0 
(2.3) Q=| © 0 ; 
0 — 0 0 
(0) 


and 1/o is the curvature, 1/o the torsion of the world line of P, 9 and o being 
constant. From (2.2) and (2.3) we see that the world line of P has zero tilt. 
The Hamiltonian is defined by 


m 
Tu TL 
(2.4) H= GI J al a ° 
Poisson brackets by 
ou dv Ou Cv du Ov 
2.5 mu, v) = — — — + Quasar aor) 
et (ra) dn OV OV, OX. OV", OV'« 


() They were that (a, ©); (Px» Pe); (dx Po); (Mag, Mu») should have the same 
values as in classical point mechanics. 
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— where x and v are any functions of the dynamical variables—and the general 
equation of motion is 


From I (eq. (5.1), (5.3), (6.23), (6.24) and (9.8)) the momentum p, and spin 
angular momentum S,, are given by 


Oh 
(2-7) Da = mV, +m—Ve; 
Q 
(2.8) Sap = mo(V'nV *2— Vis) 


and the spin pseudovector s, is given by 


x 
9 E TI aya) aya 73 
(2.9) $= MOE pV aV pV y = EmoV",, 


where e is a pseudoscalar of magnitude unity. 
The equation of motion for the point P (I (7.2)) becomes 


> : fOr, AL ff il ae 
(2.10) D, = A, cos (2€) + B, sin (AE) 4 P | ) = 
NUO? 0? 
= A, cos (AE) + B, sin AE + mle | È 
m \o* 4-02 

where 
(2.11) (ae re 

010" 


Clearly it is the sum of an uniform rectilinear motion and an oscillatory motion. 


3. — Some geometrical results (*). 


The Frenet formulae show that the world line of P is a curve whose curvature 
and torsion are constant and whose tilt is zero. Since the tilt is zero it is a 
3-dimensional curve and lies in a flat, called the oscillating flat, to which the 
trinormal V? is perpendicular ($$ 128, 139, 146, 166). Hence the spin pseudo- 
vector is orthogonal to the flat in which the world line of P lies. 


() These come from the indicated paragraphs of A. R. Forsytn: Geometry of Four 
Dimensions, (Cambridge, 1930). ; 
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The world line of P, being a curve with constant curvature and torsion, 
and zero tilt, is a geodesic in a sphero-cylindrical region and is the extension 
to four-space of the cylindrical helix in three-space (§ 302, Ex. 2). 

When there is no tilt, the infinitesimal transformation connecting the V" 
at consecutive points on the world line of P can be effected, in the flat to 
which V3 is perpendicular, by a rotation round the line V*, = (0/0)V®, in the 
plane V*,=0 ($ 166). 

The quantity À appearing in eq. (2.10), and which gives the frequency of 
the oscillatory part of the motion of P, satisfies the equation 


iG el _ de 
Gide 


where . is the inclination of two consecutive vectors in the direction of V*, 
(§ 131). 


4. — The equivalent two mass-point representation. 


Eq. (2.2) with u=1 gives 


= 1 r i! 79 
(4.1) VT 4e Le Vas 
O o 
so that 
Ve "i cl nt VO, 
Q 


Therefore (2.7) can be written 


= mV, +m [Ve + eV al 
0° 
da o? d 
= EI ie = 02%, 
m az * m of dé (Da + ); 
or 
> da, i dy, 
y = = m : 
(4.2) (Oe = dé dé 
where 
(4.3) You = Wy, + 0% Ly, = Ly + OV x 5 
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and 
o? 
(4.4) MEO E 
02 


We will call the point, whose co-ordinates are y,, Q. 

Expression (4.3) for p,, the momentum, could be regarded as the sum of 
the momenta of two point particles, one of mass m whose co-ordinates are x, 
and the other mass m' whose co-ordinates are y,. 

Consider now the expression 


, dy 
dé 


oder ea BERE, ROLO cura lac 
= p°,m Al 34 al Vv“ o V è) om |! val ol n po = V o 


O i 


dé 


(Yu— 2,)m — (ys — ©g)m 


leg 


= mo Va Vs = VV) == Sop a 


Thus the expression for the spin angular momentum is equal to the expression 
for the moment about x, of the momentum of the mass m' whose co-ordinates 
are y. 

Moreover the point y, is the centre of curvature, corresponding to the 
point «,, of the world line of P (see FORSYTH, § 130). 


5. — Introduction of internal co-ordinates and their conjugate momenta. 


Let us now transform the quantities V“, to quantities W given by 


dl Il 
D caz ni er 0 
Wa cou 1 eri 
Wipe 
(5.1) 
il ni 
Wi, persa 
€ oÀ (o Î ah y ay 


Ven 
It is seen immediately that the vectors W“ are orthogonal 
(5.2) | DEEE 
and related -to the V# through an orthogonal transformation 


(5.8) We. = 0,,V",. 
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COUPE. rat” fr Le el La MEL PTE Pike My Su eh Zu r 7 RI} 17 Sort QUE yes pe va vert NY % 1 + 
pe er i ti GaN pe x i: eh rd ti. Le fi , A ah i 4 2 { Le mi iv 
SU AE LS RIT ir È “© nd a oh Ale pi je 
a a x a Pi ni x PS F = È 
BE ‘ ) ; vo 4 woe Arne 
Ger’ A GENERALIZED HAMILTONIAN DYNAMICS ETC. - II 153 ee 
The essential form of the dynamics is unchanged by this transformation ie 
but note that è, and W° are not equal. È 
The derivatives of the W”, with respect to & are sita 

170 ° 

(5.44) W”, 0, a 
4 
9 rr 2 x i x | 
(5.4) Wi= Aw, à 
(5.46) w= AW, È 
_(5.4d) We re 0. ER 
or ayy 
ve — O' wr ie 
W ra" QW a? is 
where ; pi 
0 0 0 0 | 
0 0 A 0 | 
i 
Quy =“ ’ 
0 —A 0 0 
| 0 0 0 0 
W°, is a unit vector parallel to the total momentum p,, 
W*, is a unit vector parallel to PQ, 
W*, is a unit vector parallel to the velocity of the point Q relative . 
to the point P, 
W?, is a unit vector parallel to the constant spin pseudovector. it 
; tn 
Using (4.3) and the result h è 
De 


ou 


0 
0 ou 


Ò Ve 


the Hamiltonian H is given by 


OW", 


(5.5) H=T wwe, 


» 19 SA 9 a 
where 
(5.7) OO nO; 
ñ 8, 11 - Supplemento al Nuovo Cimento, 
i - cs 


De 


du Ov 
ow", OW’ 


(I(10.2)), à eee 
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Thus 
1 [ou dv dv Ou ope du ov dv du | 
n° TERI 1 E 2 z r I 
a ee) os = lee dW. dn. dW*,| | 0A \00, oW ox, OWS 
: ou dv ou ov | 
AGE wi awe OWA) 
Then 

70 i 
(5.9) miss We) = a Ong > 
(5.9b) m(x,, W 7) 0 

sa ee 
(5.9c) Mita, W's) = = ÒxB 9 
(5.9d) MW W)=0, 
(5.9¢) m(W",, W°s) = 0, 
(5.9f) m(W*, Wi) = Xò, - 
If we define the point 2, by 

1 

D Ba Ba + RM 
(5.10) a Ba + Ne 


we find the following set of P.b. relations: 


(5.114) mes, Wg) == = Ox > 
(5.115) me, W,)=0, 
(5.11c) m(z,, W*,) = 0, 
(5.114) miW°,, W?,) =0, 
(5.11e) MW We.) = 0, 
(5.11f) MW W*,) = 6,8. 


We have already seen that p, can be considered as the sum of the mo- 
menta of a particle of mass m whose co-ordinates are x,, and a particle of mass 
m'=m(0*/0*) whose co-ordinates are y, =x,+oW',.. It is easily seen that 2, 
is the weighted mean of #, and y, with weights m and m’. 

Moreover 


: i Il lo) a 
5.12 = II AD, WI 
( ) ) ae on? a À: 02 L a 0 J ho W x 9 
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and hence è, is constant. We can regard 2, and moAW® =p, as conjugate 
co-ordinates and momenta. 

The quantities Wi appear as co-ordinates. We could regard them as in- 
ternal co-ordinates for the particle (with +, the external co-ordinates) in a 
classical bilocal theory corresponding to that introduced by YuKawa (2). Then 
eq. (5.4b) and (5.4c), which can be written as 


| Wise À 0H" 
= | ns OW 
(5.13) Ss 3 on 
|" “in OW,’ 


can be regarded as the Hamiltonian equation for the internal motion and 
suggest that we regard W? as the momentum conjugate to the internal co 
ordinates W!.. 

It will be seen from the (P.b.) relations (5.11) that the (P.b.) relation be- 
tween an «external variable » and an «internal variable » is zero. 


6. — Concluding remarks. 


The formalism developed in Section 5 suggests how the dynamics could 
be quantized. 

Suppose that, in a quantum theory, the particle is described by a wave 
function y= y(«,, W£,) satisfying a differential equation 


(6.1) (H —h)p =0, 


where H is the operator corresponding to the Hamiltonian, and h an eigen- 
value of H. The (P.b.) relations (5.11) suggest that the operator form for H 
is to be obtained by replacing W°, and W*, in the expression (5.5) for H, by 
operators (—1/moA)0/dz, and (—2/m)0/0W*,, respectively. 

Because of (5.10) we can consider p=y(z,, WF“); and if we assume that 
a solution of (6.1) exists in the form 


(6.2) y = v(z,) 0(W",) 


eq. (6.1) reduces to the pair of equations: 


02 
(6.34) | eee n) p=0, 
0? a ; 
(6.30) (RESA + WLW, + WWE, +e n) 0=0. 


(2) H. Yukawa: Phys. Rev., 77, 219 (1950). 
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The former is an equation of the Klein-Gordon type but with the mass 
term x determined by the solution of the eq. (6.3b) for the internal motion. 
This suggests that a mass spectrum will be obtained corresponding to various 
solutions of the equation for the internal motion. Also from (6.3b) there will 
appear extra quantum numbers in terms of which a possible particle classi- — 
fication could be obtained. 

The quantization outlined above will be discussed in detail in a further 
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Introduction. 


The success met with by the ordinary Schrédinger equation in the study 
of nuclear phenomena at not too high energies supports the reasonable view 
that the interactions between nucleons may be described, at least within certain 
limits, by an ordinary potential. Besides this we have few pieces of information 
on the nuclear interactions and these essentially of a qualitative character 
only; we know that nuclear forces are short-range forces and have saturation 
properties; further we know that they depend on spin, on spin-orbit coupl- 
ings, both static and velocity-dependent, and, perhaps somewhat less cer- 
tainly, that they are not dependent on the charge state, once the exclusion 


principle is taken into account (1). 


(1) Books which may be consulted with profit in connection with the nuclear torce 


problem are the following: J. M. Bart and V. F. WEISSKOPP: Theoretical Nuclear 


Physics (New York, 1952); R. G. SACHS: Nuclear Theory (New York, 1953); A. 8. 


Davrpov: Teoriya atomnogo yadra (Moscow, 1958). 
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These qualitative facts together with general principles which we reason- 
ably assume to hold in nuclear interactions (2), allow us to write a general ex- 
pression for the nucleon-nucleon potential. (We shall restrict ourselves to this 
case, in the hope of losing no essential feature of the nuclear interaction.) 

The general principles which we have mentioned are the following: 


1) hermiticity V= Vi; 


2) translational and Galilean invariance which require, respectively, that 
the potential is a function of the relative co-ordinate rand the relative mo- 
mentum p 


V= Vio, oi, r, P) s 


where 


o, o® are the spin operators for nucleon 1 and 2, respectively; 


3) rotation and space reflection invariance, which requires the « scalar » 
nature of V; 


4) particle exchange invariance, which requires 
T (1) (2) = 7.4 2 . 
J (o ,0 » TP) = V(o®, CRT: —p) ’ 


5) time reversal invariance which requires 


(2) 


Pos O ,r,; P) FAO lio o r, CD) ’ 


[(A, 4, ODO A IA DI = AAG, see A, Aj] ; 


6) charge-space rotation invariance (charge independence). 


The quantities involved in the problem are the vectors r, p, o, of, {™, 
z®, where 7%, 7 are the isotopic spin operators of nucleon 1 and 2, respec- 
tively (*). 

Then the most general expression for the potential, if we allow a quadratic 
dependence on p at the most, depending on the tensor invariants which 
may be built up with the vectors above and satisfying the aforesaid condi- 


(*) L. ErseNBuD and E. P. WiGNER: Proc. Nat. Acad. Sci. Wash., 27, 281 (1941). 
(*) We shall denote by heavy type the vectors of ordinary space and by a lower 
arrow the vectors of charge space. } 
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tions, is (3) 
V = Vo(r, p) + Vilr, p)(o®-0) + V(r, p) Si. + Valr, PUL:S)+ 
+ Valr, PR (6% L)(o-L) + (6 L)(ePL)] + V(r, p)(o-p)(o®: p) + 
+ {V.(r, p) + Vilr, p)(o-o®) + Vi(r, p)Si, + V(r, p)(L, S) + 
+Violr, p) £[(6%: L)(o®- L) + (6 L)(6®- L)J+Vis(7, p)(c-p)(c®-p)}(t®-12), 
3(6-r)(6 +r) il 


S (oo); i — — (rXp), S = 1(6V4+ o). 


Di h 


y? 


The functions V(r, p) are real functions of the moduli of r, p and this is 
all that can be said of these functions on the basis of the aforesaid principles; 
the information of a general character which we mentioned at the beginning 
may give us certain qualitative indications (for instance they must practically 
vanish at a distance of the order 10-!* em, the relative weight they must have 
in order to explain the data of the bound state and the scattering states, etc.) 
but no quantitative or detailed indication. 

Three approaches are possible to the search of a form of these functions: 
a purely phenomenologic approach, a purely theoretical one and a mixed one. 
In the first case a potential is built up phenomenologically on the basis of the 
analysis of the experimental data: one example of a potential of this kind is 
the potential of Gammel, Christian and Thaler (*); in the second case the 
problem of the two-nucleon interaction is treated according to the meson field 
theory, from which an expression for the potential between two nucleons is 
derived; in the third case a kind of fusion between the two approaches is done 
by taking what of «certain » is given by meson theories and by supplementing 
it by phenomenological information whenever the meson theories give un- 
reliable results: one example of this kind is the Signell-Marshak potential (5). 

In the past few years the two-nucleon problem has been the object of ex- 
tensive study and widespread interest (°); we think therefore that a critical 
account of the results reached up to now in the framework of meson field 
theories (*) might be of use (’). 


n J. L. GammEL, R.S. Cristian and R. M. THALER: Phys. Rev., 10551 (L90708 
(5) P. S. SIGNELL and R. E. MARSHAK: Phys. Rev. 106, 832 (1957); 109, 1229 

i A detailed survey of the phenomenological aspect of the problem may be found 
in the article by L. HuLrHÉN and M. SUGAWARA in: Handbuch der Physik, vol. XXXIX 
(Berlin, 1957). More recent results can be found in: Nuclear Forces and the Few-Nucleon 
Problem, vol. I (London, 1960). 

(*) Within that particular approach called «weak coupling approach ». 

(7) Review articles of wider subjects in which this problem is treated are; R. J. N. 
Puituips: Rep. Prog. Phys. 22 (1959); R. H. DaLitz: Prog. Nuel. Phys., 4 (1955). 
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1. — Conditions and limits of a theoretical potential; formulation of the problem. 


1'1. — Before formulating the problem of the construction of the potential, 
in the framework of meson theories, it is expedient to point out certain con- 
ditions which must be met by the potential and which are essential to the 
concept of potential itself as a description, valid only within certain limits, of 
the complete interaction between the two-nucleons (7): 


1) it must be hermitian; 


2) it must be defined only for energies of the nucleons lower than the 
threshold for meson production; 


3) it must be independent of the total energy W of the system; 


4) it must be such that the solutions of the Schrédinger equation give 
for the binding energy of the two bound nucleons and for the phase shifts in 
the scattering processes the same results as those derived directly from field 
theory (*). 


These conditions impose on the potential operator restrictions both in the 
r dependence and the p dependence. 
The condition 2) may be written, in the center-of-mass system, 


DIE 
ZE 


M 


(M is the nucleon mass and yw, is the 7-meson mass). 

Considering p as the nucleon recoil momentum and applying the uncer- 
tainty relation, we obtain that a theoretical definition for the potential shall 
be possible only for nucleon separations 7> un with 7. & 0.55-10- em, 
i.e. for a separation larger than about one third of the meson Compton wave- 
length. (We shall see later other reasons supporting this limitation.) 

The region r<y7,,,, therefore must be excluded from a theoretical treatment 
and, if one wishes to retain the description of the interaction within this region 


min 


(§) A. KLEIN and B. H. Mc Cormick: Phys. Rev., 104, 1747 (1956). 

(*) K. NisHisima: Suppl. Progr. Theor. Phys., 8, 138 (1956). 

(9) M. TAKETANI, 8. NAKAMURA and M. SASAKI: Progr. Theor. Phys., 6, 581 (1951) 
(reprinted in Suppl. Progr. Theor. Phys., 3, 169 (1956). 

(") Phase shifts have in the framework of field theory a well-defined meaning. 
Bound states may be introduced in the following way. The variety of the solutions 
of the field equation contains a subvariety of solutions corresponding to a value of 
the mass of the system of two particles which is less than the sum of the masses of 
the two particles: it is assumed that this subvariety describes bound states. 


eee See US, ee Cee ee, of > ÉD a A A ee PRI we a - 1 
rt Ve MG va (E. TE A i bi Y 


THE NUCLEON-NUCLEON POTENTIAL'IN QUANTUM FIELD THEORY - I 161 


by means of a potential, this can be done only on a phenomenological basis. 
For this purpose the assumption of a hard core is at present the most common 
attitude. 

As for the p dependence, the conception of a theoretical potential, as has 
been outlined above, is valid and useful for small relative nucleon momenta. 
This implies that the retardation effects in the transmission of the interactions 
between the two nucleons are small, thereby suggesting the inclusion of quad- 
ratic terms in p at most. 


1°2. — In the framework of field theory the problem may be formulated in 
the following way. The field equation 


(1.1) ih = |P(t)> = H|B()>, 


or the corresponding stationary state equation 
(4715) (H—W)|®>=0 
must be reduced, in an appropriate variety of state vectors (see note p. 110), 


to a Schrédinger equation which should, within certain limits, be physically 
equivalent to each of them respectively: 


a p? 2 
2 = |o(t) => |= + lœ(t)> , 
(1.2) nz Loto) (+) 
Pp e : + 
ai SEAN] > = Wil 
(1.2") (E TL Jie 1Œ 


Here the ket |®» is the field functional, H the total field Hamiltonian, sum 
of the free field Hamiltonian H, (H, = H,+H,,) and of the interaction Hamil- 
tonian H' 


H= H;XH&'H,;4H,+H', 


W is the total energy of the system and the ket |g) is the ordinary Schrodinger 
ket describing the state of two physical nucleons. 


1°3. — The question may be posed if a potential meeting all the requirements 
stated above can actually be found. The problem may be conveniently dis- 
cussed in terms of S-matrix quantities and, as it will be seen from the following 
discussion, the answer will turn out to be, within certain limits, positive (11). 


(11) M. A. Braun: Soviet Phys. J.E.T.P., 10, 582 (1960). 
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Let us treat the two-nucleon system in ordinary wave mechanics and let us 

denote by |g,) the eigenvectors of the unperturbed Hamiltonian h, by |pS” 

the eigenvectors of the total Hamiltonian A,+V with outgoing waves, by 

lp,> the eigenvectors of the total Hamiltonian corresponding to bound states 

(it will be W,,< W, when the total momenta of these states are the same). 
An element of the t-matrix is defined by the relation 


(Hr 


(1.3) top = SPa|V lve > 
from which it follows that 

Psv") = dg — tapi (LH, — Hg — ™) . 
The link with the s-matrix is given by the relation 
(1.4) Sp = Og — 2700 (E, — Ete - 


Matrix elements #,,, may be defined in an analogous way by the relation 


<q a LR = Rates — E,) à 


The matrix elements ¢,, are connected to the matrix elements of the potential 
operator by the well-known Lippmann-Schwinger equation: 


Votre 
1.5 EI 7 En 
à a 2 pr Hg — in 


This may be solved with respect to the potential, making use of the com- 
pleteness of the set |p; >, |y,,> (1°): 


+ 
t x mere 


LB 
ge t i DS AA LUE I 
VE o 75 | = 5 Î . 
25 E i 2 i 


m 


(1.6) V 


The conditions of orthonormality and completeness of the sets |wS>, |pm> 
. x . . : . . 
on one side and |g;) on the other impose on the t,; and t,,, certain restrictions 
which are typified by the following equation 


RIE pr bant55 
(1.7) e een 
3 2 (i, — By ip Bi 0) sd 


(2) Yu. V. NovozniLov: Soviet Phys. J.E.T.P., 8, 515 (1959). 
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where the summations are over the spins and momenta of the two particles 
and their bound states. 

These conditions are very important and assure that (1.5) with the potential 
V will have as its solution the same ¢,, as appear in (1.6), and that V is 
Hermitian. 

In quantum field theory the elements of the 7-matrix are defined by a 
relation analogous to (1.3): 


(1.8) MAI EIA 


where |®,> are asymptotically stationary states, which are the analogous of 
the states of non-interacting particles in ordinary wave mechanics (*), |?) 
eigenstates with outgoing waves, which develop from the asymptotically sta- 
tionary states |®,> with energy Æ,, H is the total field Hamiltonian. The 
S-matrix is expressed in terms of the 7,, by a formula analogous to (1.4). 
Beside the 7,, we may define the analogous matrix elements 

Tru = (PH — 2; Pin ; 

where |¥,,> are bound states in field theory. 

It might be thought to assume the field-theoretical T,, and T,,, as wave- 
mechanical transition amplitudes #,; and #,,, to be inserted into eqs. (1.3) 
to (1.6); the potential would then be defined by (1.6). 

This identification is prevented by the following argument. 

The quantities ¢,, and t,,, which one inserts in (1.6) must satisfy conditions 
of the type (1.7); now the quantities T,, and T,,,, do satisfy equations of the 
type (1.7) but with summation taken over all the field-theoretical states, including 
therefore, as well as states with T=y and M =m, (corresponding to the two 
particles we are treating), also states with /"4y and Mm (corresponding 
to other sets of particles). 

The relations satisfied by the T,, and 7°_,, are therefore different from those 
prescribed for the #,; and f,,,. ers 

The attempted identification may be reached for other quantities T,,, Tam 
which we will define presently. Let us separate from all the states |[#,7> the 
states with [= y: let P be the operator for projection onto these states, so 
that 
Qu) PPP") = SAY CROLEY + DPL. 


“1 n 
/ 


(*) These asymptotically stationary states may be thought of as those introduced 
by Van Hove (13) or by NovozuiLov (?); in most common S-matrix considerations 
they are often assumed to be eigenstates of the free Hamiltonian H). 


(3) L. Van Hove: Physica, 21, 901 (1955); 22, 343 (1956). 
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We define the state vectors 
(1.10) CODE Am 
7 

where N is given by 
(105100) CB IPB, = bp + Nap: 
The states | PD) are orthonormalized 

(OPODO 0: 
and are assumed to be a complete set in the subspace determined by the 


projection operator P. 
If we define matrix elements 7',, and 7°,,, through the following relations 


an 


=> Sl 
12 CEI DELA Rene) 
(il ) D, | B ap FREE 
b i 
(1.13) CAR EEE 


the a and Da satisfy, in virtue of the orthonormality and completeness of 
the states |», |Y,,> on one hand and of the | P®.} on the other, the pre- 
scribed conditions (1.7). Then, as far as conditions (1.7) are concerned, nothing 
prevents us from inserting the Do in place of the ¢,, into the definition of 


‘the potential. 


All this is, obviously, up to now purely formal; but it can be shown that 
the Ts, are identical with the 7,, on the energy shell (H,=,), so far as 
processes of creation and annihilation are impossible. 

Indeed from (1.11) we have 


Ly 


This expression has no singularities for E,= E; if processes of creation and 

annihilation do not occur; therefore, as can be seen from (1.10), the singularities 

of LOA EO’ and <®,|[W°°> are the same under these circumstances. 
“Then 


lim (E,— #,)®,|P>= lim (H,— B®, |P , 


Ey—EB Ex>Eg 


which indeed prooves the assertion we have made. 
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The «potential » is so defined by the following relation: 


~ Tors Pavia 
LL) SV N E 
| i 2 o, Ein 25, E; 


SOT == HPO SE > (O,|H — B,| VP) POW _ POST 


en > Tat 
SN ee ee a PO,’ ae ye 
ne y RENTE v7] 


VA, Pure NY dar 
> 2 ea PH — E, PD 
m 1m B 
We see so that, for energies lower than the creation threshold, it is possible 
to define, in the framework of field theory, a «potential» which correctly 
predicts the field theoretical cross-sections. 
However, two important remarks must be added to the discussion out- 
lined: the first regards the uniqueness of the definition, the second the non- 
local character of the «potential ». The «potential » defined above is not unique: 


the discussion is unchanged if in place of the set @, another set is used: 


I A | gum TT 
(1.10’) DID 
B 
2 5 3 3 | ; Fe 
where U is an arbitrary unitary operator, which can be written in the form + 
U=1+0,, with U, a non-diagonal matrix which goes to zero when N —0 
and which has no singularities at £,= £; when N_; has no such singularities; ; 


in this way one obtains a « potential » 


(1.14’) oa Di KO, FR ah eO TRE 


V1Y3 


which is not a unitary transformation of V and which may describe scattering | 
processes and bound states equally well as the original V. 

As for the second feature it is easy to see that the « potential » which we 
have defined is a non-local one; the equation equivalent to (1.1’) reads 


1 i di V4 ! 
(1.15) (2E(p) —W) (p) = 2 (p; p'}g(p) ap" , 
with 


E(p)= ¢Vc? M? + p?. 


Notice that the operator /(p, p') depends on the momenta not through 
p—p’ but through a more general dependence; accordingly its Fourier trans- 


o] 
Uci 


e. 
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form in co-ordinate space is of the form 7 (r'?—r™, r—r®”) and is there- 
fore a non-local operator. 

Strictly speaking a local potential cannot be defined, if it is to describe 
the scattering correctly. The scattering amplitude, as a function of energy is, 
for a local potential, analytic in the upper half-plane and on the real axis. 
If it is identical with the field theoretical amplitude, which is also an analytic 
function in the same region, on a segment of the real axis, it is also identical, 
on account of analiticity, on the whole real axis. 

But this cannot be true, since above the threshold energy, creation and 
annihilation processes do occur, and these, clearly, cannot be accounted for 
by the wave-mechanical amplitude. 

In simple terms we may say that, even under the energy threshold the 
scattering amplitude «feels » the presence of singularities due to creation 
and annihilation processes; at an energy under the threshold far enough from 
these singularities, however, it may be expected that their effect on the scat- 
tering amplitude is small and the scattering may be described with a good 
approximation by a local potential. That this statement is plausible can be 
seen from the treatment of spinless particles by CHARAP and FUBINI (14) who 
have shown that, for this case, the field-theoretical amplitude and the ampli- 
tude derived from a local potential are the same under the threshold when 
neglecting terms which are 0(p°/pù,) (*). 

In view of this fact the possibility of attaining the construction of a po- 
tential with a local character in quantum field theory, for separations 
larger than about one third of the pion Compton wavelength, appears somewhat 
impaired; it seems likely that it will be possible to do so only for larger sepa- 
rations. 

In general the potential V is defined as the limit of the operator Y° when 
p/Me <1 (this is called the « adiabatic limit »). 

For this discussion we have used concepts and quantities belonging to the 
S-matrix formalism; other methods for the construction of the potential, as 
it will be seen, may be used, still in the framework of field theory, which yield, 
as an equation equivalent to (1.1’), an equation of the form 


> € T L D If Je I ij 
(1.16) (2E(p) — W) ¢(p) = 55 | (Pp, p's W) o(p') dp’ , 
where the operator 7° explicitly depends on W, the total energy of the system. 


(“) J. M. Cnarap and S. P. FuBINI: Nuovo Cimento, 14, 540 (1959); 15, 73 (1960) 
see also M. Cini and S. FuBINI: Ann. Phys., 10, 352 (1960). 

(*) A similar proof for particles with spin, independent of perturbation theory, 
has not been possible up to now because the dispersion relations, on which the treatment 
of Charap and Fubini is based, have not yet been proved for this case. 


> 
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Several ways have been devised to eliminate this unwanted dependence, all 
of which, however, involve a large amount of ambiguity, 

All these methods, which fully employ quantum field theory, are called 
«non-adiabatic methods », in contrast to the old-fashioned methods, which 
were called « adiabatic methods ». 

These methods assumed the period of the nucleon motion to be large rela- 
tive to the time needed to exchange mesons; the two nucleons were then re- 
garded as approximately at rest during the exchange of mesons; the problem 
with fixed nucleons was then solved as a first step and the eigenvalue thus 
obtained, which depended parametrically on r, was adopted as a potential 
for the problem of the nucleon motion. In more precise terms the idea under- 
lying these methods was the limit w_/M — 0 (this is called the «static limit »). 
Now this limit is in general not justified, as CHARAP and FUBINI have shown and 
therefore the adiabatic methods (like any fixed source approach) can only 
have the meaning of a model approach, not of a limit of a complete field 
theory. 


1°4. — The basic attitude regarding the eq. (1.1') is well known: the existence 
of an analytic solution is assumed and this is built up by means of some series 
expansion. This is the typical attitude underlying the so-called « weak coupl- 
ing approach ». 

The potential shall therefore have the form of a series expansion. It is 
clear (as comparison with electrodynamics shows) that an expansion in powers 
of the mathematical coupling constant is meaningless since the terms of the 
expansion exhibit the well-known divergences. On the other hand an expan- 
sion in the number of exchanged mesons may be significant. It is hoped that 
such an expansion, at least in the adiabatic limit, is equivalent to an expansion 
in powers of the renormalized coupling constant and with all other quantities 
properly renormalized; then, apart from questions of convergence, one has to 
do with a «clean » series. It must be pointed out that in a series expansion 
of this type the problem of convergence appears, in the case of the nuclear 
potential, in less critical terms than in other applications of field theory. Wick’s 
classic argument (1°) in fact shows that the terms due to exchange of 3, 4, ... 
mesons will be important at separations between the two nucleons less than 
about #/3u,0, h/4u,¢,.., i.e. in a region which, as we have seen, must be 
excluded from our theoretical consideration. It is therefore reasonable to think 
that the first two terms in the expansion in the number of exchanged mesons 
give an adequate description of the interaction between the two nucleons in 
the region where the interaction may be described by a theoretical potential. 


(15) G. C. WicKk: Nature, 142, 993 (1938). 
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15. — In order to write the total field Hamiltonian we must take into ac- 
count on one hand the already stated general principles (see 1°1), except that 
the principle of invariance under Lorentz transformations must be substituted 
for the principle of Galilean invariance; on the other hand we must take into 
account the available phenomenological information. The range of nuclear 
forces and our qualitative knowledge on the interaction of mesons with nu- 
cleons compel us to exclude u mesons and to regard 7 mesons as those mainly 
responsible for nuclear forces. As for the heavy mesons, which we collectively 
denote as K mesons, it seems likely that they play an important role in 
nuclear forces, although their part is not yet quantitatively known. From the 
uncertaintly relation it is easy to see that if K mesons play a part in nuclear 
forces they do so to an appreciable amount within a radius of the order of 
r.—=thlu,e > luc. 

Meson theories have so far only taken into account x mesons; it is obvious 
that such theories can claim no validity for r< rx (*). 

Confining the consideration to 7 mesons (*), two other essential proper- 


- ties are given by experiment; spin (which is equal to 0) and intrinsic parity 


(which is negative). These bring us to regard the meson field as a pseudo- 
scalar field. Finally, if we add the assumption that the interaction is linear 
in the meson field variables and in their first derivatives, we can use two types 
of interaction Hamiltonian 


(1,17) Hi. = in | Pevanae, (ps-ps theory) , 
La TG Pol OD 

(1.18) ES © e fera Pa da, (ps-pv theory) (*), 
© 1% le 


where y, Y=y'y, are the nucleon field operators, D (a vector in charge space) 
is the meson field operator, 7 is the isotopic spin operator for the nucleon, 
z= pelh, g, and g, are the coupling constant for the pseudoscalar and the pseudo- 
vector coupling, respectively, y, =— 8x, (r=1, 2,3), va=f, Y5=V1YaV:Va (% 
and B are the usual Dirac matrices). 

The Hamiltonians (1.17) and (1.18) are in accordance with the so-called sym- 
metric theory of Kemmer (6), which ensures charge independence. Another 


(*) This is another reason for limiting the validity of the description of the inter- 
action between the two nucleons by a potential derived from pion theories to sepa- 
rations r>4 (flat). 

(**) Accordingly we shall hereafter write x without index for the pion mass. 

(**) A contact term may be omitted since it contributes a term of zero range to the 
two-nueleon interaction. 
('°*) N. KEMMER: Proc. Cambridge Phil. Soc., 34, 354 (1938). 
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way to ensure charge independence is to use the neutral pseudoscalar theory 
whose interaction Hamiltonians are 


(1.19) Hi, = in ] Brsp ax ; 


A 


ij CD 
(1.20) He a D =e 2 [y LP da, = dx , 


where ® is a scalar in charge space. 

Certains terms of the potential built up from (1.17), (1.18) contain the 
factor (7-7); if one substitutes this factor by unity one obtains the 
potential yielded by (1.19) and (1.20); we shall therefore always use (EG) 
and (1.18). 

We also quote the form taken by the Hamiltonian when assuming the 
source as fixed (« fixed source theory ») which has been used by many authors 
for the calculation of the potential (K= Me/h): 


; fi, ER j ; oe gi ira 
(1.21) Hy = Jp [ete VE Dax + ste [Nt Dda, 
(1.22) Hi, =" [ooo V)(T-)dx, 
x 
where 0(x) is the fixed source density, normalized in the usual way Jola x) dova: 


if, as frequently done, point sources are assumed, p(x) becomes 


It will be noticed that the forms (1.21) and (1.22) are not invariant under 
Galilean transformations. 
It is easy to see the physical meaning of the two terms of (1.21). The first 
term describes the interaction between a fixed nucleon and a p-wave meson. 
The second term has been the subject of debate (17%) a few years ago. 
It has its origin from processes of virtual nucleon pair creation and annihilation 


S. D. DRELL and E. M: HrnLEY: Phys. Rev., 88, 1053 (1952). 


7) 

ta G. WENTZEL: Phys. Rev., 86, 8C2 (1953). 

(19) K. A. BRUECKER, M. GELL-MANN and M. L. GoLDBERGER: Phys. Rev., 90, 
476 (1953). 

(20) N. M. Krozz and M. A. RUDERMAN: Phys. Rev., 93, 233 (1954). 

(21) S. DESER, W. E. THIRRING and M. L. GOLDBERGER: Phys. Rev., 94, 711 (1954). 

(22) A. KLEIN: Phys. Rev., 95, 1061 (1954). 

(23) S. F. Epwarps and P. T. MatruEws: Phil. Mag., 2, 176, 467 (1957). 
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and would seem to dominate the interaction. Arguments have been brought 
forward, however, which support the belief that the effect of these processes 
is largely overvalued by the second term of (1.21) as compared with their treat- 
ment by means of the complete Hamiltonian, at least at not too high energies. 
In other words these processes must be strongly damped (perhaps by a factor 10), 
if the description of the interaction by means of the Hamiltonian (1.21) is to 
be retained. The arguments in favour of a strong damping are of a qualitative 
rather than a quantitave character and, although supporting the belief that 
these precesses are not of decisive effect in the mechanism of nuclear forces, do 
not supply a definite prescription for the evaluation of their role. The Hamil- 
tonian (1.21) is therefore modified by the inclusion in the second term of a 
factor A, presumably very small: 


(1.23) H,. = L Jette: Vig dae + À ei Jeter dax. 


The contributions to the potential given by the second term will be easily 
identified by the factor A. 

In the case A=0, #.e. in case the second term of (1.24) is negligible, the 
Hamiltonian (1.21) becomes the Hamiltonian (1.22), provided g,/2K = gs/x. 

This means that up to a certain order of approximation there are processes 
which are described in an equivalent way both by the ps coupling and by the 
pv coupling; this property was known as « Dyson’s equivalence theorem » (24). 

Of these «fixed source theories » we must however remind what we said 
when discussing the static limit (u/M — 0). 


2. — Methods for the construction of the potential. 


In the past the methods more frequently employed have been the adia- 
batic methods (25). According to these methods the potential was defined as 
the lowest eigenvalue of the Hamiltonian 


Jah di Bagle He ’ 


where H,, is the free meson field Hamiltonian and H/, is the fixed source 


(24) F. J. Dyson: Phys. Rev. 73, 929 (1948). 

(25) N. KEMMER: Proc. Roy Soc. of London, A 166, 127 (1938); H. FRÔHLICH, 
W. HerrLER and N. KEMMER: Proc. Roy. Soc. (London), A 166, 154 (1938); C. MOLLER 
and L. ROSENFELD: Kgl. Dan. Vid. Selskab., 17, 8 (1940); G. WENTZEL: Quantum theory 
of fields (New York, 1949), p. 63; J. IWADARE: Progr. Theor. Phys., 18, 189°(1955); 
14, 16 (1955); C. MoLLER: Pseudoscalar meson theory; CERN, 55-19 (1955). 
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interaction Hamiltonian (given by (1.21) or (1.22)). The calculation of this 
eigenvalue was attained either by means of ordinary perturbation theory or 
by means of canonical transformations. Since the static limit was inherent 
to these methods they are now felt to be inadequate and obsolete. 

There are various non-adiabatie methods, which may be classified accord- 
ing to four approaches: a) Tamm-Dancoff method; b) S-matrix methods; 
c) method of Brueckner and Watson: d) method of the Bethe-Salpeter equation. 


21. Tamm-Dancoff method (26:27), — In the Tamm-Dancoff method the 
equation (1.1') is transformed into an infinite set of simultaneous integral 
equations for the amplitudes a” (A) of the free Hamiltonian eigenstates (where 
m is the meson number, # the number of nucleon-antinucleon pairs, À refers 
to the momenta, spins and isotopic spins). 

The eigenvectors |m,n, A> of the free Hamiltonian are defined by the 
equation 


(2.1) H,|m, n, A> = H”(A)|m, n, 2) 


and form a complete orthonormal set. The state functional |D) may then 
be expanded into a series of the Im, n; A>: 


(2.2) |B> = dar 


Mm 


(A), 2, A>. 


mn; À 


Substituting this expansions into the equation (1.1’) one obtains 


(2.3) > a®(A) En (A)|m,n, D + Y a" (A) H'|m, n, A’ = W > a2(a)|m, n, a> 


m,n À m,n, À mn À 


and from this the following set of integral equations : 


(2.4) [W—2,(0)]a,() = ZX <m, n, AH" |p, q, pat) . 


Pg: 4 


If the interaction H' is assumed to be the ps-interaction (*) 


H' = in| Pryay@ax, 


(25) I. Tamm: Zurn. Éksp. Teor. Fiz., 9, 449 (1945). 

(2°) S. M. Dancorr: Phys. Rev., 78, 382 (1950). 

(*) The ps coupling has the advantage that it allows the renormalization to be 
carried out with Dyson’s techniques, what cannot be done for the py coupling theory. 


a 
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this set takes the following form: 


n+1 


"(4) = > D D <m, n, À| H’ |p, q, ur au) : 


m 


(2.5) W— E"(2)|a 
p=mEl q=n—1 4 


The set (2.5) is a mathematically exact formulation of the problem, but it 
involves an infinite number of equations. An approximate way to solve it 
is that introduced by TAMM and DANcorr. The approximation introduced lies 
in putting equal to zero the probability amplitudes a” for n> N and m> M 
(with M, N fixed numbers). A set of a finite number of equations is thus 
obtained which may, in principle, be solved in an exact way. 

For the consistency of the method it is to be hoped that the solutions con- 
verge for N, M — co; in other words that the results of the calculation are 
insensitive to the number of particles N, M when N, M are large enough and 
tend to a finite limit when N, M — co. Unfortunately this is only a hope 
because nobody has ever solved the set of equations but for the case of very 
few components and the assumption therefore lacks any support. 

A second method to solve the set (2.5) is that used by Livy (28). This 
consists in the successive algebraic elimination by induction of the amplitudes 
of the states containing mesons and nucleon-antinucleon pairs. What remains 
after this elimination is an integral equation for the amplitude of the zero-meson 
zero-pair state, in which the kernel has the form of a power series expansion 
in the coupling constant. In the actual solution the kernel is then cut at a 
certain order of the coupling constant; this implies the assumption of the 
convergence of this expansion. The equation thus obtained has the following 
form: 


(2.6) [W — BU) () = D KA, 4; W)ad(u) . 


The kernel K as the form 


(2.7) K= > K, (A, uw; W), 
with 
(2.8) EAS APCE 


Eq. (6) appears like a Schrédinger equation; there are reasons however which 
prevent it to be interpreted as such. In fact the a? appearing in the equation 
cannot satisfy the normalization condition, since this for the system in ques- 


(5) M. M. Lévy: Phys. Rev., 88, 72, 725 (1952). 
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tion is 
n E 
(2.9) MILA a 
min, À 


and the condition 
(2.10) D|ag(4) ress 


cannot in general hold. 

Furthermore the operator in the integrand depends on the total energy W. 
The procedure of Lévy and Klein (?°) to eliminate the energy dependence brings 
about an interaction which is in general non-Hermitian. This can be seen for 
instance in the lowest approximation; in this case the integral equation for 
the aj(A) is 
pene DV 2 AYA) = > <0, 0; Al A’ 1, 0, 4 1; 0, | 2°10, 0; w au). 

W— Et(u) 


7 
pupi 


The energy denominator is expanded in the following way: 


; 1 eo Eu) 
(2.12). [WET = per EIA a Ca a È | 
Eu) — Eu) no \ Eau) — Hole ) 

Substituting this expansion into (2.11) we can use à method of successive 
i racti thie à iste | a ti 0 7 0 0 : : RA 
interaction which consists in substituting (W — Bu'))at(2') in the second 
member of the equation with its value given by the equation itself in the preced- 
ing iteration. One obtains then an integral equation for the a} with a kernel 
no longer depending on the energy and having the form of a power series ex- 
pansion in the coupling constant. The second and fourth order interactions 
are, written in a symbolic form, 


== Hd, ’ 
(2.13) ; 
Va = HAE He He) oy 


where H,,, Ho, couple the zero-meson and one-meson subspaces by creating 
and annihilating one meson respectively; clearly it is 


(2.14) E 

H,, and H,, are defined by 

(2.15) | <a Ha lu) = <p! [Hm | > (BCH!) — E), 
2.15 

| cu x |) = <a | Holu (Pol) — E), 


(29) A. KLEIN: Phys. Rev., 89, 1158 (1953). 
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(O*, with O any operator, has an exactly similar meaning). V, and V of (13) 
are not Hermitian, since 


(2.16) (H*) = — H} 


10? 


Hi HE 

They turn out to be Hermitian when retardation effects are neglected, at 
least in cases of particular interest. Nevertheless the objection of their non- 
hermiticity in principle remains. 

FELDMAN (3°), making use of the formalism of canonical transformations, 
showed that it is possible to decouple the amplitude a} from all other ampli- 
tudes up to any desired order of the coupling constant by means of a sucées- 
sion of canonical transformations. This procedure will have physical meaning 
when the conditions, stated in the introduction, which make valid the de- 
scription of the interaction by an ordinary potential, are verified. One obtains 
then an equation for a} with a Hermitian kernel which may be interpreted as 
a Schrôdinger equation for the two nucleons. 

The results of Feldman’s method with the ps coupling Hamiltonian in the 
fixed source model are, with terms up to the order (u/2M)?, the following: 


REI D 7 ny 
(2.17) ne 2 (7) (c+ 7) (6+ V)(c®-V) exp [— ar] La 


O Polo e one 
Pale panera (tro) Cer Oe EEE 


(39) D. FELDMANN: Phys. Rev., 98, 1456 (1953). 
(*) In addition to the terms listed there is a term 


gi fu \ ee 
ue —— | —— | (7-72) (g@™ +g) —— br 3 
{ el PR rotte 


obviously it will be omitted since the validity of a theoretical potential is limited to the 
region x > Xin» 

From an historical point of view it must be pointed out that the singularity 1/7? 
of the potential V, was viewed as a serious difficulty and various devices were con- 
trived in order to overcome it. Among these we quote the mixed theory of Rosenfeld 
and Moller (31), Bethe’s cut-off (32) and the non-linear theories (33). In the light of the 
arguments stated at the beginning it appears clear that this difficulty fades when one 
bears in mind that a theoretical potential cannot hold for nucleon separations x < pe 

(*) €. MeLLER and L. ROSENFELD : Kgl. Dan. Vid. Selskab., 17, 8 (1940); C. MoLLER: 
Kgl. Dan. Vid. Selskab, 18, 6 (1941). See also J. M. JAUCH and N. Hu: Phys. Rev., 65, 
289 (1944). 

(2) H. A. BrrHz: Phys. Rev., 55, 1261 (1939), 

(8) L. I. Scurrr: Phys. Rev., 84, 1, 10 (1951); 86, 856 (1952). 
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Or 


(2. 1 8) Vi = Ve sù Va , 
Sata 2 oo eS ten) 
2.184 ie Soap es | La 
| a) - 5 (£ A a È Ge ie 

2 r Lai È 1\? exp [— 2a] 
2.185 Vis = Gye? 2 ( Le dela dl 
| i 4 i) 2M a x? È 


where 2 = xr; K,(x) is a modified Hankel function of the first kind: K,(x)= 

= (x/2)i" 1 H(ix). The potential V, and V, of (17) and (18) are in agreement 
with those previously obtained by LÈévy and KLEIN. 

The term V,, has its origin in the two-pair processes; the term V, in 

the one-pair processes; the no-pair processes bring about a term of the order 

(u/2M)* which is therefore here neglected. 

The potential V=V,+ V, given by (17) and (18) is wholly inadequate, on 
account of the near-cancellation between the attractive term V,, and the 
repulsive term Vy. 

The pv coupling Hamiltonian yields in the fixed source model the T.M.O. 
potential (see (2.43), (2.44)). 

An equivalent method to overcome the difficulties above has been used 
by FUKUDA, SAWADA and TAKETANI (**) (*). 

Denoting by P the projection operator on the two-nucleon subspace, we 
may write 


(2.19) |®> = P|®> + 1— P)|®>. 


Projecting the eq. (1.1’) on the two subspaces one obtains the following 
two equations: 


(2.20) (W— H,)P|®> = PH'P|®) + PH’(1— P)|®), 
(2.21) (W—H,)(1—P)|®) = (1—P)H'(1— P)|®> + (1— P)H'P|®). 


A formal solution of (1.1') may be written 


(2.22) [Dyer WP D, 
with 

1 
(2.23) J(W) = 


LUNA) PH" 


(34) N. Fuxupa, K. SAWADA and M. TAkeTANT: Progr. Theor. Phys., 12, 156 (1954). 
(*) A very similar procedure was devised by OxusBo (°°). 
(35) 8. OkuBo: Progr. Theor. Phys., 12, 603 (1954). 
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The energy dependence of J(W) may be eliminated making use of the equa- 
tions of motion. 
It is evident that 


(2.24) <OPJ'|JI PD) = (O| Oy =1. 

This suggests a linear transformation 

(2.25) Lo = JJP|E) 

in the two-nucleon subspace, assumed to be complete; the |y) are correctly 


normalized 


<tal 45 == Ô ad 
and satisfy the equation 
(2.26) (W—-H,)|x>=[P(J*7)*PJI*H,+H')JP(J'J)*P—PHP]}}|x}; 


which may be interpreted as a correct Schrôdinger equation for the two phys- 
ical nucleons. At the static limit (H,=0) the interaction kernel yields the 
static potential 


1 


(2.27) V= Pag 


PPIH E HWP s 


The factor P(1/J'J)P is called the «probability operator» inasmuch as its 
expectation value represents the probability for two «physical» nucleons to 
be bare 


bo 
(o ©) 


1 i o 
(2.28) GP REATI) ID RE DIE 


It is subject to the physical condition 


1 
2.29 0< <y|=— |?) < 1. 
The factor 
PJ'(H,+H')JP 
is called the «normal part ». 


In order to obtain the potential up to the fourth order in the coupling 
constant one must expand (2.27) into a power series in the coupling constant. 


1 
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Different forms of the potential obtain depending on the way of treating 
the two factors. Putting 1/(J*J)=1 and expanding the numerator up to the 
fourth order one obtains the Brueckner and Watson potential (see (2.74), (2.7 5)) 

Expanding the factor 1/J*J up to the second order and the normal part 
up to the fourth order (in the expansion of 1/J'J the first term is 1, the second 
term is of second order in the coupling constant) one obtains the T.M.O 
potential (see (2.43), (2.44)). In this case however the condition (2.29) be- 
comes 
1— P 


(— H,)? HEP | <a il . 


(2.30) EER H' 

Now while (2.29) is automatically satisfied, (2.30) like any condition where 

(J*J)-1 is approximated by a series expansion cut at a certain order, must be 

explicitly imposed as a supplementary condition and can thus be a useful 

guide in the discussion of the region of validity of the potential obtained. 
F.S.T. propose as static nuclear potential 


ghee la e er rg Lesbo 
(2.31) V= PET ates Hoey et hs mire 

7 re 1-P H' \ 

{ 

\ (Ho)? ff 
(<...> means meson vacuum expectation value), 4.6. with the denominator 
expanded up to second order and the normal part up to fourth order. The 
argument for this prescription is the following; for the fixed source model of 
the neutral scalar theory, the exact solution is known (**); this solution is 
reproduced by the method above only if the prescription is adopted. Further- 
more the presence of J'J in the denominator reduces the singularity of the 
normal part, as can be intuitively seen considering that a strong singularity 
arises when the meson cloud around the nucleon is intense and, accordingly, 
the probability for the nucleons to be bare becomes much smaller. 

It must be pointed out that the explicit evaluation of the potential taking 
into account the probability factor, is very sensitive to the renormalization 
procedure (cut-off or other); the validity of condition (2.29) (or an analogous 
one in which an approximate form of (J*J)-1 is used) requires the inclusion 
of the self-energy terms which cannot be simply ignored. 

The detailed calculation along the lines set up by F.S.T. has been carried 
out by INOUE, MACHIDA, TAKETANI and Toyopa (37), starting from the pv 


(35) See e.g. WENTZEL: Quantum theory of fields (New York, 1949), p. 47. 
(37) K. Inoun, S. Macuripa, M. Tagerani and S. Toyopa: Prog. Theor. Phys.. 
15, 122 (1956). 
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Hamiltonian (1.22) and choosing 


CR) = exp | ca | 


32 
Menace 


for the Fourier transform of the source function describing the nucleon. 
The resulting potential is 


1g? 12 F | k 
(2.32) V=3 da — (TMT) | 6? (k,)[(0 6) (2A, + B,)—S32(A,— Bi) ] — dk, 
Sea aE Wi 
0 
ga\? de f2\2 ar 2 Ae 
_ ( 2 = 2) SPO [2(T®-7®)(B,B, + 2A,A,) + 
47) hier \n DIRLO AA 
0 0 
+ 2(6-o)(A,B, + BA. + 4,42) — Sig(4:B,4+ B,A,— 24:4:)]° 
02(k,)c' (ky) ky ks [ SE SRI ; 
—, he (RE 3A : (1). 7) 
WOW + ali eo | | VBA MORE ap LE 
à a l ki | 
| [e {(o-0)24, zie B,) = Sp(Ai A B,)| gi ’ 
; 1 
where 


1 : 
A = x (k,r cos kyr — sin k,r) , 


B, = n (2 sin kr — 2k,r cos k,r — kîr? sin kr) , (Ciel Zee 


Oe yey + ke : n 


Viù a . . 
and A, is a renormalization constant (*). 


(*) Following the method of Chew (‘8) I.M.T.T. carry out the renormalization by 
the formula 


7-1 
J2= G27122 , 


where G, is the renormalized coupling constant and Z,, Z, are renormalization constants 


for the vertex and the nucleon propagator, respectively; up to the second order in the 
coupling constant one obtains 


x 


Lan 3 GE I ° a2( Io) 2 
7 il SU: 7 plea le “itl pee Ta c?( Li 
i GS 3 AA A, hol Qa Ro: dk. 


(5) G. F. CHEW: Phys. Rev., 95, 1669 (1954). 
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The most serious objection which can be moved to the T.D. method is that 
within its framework the renormalization program cannot be carried out in 
an unambiguous way. Even Cini’s (°°) covariant renormalization method is 
not free of ambiguities (*°); nor completely renormalizable is the «new Tamm- 
Dancoff method » proposed by Dyson (4!). The essence of this method is the 
re-establishment of symmetry between emission and absorption (or between 
particles and antiparticles), which is lacking in the old T.D. method, through 
the definition of a new probability amplitude 


ao(N EN) == <0, CIN) AN) ®O>, 


where |P) is the actual state of the system, |.) the true vacuum, A(N) is 
the product of free particle annihilation operators which annihilates the par- 
ticles specified by N, and C(N') is the product of the corresponding creation 
operators. 

The new T.D. method has been applied to the nuclear force problem by 
TAYLOR (22) and KURSUNOGLU (4%) and at the static limit yields the same 
results of the original T.D. method. 


22. S-matrix methods. — The S-matrix formalism is one of the most widely 
used for the construction of the potential. In the past this formalism was 
favoured because it supplied the tool for the treatment of renormalization 
problems and connected subjects and more recently because it has opened 
the gate to the introduction of dispersion techniques in this field. 

This shift of emphasis is reflected in the use of the formalism by the various 
authors in the course of time, as it will be seen from the following account, 
in which a historic line will be followed. 

We point out from the outset that the potential yielded by these methods 
is a potential for scattering states. 

The common feature of the definition of the potential in the framework 
of the S-matrix is the identification of certain quantities describing the scat- 
tering of two particles in particle wave-mechanics, such as the s-matrix ele- 
ments, the f-matrix elements, the k-matrix elements, with the analogous 
quantities such as the S-matrix elements, the T-matrix elements, the A-matrix 
elements, describing the same situation in quantum field theory. 

The actual construction of the potential is carried out in different ways 


) M. Cini: Nuovo Cimento, 10, 526, 614 (1953). 

) F. J. Dyson: Phys. Rev., 91, 421 (1953). 

(41) PF. J. Dyson: Phys. Rev., 90, 994 (1953); 91, 1543 (1953). 
) J. C. TAYLOR: Phys. Rev., 95, 1313 (1954); 96, 438 (1954). 
) B. KursuNoëLU: Phys. Rev., 96, 1690 (1954). 
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by the various authors, depending on which quantities aré chosen for the 
identification and how the latter is performed. 

We briefly list the quantities of interest and the relations establishing a 
link between them, which are written in the center-of-mass system. 

From an S-matrix element between nucleon states |a) and |b> a T-matrix 


element is defined by: 
(2.33) <a|S|b> = dan — 27 6(H, — E,)<a|T|b> . 
In virtue of the unitarity of the S-matrix, 7 satisfies the following equation: 


Pr 


Na 
CC 


(2.34) 27 > ÔE,— E,)<a|T* 


Scattering processes may also be described by the hermitean operator .X which 
is linked to the S-matrix by the equation 


TEA 
(2.35) eee ee 
1+3iX 


Besides the #-matrix, the reaction matrix K is also used: 
(2.36) {a|K|b) = 2x Ô(E,— E)Ka|X]|b). 


Completely analogous quantities appear in particle wave mechanics and 
will be indicated by small letters. 

We finally note that ¢ satisfies the Lippmann-Schwinger equation (44): 
(2.37) TESTI ke AMATO 
2% | € E, = E, + 17] 

The old-fashioned construction of the potential was carried out according 
to the following procedure; the Lippmann-Schwinger equation was solved by 
a perturbative method: 


(2.38) t=t,+%+..., 
where 
(2.39) <a|t,|d) = al)", 


4 | (as \ ) / yee | 
(2.40) <a|t,|b> = <a|V,|b> + > else SEE » 
e HE = E, + V7] 
(85) See e.g. L. HULTHÉN and M. SuGawara: Handbuch der Physik, vol. xxxIx 
Berlin, 1957), p. 135. For the formal theory of scattering see also: A. S. DAvYDov: 
Teoriya atomnogo yadra (Moscow, 1958), Chapt. 1x. 
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ete.: each of these terms was then equated to the static limit (u/M — 0) ot 
the field theoretical terms 7, T,,..., arising from the one-pion, two-pion, ete., 
exchange graphs: 


(2.41) <a | ty | Da IE | D> stat ; 
(2.42) <a |t, | b> = <a Ta | D> stat 9 
ete. . 


A few remarks are in order at this point. We note first that the T-matrix 
is defined only on the energy shell. While the matrix elements on the energy 
shell are sufficient for the definition of the one-pion exchange potential, matrix 
elements off the energy shell are necessary for more than one-pion exchange, 
as can be seen for the two-pion exchange from eq. (2.40). The method can 
therefore yield a whole class of potentials arising from the same S-matrix. 

Moreover the crude identification assumed by the method ignored certain 
difficulties: in particular the static limit was carried out on the T-matrix ele- 
ments and these elements were then equated to the t-matrix elements without 
examining if the necessary conditions were satisfied. Now, as we have noted 
in the introduction, the static limit is not justified; moreover the t-matrix ele- 
ments satisfy certain conditions, which are not satisfied in general by the 
T-matrix elements. 

In spite of its shortcomings the method has been widely used in the past 
by various authors, some of which have also undertaken a calculation of ra- 
diative corrections of various orders and an assessment of their effect on the 
one-pion and two-pion exchange potentials. 

As examples of these researches we quote the contribution of WATSON 
and LEPORE (4), who by means of an evaluation of the lowest order radiative 
corrections to the one-pion potential, showed that these reduce to coupling 
constant and mass renormalization and that of HirpA, IwADARE and MA- 
CHIDA (4), who reached the same conclusion but taking into account radia- 
tive corrections of all orders; the calculation of the potential by I. SATO (4) 
for both the ps-ps coupling (giving the Lévy-Klein potential) and the ps-pv 
coupling (giving the T.M.0. potential) and the compact formulation of the 
problem by IMAMURA (48). The best known result of this group of research 
was the T.M.O. potential (*) which enjoyed a certain fame for several years 
and which results from the explicit calculation of the expression of NAMBU (°°), 


(45) K. M. Warson and J. V. Lepore: Phys. Rev., 76, 1157 (1949). 

(46) K. Hirpa, J. IWADARE and S. MacHipa: Progr. Theor. Phys., 15, 189 (1956). 
(47) I. SaTO: Prog. Theor. Phys., 10, 323 (1953). 

(48) T. IMAMURA: Progr. Theor. Phys., 18, 183 (1955). 

(49) M. TAKETANI, 5. Macuipa and 8. ONUMA: Progr. Theor. Phys., 7, 45 (1952). 
(50) Y. NamBu: Progr. Theor. Phys., 5, 614 (1950). 
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using the ps-pv coupling: 


Go 
9 4° = ee (TETE) 
(2.43) Vee ei 
= 1 3 BO\ (oe |= a 
(REA (1 : aa = | 
(2.44) a = eH [TP 1) U(æ) + (6P-0@) U(2) + Sr Ur(t)], 


(2.44a) PACA =\(; + =] K,(2x) + (È at =) Ke. 


1 4 G2 4x 
pnt ie one 2 3 La) 
(2.44b) U,(0) = x | K(2@) + È i SI 
813 15 
(2.44c) U(x) = — aE Kn) (= 23 ie K, (£ n) 


A more elaborate version was given by HENLEY and RUDERMAN (51), with an 
attempt of a calculation of multiple scattering corrections; they showed that 
it is possible to write the S-matrix element for the two-nucleon scattering in 
terms of matrix elements for meson-nucleon scattering. 

Using the same idea, MIYAZAWA (5°), treating nucleons as fixed sources 
coupled to the meson field by pv coupling, was able to express the potential 
in terms of the experimental quantities for meson-nucleon scattering, making 
use of dispersion relations to define the meson-nucleon scattering amplitude 
for all values of the energy from minus infinity to plus infinity (real scattering 
obviously occurs only for energies larger than the pion rest mass). 

The potential obtained by Miyazawa, which was explicitly calculated by 
Konuma, MIvAZAWA and OTSUKI (5), is remarkably different from those pre- 
viously obtained; we do not quote it here on account of its length. 

More recently the problem has been formulated on a more critical basis. 
One point which has been definitely clarified is that the static limit is comple- 
tely unjustified. In spite of this, it has been shown that the construction of 
a potential up to two-pion exchange is nevertheless still possible. The dis- 
cussion of the possibility of the definition of the potential and of its limits and 
the actual construction of this can be done with various techniques; in the 


(5) E. M. Henry and M. A. RUDERMAN: Phys. Rev., 92, 1036 (1953). 
(°°) H. Miyazawa: Phys. Rev., 104, 1741 (1956). 
(93) M. Konuma, H. Miyazawa and 8. Orsuxt: Progr. Theor. Phys., 19, 17 (1958). 
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introduction we have given an account of one of these more modern approaches 
due to BRAUN (#). The construction of the potential is based, following 
the suggestion of NovozHILov, instead of on the perturbative procedure on 
the Lippmann-Schwinger equation, on the use of a non-linear equation which 
results from a formal solution of that equation, which gives the potential ex- 
plicitly in terms of the #-matrix. As we have seen, it is actually possible, 
although not in a unique way, to obtain from quantum field theory certain 
T-matrices which are identical, under the creation threshold and on the energy 
shell, with the T-matrix and which possess the necessary requirements in 
order to be inserted in place of the t-matrix in the non-linear equation, thus 
yielding an expression for the potential. We record here the formal expres- 
sion of the «potential » thereby obtained: 


~ 


(2.45) Vas = {1 + NV (By, Eo)1 + NT the, 
= E,+ E 
(2.45a) Ve he, By) OMR 


Ls rh, (EL <-> 3 (By + E;)) 


eine, 
We recall that this is a non-local potential; we refer to the introduction for 
a discussion of this feature. 

An alternative formulation of the problem is due to Gupra (54), who used 
the #° operator linked to the S-matrix by eq. (2.35). If X is expanded in the 
following way 


(2.46) KH =X, 


and this expansion is introduced step by step in (2.35), one obtains a « uni- 
tary expansion » of the S-matrix, since it is unitary at any stage, unlike the 


usual Dyson expansion 


(2.47) ISIS 


which is unitary only as a whole. 
The #, can be easily expressed in terms of the S,; we record here only 
the expression for #, and #, tor the two-nucleon scattering: 


HEI a 


n . vt 
A,= 48, 8). 


(54) S. N. Gupra: Phys. Rev., 117, 1146 (1960); Nuovo Cimento, 18, 823 (1960). 
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Taken the unitary expansion of S in first approximation an effective inter- 
action Hamiltonian H,, can be defined such that the contribution of H., to 
the unitary expansion in first approximation is equal to the total scattering 
operator 


1+4(—iff)[Hudxdt 1-210 
1—4$(—i/h)(Hudedt 1+ 41%” 


which immediately yields 


(2.48) the (1/0) [Hada 


The «potential » may be defined by expressing the effective interaction 
energy of the two nucleons, in the center-of-mass system, in the following way: 


(2.49) <a | Hace dt|b> = 27h0(E,— E,)<a|V |b> . 


It is easily seen that its definition in terms of the “-matrix is 
(2.50) {a|X |b) = 2x Ô(E, — #,)<a|V |d>. 


The one-pion and two-pion exchange potentials can be easily obtained from this 
taking the adiabatic limit of #2, #,. 

The result of the calculation with ps-ps coupling for the one-pion exchange 
is in agreement with the usual expression; the two-pion exchange term is 
given by 


(2.51) vr=— (i uo[P(a)— (12-12)Qa)], 
where 
3M 1 u 
9.51 ; — SM? Oe at y 
Re) P(æ) = sa [av far (u — v)2(1 — uw)? FP [S72], 
ie ai ; v : 2 
(2.51b) ala 3 fo far PR en exp [— fa] 3 P(e) ; 
with 
(2.51¢) j= (ert CRS 
(uv) 4%) 
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One might think of following a procedure analogous to that described above 
for the T-matrix consisting in identifying step by step the terms of the k-matrix 
of particle wave mechanics with the terms of the Æ-matrix of field theory. 
This would obviously result in a potential identical with that defined by means 
of the T-matrix but different, in general, from that obtained from Gupta’s 
procedure. While the one-pion exchange term of Gupta is still identical with 
that of the previous method, the two-pion exchange term differs from that 
of the previous method by a principal part. It is easy to see that the two-pion 
exchange term in the potential defined by identifying K, with k, is 


<a|V4|b> = <a|K,|b> — > <a K, ©? ( ) cel |0 E 


if 
E ra E, 
while Gupta’s two-pion exchange term is 


<a 


Vilb) = <a| Ky] b>. 


On account of this fact Gupta’s two-pion exchange potential contains, 
unlike the potential defined in the preceding way, only 7-matrix elements 
on the energy shell, as can be seen by comparing the expression of the two 
potentials in terms of 7-matrix elements: 


{a|T,|c) e] T,|b) 


ANA ATO) ee rr 


e 


<a|T,|b> + ix Y <a|T,|e> (E, — E,)<o|T, |b» . 


I 


CAD 


In the foregoing expressions the adiabatic limit is to be intended whenever 
T- or K-matrix elements appear. In the old-fashioned methods use was made 
as a rule of the static limit; it can be shown that this is needless because the 
potential can be constructed up to two-pion exchange making use of the adia- 
batic limit only. That this is unjustified has been repeatedly emphasized and 
is shown in detail in Gupta’s paper. 

The recent interest for dispersion relations and techniques in theoretical 
physics has suggested a promising and powerful approach to the problem of 
the nuclear potential. 

Although persisting difficulties still prevent the integral extension of the 
method to the case of nucleons with spin interacting through pseudoscalar 
symmetric meson field, the cases of the scalar nucleons interacting through 
neutral scalar meson field or charged scalar meson field can be exhaustively 
discussed (11). 

We give here a brief account for the case of the neutral scalar field. Let 


a 13 - Supplemento al Nuovo Cimento. 
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us assume that a dispersion relation at fixed momentum transfer 


dk’, 


< : 1 f Imt(E', t) 
(2.52) t(E, T) = } (tT) nr | E'— = in 


0 


(where 7 is the momentum transfer) is available. From this the following re- 
lations are immediately obtained : 


foe} 


a E, t) = t(0 ss nes 
(2.53) t(H, t) = (0,7) + af EEE =ir 


dE’ , 


0 


dk’. 


- fl (E t(E',t) 
9 À Vi = n 
(2.54) V(t) (Oy we) a El 


0 


Beside these we consider the condition of unitarity of the s-matrix: 


(2.55) icp'le- tp) + ake tp") ò (A <p"|t|p> dp = 0. 

ì ; hs Mat 
On the other hand for the field transition amplitude G, which differs from the 
previously defined 7 by a factor M?c‘/E?, which is about 1 for p?< M?e?, 
a dispersion relation holds, and a unitarity relation, which for energies lower 
enough than the threshold, may be written in the following way 


: : L CH Im “G(E, t) E 
2.5 ‘E, t) = | di 
(2.56) CE, t) G(0, t) 4 2/7 L'-E—-in Le Fee 


0 


12 


a . 1er "I 27 pena " È I 
(2.57) | i<p'iG’— CIp\at = |*|p (E) <pri@|p>ap’=0, 


where Im “G is the contribution to the imaginary part of G coming from 
elastic processes. 

If one identifies ¢(0, 7) with G(0, 7), it is seen that ¢(#, 7) on the one hand 
and G(E, 7) on the other satisfy (when neglecting 0(2/2,,.a) the same dis- 
persion and unitarity relations. This identification provides therefore a rea- 
sonable basis for the definition of a potential from quantum field theory. Sub- 
stituting G(0, 7) for #(0, 7) and solving the non-linear set of eq. (2.53) and (2.55), 
eq. (2.54) provides a definition for the potential. 

This program may be carried out for scalar nucleons interacting through 
a neutral scalar meson field. In fact for this case a dispersion relation of the 
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type (2.52), known as Khuri’s relation, is available and a method for solving 
explicitly the set of eq. (2.53) and (2.55) may be devised; the result is a po- 
tential which is a superposition of Yukawa potentials with a spectral function 
giving the mass spectrum, when neglecting the so-called rescattering correc- 
tions (*). The one-pion exchange potential is the well-known Yukawa po- 
tential; the two-pion exchange potential has an expression such that it would 
vanish in the static limit, if this were justified. 

The program may also be carried out when the scalar spinless nucleons 
interact through a charged scalar meson field. 

In fact t(#, t) may then be split into two parts, each of which satisfies a 
Khuri relation; an analogous splitting for the field theoretical G(Æ, 7) enables 
an identification to be performed in a way similar to that described above. 

When the spinor character of nucleons is taken into account, the program 
outlined cannot be carried out because in this case the necessary dispersion 
relations are not yet available. One might be confident in the eventual avail- 
ability of these relations and therefore one might postulate them; however the 
transition amplitude is made up of five terms corresponding to five spinor in- 
variants, the choice of which is however not unique; while one might expect 
the dispersion relations to hold for a particular choice of these only and not 
for an arbitrary one. Under these circumstances CHARAP and TAUSNER (56), 
returning to the old definition based on the step identification of the t-matrix 
elements and the 7-matrix elements, have shown that the condition that the 
two-pion exchange potential be independent of the energy in the adiabatic 
limit leads to restrictions on the choice of the spinor invariants. This con- 
dition may therefore be deemed a useful guide for the choice of these. 


23. Method of Brueckner and Watson (57). — The method employs the Lipp- 
man and Schwinger integral equation and the technique developed by CHEW, 
GOLDBERGER and WATSON to treat this equation (58). 

Let |y> be an eigenstate of the free Hamiltonian. A scattering solution 
of eq. (1.1’) may be written 


(2.58) ID) = 21»), 


with © the Moller wave matrix, satisfying the Lippamnn and Schwinger 


(*) For recent progress in the treatment of these corrections see CINI and FUBINI (9°). 
(55) M. Cini and S. Fusini: Ann. Phys, 10, 352 (1960). 

(56) J. M. Cuarap and M. J. TAUSNER: Nuovo Cimento, 18, 316 (1960)). 

( 


58) G. F. CHEW and M. L. GOLDBERGER: Phys. Rev., 87, 778 (1952); K. M. WATSON: 


Phys. Rev., 89, 575 (1953). 
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equation 

> 59 O fe I ¥ Ù 
2.09 (== Sem 
Wop oy 4 


whose formal solution is 


(2.60) Q=14 Wen a 

In the case of nuclear forces we are interested in the diagonal part of ©, 
which refers to the scattering of two-nucleons, without creation or annihilation 
of real mesons, and we need to express this part in terms of operators which 
are diagonal in occupation number. 

This expression for 2, which we denote by Q,, is reached by an inductive 
procedure which we relate for the case of an interaction Hamiltonian assumed 
to be diagonal in nucleon occupation number and non-linear in meson field 
variables ®. This Hamiltonian can be written 


(2.61) HE te ap 


where H, is odd in ® and H, is even in D. We then obtain for © 


2.62 Q=1 Sb ee apt 
(2.62) gh 
sul : (H, + Do) 4 MEET AT de; ; 
NE Te ae Dl gia Fag hee ss a sù) 2| (*) 
where 
. Il 
(2.63) ID = Eb tale > A = W = He (+) : 


DETTE 


The last term of (2.62), odd in H, (and thus in ®), is non-diagonal in meson 


(*) From the operator relation 


1 l 
RS ay LE TEA io + Ho) - A HI, 


1 
L= (H + H;) A (H, Hi Hy) ’ 
using the identity 
Pre) (DES 


the last expression of (62) is obtained. 
(**) in may be omitted when not treating processes involving real mesons. 
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occupation number and therefore does not contribute to the diagonal part 
of 2, which we require. Accordingly it may be dropped. 
We have then 


1 
2.64 QO = à 4 EA) 
| È | i A ta He D, (Ha | ae 
We define 
| Vo = (A, + Do)aiag 


U(e) = (H, + Dy) 


non diag (e2) 


| U,(0) = (H, E Dy) non diag (02) ? 


U,(e) creates or annihilates an even multiple of 2 mesons, U(o) creates or 
annihilates an odd multiple of 2 mesons. 
By induction from (2.65) we define the following quantities: 


Un i (Un(e) an DI ’ 
U,,(e) == (U,-4(e) ns Bie at. diag (e2"+1) 9 
2.66 i 
(2.66) Uso) = UO PARERE 
1, fi è ISP, == Ur > v; b) 
i=0 
with 
n 1 
(2.67) DL, (0) MEA 


A ee) 


U,(e) creates or annihilates an even multiple of 2"*! mesons, U,(0) creates or 
annihilates an odd multiple of 2”*! mesons. 
At this point it can be shown that, if Q% is of the form 


Il 
2.6 Oy = PTE TT È ER ve Lia 7 Us, 
ee) da I i A = Ye. a UE) a> (GL (0) ( 7 (e) a (0) È 
the Q*P is of the form 
] 4 
2.69 AN — VU te) LU (0) | ; 
(2.69) © 00 î feno ( 11 + Unsi(0) 4(0)) a 


where Upis(e), Unislo), Yan have an exactly analogous definition as (2.66), 


a 
Uni 
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(2.67). Proceeding in this way for n large enough, if the matrix elements of C 


foe) 
become negligible, so that the series > v, converges, one obtains 


i=0 
= c ] : 
(2.70) Cal rier pee ; 
where Y is the sum of the series Vo. 
i=0 


The quantity Y is diagonal in meson occupation number and is assumed 
by BRUECKNER and WATSON as the definition of the two-nucleon interaction 
operator. 

The method here displayed can be applied to the case of the Hamilton- 
jan (1.23) 


Se; | oe RUE one 
Hy, = Me [etsite VEN ren she] CE es 


(the case of a Hamiltonian which is non-diagonal in nucleon occupation number, 
for instance the Hamiltonian (1.17) 


He = in| wY; tyPda 


requires a few small changes). Denoting by H, the meson-nucleon p-wave 
interaction term (which is linear in the meson field variables), by H, the meson- 
nucleon s-wave interaction term (which is quadratic in the meson field var- 
iables), the results are 


(2.71) = [AT ir HAE H, Wo HH lie ? 
(2.72) D = [Uo(0)(W — Hy — %— U(e))U(0)];ag + 


Restricting ourselves to the contribution to the potential arising from the 
exchange of at most two mesons, these two expressions are sufficient; regarding 
H, as a perturbation and expanding in powers of H, the result is 


(2.73) 7 = HOS HY? 4 HOS HS ee 


Il 
A 


(— (— ++ _ Cei + 
| A j Ave 7 Le elia A— HS He Hee 


A A 


The meaning of HT, HS, HM”, HS”, ete., is obvious. The last term 
exhibits the effects of multiple scattering. This term may in fact be written 
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in the form 


Ha 0) ~ He 4 


where @ is a solution of the equation 


o=l+-u@, 
4 


in which 


VA He yi HE 
À 


acts as a «potential » for the scattering of the meson emitted by HY and 
absorbed by H‘”. 

The evaluation of the multiple scattering effects is the main advantage 
of the Brueckner and Watson method; according to the evaluation of B.W. 
these effects are negligible for «> 0.5 and therefore the last term may be 
written with a good approximation 


I i neue 
= HS TWP + HSS SH Hel A. 


In the static limit in which A ——H, the expression for the one-meson and 
two-meson exchange potential respectively are 


ree 

9 =" 2 : (1). 7@)). 

(2.74) Ve = ee (3 M) ae (ne) 
il Les 3 3 \| exp [— x] 
ICE (i | x | =| pa È 

i IO TV 
D} î = — 22 2 = = Vi 20 n 
2 Pi dir Ch) (4) na? ca 


2 ] æ\? : 
+ pe? a (; La) (4) 6 | Li exp [— 2a] — 


4 2 \2 
er (4 Vie ) {R, (0) + (0-6) R(x) + SiPa(e)} ; 


2 LO dn 23 va 4,029 
(2.754) R.(0)= = {ce To) È = =) (20) + È Er =] K(20) 4 
3 (D. (2) ! 3 dI K,(2) + SE bu Ki, (a) exp[- 2]} 
+ ( (T 14 )) Pr | 73 + ae he phe E, a) pe, 0 ’ 
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TT Bra ne 
Ct 2 > (1) (2) 2 2 K T 1 1 K (a) ex [— 
soa (ac) er. K;(x) + Ta gi 0 xp "216 
E 2 io DIE 
(2.75c) R;(x) = AE ot 3 K, (2x) + ES K(2%) 


5 1 . 
Lian (3 — 2(¢- 7) IE + 5 È. À) K,(x) + (= + = Ki(r)| exp [— all ; 


3 x? Nhe Hp a % 


It must be noticed that, as can be seen from eq. (2.70), the interaction 
operator Y° of B.W. acts as a «potential» in a Schrédinger-like equation 


(W—H,)P|®> = ¥ (W)P|®> 


for the two bare nucleons (*) and depends on the total energy W of the system. 

To the procedure of B.W. therefore the same objections must be moved 
as those stated when discussing the Tamm-Dancoff method (*). 

GARTENHAUS (*) obtained expressions for the potential similar to those of 
B.W. using pv coupling and a fixed and extended source, 7.e. the Hamilton- 
ian (1.22). The expressions of the Gartenhaus potential are in agreement 
with those of B.W. for À — 0 except for the factor v(k), the Fourier transform 
of o, which is chosen, 

k2 
eral 


24. Method of the Bethe-Salpeter equation (28:61). — This method is based 
on the relativistic two-body equation derived in the framework of field theory 
by Bethe and Salpeter (92). The problem lies in the transition from this equa- 
tion to a three-dimensional equation, which may be interpreted as a Schrô- 
dinger equation for the two nucleons. This reduction, which involyes a tran- 


v(k) = exp 


(*) This is a consequence of (2.58) where |y)is an eigenstate of the free Hamil- 
tonian and of the definition of Q,. 

(**) Its actual equivalence with the Tamm-Dancoff method has been proved by 
KLEIN (59). 

(°°) A. KLEIN: Phys. Rev., 91, 1285 (1953); 94, 195 (1954). 

(6°) S. GARTENHAUS: Phys. Rev., 100, 900 (1955). 

(61) A. KLEIN: Phys. Rev., 90, 1101 (1953); 94, 1052 (1954). In the following two 
articles the calculation of higher order potentials and the question of convergence are 
discussed: A. KLEIN: Phys. Rev., 91, 740 (1953); 92, 1017 (1953). 

(92) E. E. SALPETER and H. A. BETHE: Phys. Rev., 84, 1232 (1951); E. E. SALPETER: 
Phys. Rev., 87, 328 (1952). 
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sition to a one-time equation and, subsequently, a transition to a four- 
component wave function, has been examined in detail by MACKE (°°). 


Let 
(2.76) P= RP 


be the symbolic form for the B.S. equation for bound states, where ‘7 ='7 0,6) 
is a 16-component spinor and Æ is an operator which must be built up with 


the Feynman graphs, according to the well-known rules. 


Let us introduce a one-time wave function y= #0, and let us define an 


operator J which, when applied to Y yields y: 


(2.77) y = JY 


(the effect of J is therefore to equate the two times #4, #® of Y). 


For y an equation of the type 


(2.78) ee) 


analogous to (2.76) will hold. 


The eq. (2.76), (2.77), (2.78) may be written in terms of relative co-ordi- 


nates or relative momenta. 
In the former case 


ee ORC 


Y = Po ol(%, 0) . 
In the latter (assuming P= p+ p® = 0) 


PE ao) 


ni 


p= vp W). 


A representation for J easily obtains in momentum space. 


the operator J to both sides of the equation 


] .P 
(2.79) PL Lo) zine» oo iF Dn) P(X, Do) Do 
one gets 
if 
(2.80) pla) = 5 | (x, Yo) Apo; 


(63) W. Macken: Zeits. Naturforsch., 8 a, 599, 615 (1953). 
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hence 
(ce) 
J à d 
2.8 J > — Apr - 
a 2arihi Po 
00! 
Any summable function of p, is a right inverse of J, but is not in general 
a left inverse. Let in fact f(p,) be a function of p, such that 


I È 
sain [Pol ar — 1"; 


—œ 


then 


I ao 
Jf(Po) p(x) = sain | fo) (2) A = W(x), 
while 
fot P(x, Po) = (Do) va) F P(x, Do). 


Therefore the reduction of the B.S. equation to a three-dimensional equa- 
tion is affected by the non-uniqueness of the right inverse of J. Nevertheless, 
in the case of the static interaction, 1.e. when retardation effects are neglected 
in the transmission of the interaction, the right inverse of J, J,*, turns out 
to be uniquely determined as that operator which transforms w,, into W,, 
(the meaning of w,, and Y,, will appear clear presently). 

The result is 


(2.82) Te ho(B®(p) + E®(p)— W) 
ae Om e 0 === 


Ep) — EW — om)(Ep)— 2 W + eps) | 


Expression (2.82) may be taken as a reasonable Ansatz for J7!. Once J>* is 
determined by this Ansatz, the interaction operator in the equal-time equation 
may be formally derived. The following equation for K is obtained 

(2.83) K=KJJT+4, 

thus 

(2.84) Y= KJ, ‘y + AP, 


from which 
(2.85) P= (1—A)IKJ y ; 


hence 


(2.86) y = JA — A) KJo ty 
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and the operator g of (2.78) is therefore given by 
(2.87) ga J(1 — Ad) KI". 


For the further reduction to a four component equation, eq. (2.78) is first 

transformed from the spinor indices «« to the indices ee, which refer 

to free-particle energy and spin eigenstates, using the normalized free-particle 

solutions of the Dirac equation w,.(p)=w and w.(p)=u' as unitary matrices. 
With the transformation 


peut gi= ua qu D y 
equation (2.78) becomes 
(2.88) SIE 
where 
/ ti 
D = Yet) - 


Denoting by y’, , the eigenfunction corresponding to positive energy states 
of both particles, w, the eigenfunction corresponding to a positive energy 
state of one particle and a negative energy state of the other, etc. (each of 
which in turn has four components corresponding to the four spin states of 
the two particles), the eq. (2.88) breaks into a set of four simultaneous equa- 
tions: 


I [ fr I I 
= i 
Wie mm et Vi > I+ +,nVn 9 


7 


(2.89) ; PU 
Yy > Gy, + +Ÿ, ok i > Grin’ Pn! 5) 
n 
(with an obvious meaning of 7). 
By successively inserting y, from the second equation into the first, an 
. I . . 
equation for wy, , 18 obtained 


(2.90) TU 

with 

(2.91) f = 4. RIPARA. + 5: 9, ADE It E > g + +n Inn’ In, ++ A ee 
n ULI 


Using the appropriate well-known analytical forms of the graphs, f is found 


Uni 


Uci 
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to have the following form 


(2.92) = (2 E(p) — Wij ay de 
with 


B(p) = Va + pi. 


Eq. (2.90) thus takes on the form of a Schrédinger-like equation in momentum 


space: 
if ] 7 if x I I ti 
(2.93) (2E(p) — W)y,.,(p) = aI? (p, p': W)w.,(p') dp 


The operator f (and therefore 7) is built up by summing the contributions 
from all the permutations of points of all Feynman graphs, omitting all « prod- 
uct graphs » (*), de. 


(2.94) je Dalene 


where f,. is the contribution from a Feynman graph with a given topological 
structure and the points ordered in a given succession x; the prime in the 
summation symbol stands for the exclusion of « product graphs ». 

The potential V therefore has the form of a power series expansion in the 
coupling constant. By rearranging this expansion so as to have a sum on 
renormalized Feynman graphs, one obtains an expansion in the number of 
exchanged mesons. (This assumes that the theory is renormalizable and the 
sum is convergent.) 

The one-meson and two-meson exchange potential is obtained from the 
second and fourth order renormalized graphs respectively. 

The calculation of this potential, at the static limit, and its expression in 
terms of experimental data on pion-nucleon scattering (using the dispersion 
relations of fixed source meson theory for p-wave and s-wave meson scattering, 
along the lines established by Miyazawa) has been outlined by KLEIN and 
McCormick (8:64), 

We quote here an approximate evaluation of the static two-meson exchange 
potential (the one-meson exchange potential has the now familiar form) done 
by K.McC. for the purpose of exhibiting the order of magnitude of the « cor- 
rections » to the B.W. potential: 


(2.95) Vi == Vip + Vans + Fe 


(0) Given two graphs f, and f, we mean by «product graphs» all the graphs. 
Bos hile. ft’... which may obtained by matrix multiplication of the graphs f, and f,. 


iis Ja 
(4) A. KLEIN and B. H. McCormick: Progr. Theor. Phys., 20, 876 (1958) 
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. AR à U ‘ Gi à ) | 
(2.95a) L DD NE ue? (5) (a) {Ri(x) “ire (oo) R(x) "n Sqoft3(4)} a? 


43 (6-6) (È Br A K, (2x) + & "n =) K,(20)| — 
i: ie + Fi) fm + (a =) K,(22) 


x 


5 2 2 ax Don 
moitie u (4 \ XP 5 22] 
ut (4) Ashe we 


XL 
6 12 10 4 
+ += +1 


at ae aL? 
JON PAIS ee eNOS Dee lls 6 DURS 
== 2, ue? (4, fen È (6 D'or ane te: i — = Di x x x 
4 exp [= 2%1| 6 12" 8021 
— = (qq) 8, ome means Sand 
3 {E ) x? È tw Li | 


2 
| exp [— 2x] + 


7 TE A PS 
2.95b le |e -- 
ee) "pa Sine Gé) A | oe? 


u \*( gi \?8 Lo ANB Oi eee: 5 
no a cs = (0-2) (pr + ga) (20) + gs Ke(20)] + 


si STILI gi \? 12 
Ara we (Lr) (eA TT 


a Na (2%) 
2.95c VOTE (E LE | A 
a > opt (47) CE TNT 1 


si, onc? | pi )( gi ir (10-10) a _ Kol22)| 
ré x ; 


DV | AAC Tee a 


The expressions (2.95, a, b,c) are, respectively, the contributions from 
p-wave, from p-wave and s-wave interference, and from s-wave; in each of 
them the first term is what would be obtained treating the meson-nucleon 
scattering amplitude in Born approximation and is in agreement with the 
corresponding term of the B.W. potential; 4, and 4, are parameters which 
measure the weight of the deviation of p-wave scattering from Born approx- 
imation; 0; and 0, are s-waves damping factors, i.e. they play a role similar 
to that of the factor 7 introduced in (1.23). These parameters are determined 
on the basis of the analysis of experimental data of pion-nucleon scattering; 
using the accepted values of meson-nucleon scattering amplitudes, K.McC. 
adopt the following values: @, — 0.01, = 0.57, À =2.9(M}u), Age 10M): 


te 
+ 
Uni 
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3. — Remarks on the various theoretical potentials. 


31. Static potentials. — We turn now to a description of the behaviour of 
the static potentials obtained by the various authors. 

The potentials are plotted in the various diagrams for triplet even states, 
singlet even states, triplet odd states, and singlet odd states, respectively. 
For the sake of comparison the diagrams of Fig. 1a), b), c), d) exhibit, together 
with the complete potential V=V,+V, of B.W. the one-pion exchange po- 
tential. 

The asymptotic behaviour of the latter, as can be seen from the expression 


V AGS C(TD- TO) Le Dog) ES 1 È | 
ARE te (eae DH ES PES SRE 


exp [— #] 


HA 


? 


which is consistently obtained by all authors, is of the type e-® whereas the 
asymptotic behaviour of the two-pion exchange potential is of the type e * 
as was to be expected. The behaviour of V, is recorded in the diagram for 
G?/4ahe = 0.08, a value determined from the analysis of deuteron and nucleon- 
nucleon scattering data (*). 

From the diagrams three features must be pointed out: 


1) the central part of V, is attractive but very weak; 


2) the tensor part of V, has the sign favourable for the deuteron electric 
quadrupole moment and is sufficiently strong; 


3) the corrections which the complete potential V=V,+V, brings to 
the potential V, are relevant for 0.5<x<1 and negligible for x >1.5 (as 
was to be expected), but are not wholly negligible for 1 < x < 1.5 (*) (*#*). 


(*) This value is in good agreement with that obtained from the analysis of meson- 
nucleon scattering. The experimental verification of the one-pion exchange potential 
will be the subject of the second part of our article. 

(**) The validity of these considerations is not limited to the B.W. potential, which 
is here taken only for argument’s sake. 

(**) A preliminary evaluation of the three-pion exchange potential for the purpose 
of estimating the nucleon separation at which this becomes comparable with the one- 
plus-two-pion exchange potential, has been done by Macmpa and SENBA (85). According 
È their calculations the three-pion exchange potential is estimated to be negligible 
Oe wr Oe 7c 


(95) S. MacHIDA and K. SENBA: Prog. Theor. Phys., 13, 389 (1955). 
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T.M.O. potential. The T.M.O. (2.43), (2.44) is plotted in the diagrams of 
Fig. 2a), db), c), d). 

In the triplet odd states, besides a strong tensor «attraction », which is 
suited to account for the deuteron quadrupole moment, it yields a strong central 
repulsion; this central repulsion is an unfavourable hint on the T.M.O. po- 
tential. It should be noted however that the modulus of the central part is 
less than the modulus of the tensor part for x > 0.68, while it is larger for 
æ < 0.68, 7.e. virtually in the phenomenological region. Therefore a treatment 
of deuteron and of low energy scattering in triplet even states with the T.M.O. 
potential for x > 0.6 and a suitable phenomenological adjustement for x < 0.68 
might not be impossible. On the other hand if for the phenomenological inter- 
action in the inner region a reasonable «hard core » is assumed (according to 
the attitude most common today) the T.M.O. potential is faced by serious 
difficulties. 

In the singlet even states there is a strong attraction, and this is a favour- 
able indication regarding the low-energy singlet scattering parameters. It 
should be remarked however that the value of the T.M.O. potential around 
æ—1 might be too small to give a satisfactory value for the singlet effective 
range. 

In the triplet odd state there is only a weak potential in the region #> 0.7; 
this might be a favourable feature for the explanation of the observed angular 
distribution in high-energy n-p scattering. 


975 


B.W. potential. Its expression is given by (2.74), (2.75) and its behaviour 
is plotted in the diagrams of Fig. 3 a), b), c), d) for different values of A. It 
can be seen that in singlet even and odd states the potential is not very sen- 
sitive to the value of À, while in triplet states it is rather sensitive to the value 
of À This sensitivity would appear to be more sharp for small values of A. 

In the triplet even state the central part of the B.W. potential turns out 
to be radically different from the T.M.O. potential. A comparison with the 
T.M.O. potential can be made from the diagrams of Fig. 4 a), 5), €), d) where 
the B.W. potential is plotted for 2=0. It is seen that the potential for 
singlet even and odd states, the central part and the tensor part of triplet even 
states are in agreement with or not very different from the corresponding 
T.M.O. potentials; the central part of triplet even states exhibits a remarkable 
attraction instead of T.M.O.’s repulsion. 

This difference has been widely debated. It stems from the term 


G2, \22 exp [— æ][4 + 42 +27 _ 
= = (1). 472) 3a = © K x) 
8V, = pe? (sas) = ar) x i, (0) + 
2 3 2 2a + wu? _ ; 
+ (2 + 2 + 208) K,(a) +3 (00-00) (a + o) K,(a) + n k (2) - 


FR Ba + a? _ 
ay Si: (a eG) Pe K(2))|. 


d 
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According to B.W. and other authors (5166) this term is the first of a series 
of terms coming from an expansion of the non-adiabatic part of the one-pion 
interaction. An expansion of this kind might be expected to be not very 
fastly convergent; hence it might be unjustified to consider only the first term. 
The contribution from the whole series (without expansion) is evaluated by 
B.W. to be negligible for a reasonable deuteron wave function. Hence the 
omission of the term in question from the expression of the potential. 

Other authors (33435,51,52) object to B.W.’s point of view that this proce- 
dure fails short in the case of the neutral scalar theory, inasmuch as the exact 
solution, for the fixed source, which is known in another way (*), requires the 
inclusion of the term omitted by B.W. 

According to the point of view of F.S.T. the difference arises from the 
different treatment of the probability operator 1/J*J. The B.W. potential is 
obtained when assuming the value 1 for the probability operator, while the 
T.M.O. potential is obtained when expanding the probability operator up to 
the second order in the coupling constant. 


Gartenhaus potential. The Gartenhaus potential is plotted in the diagrams 
of Fig. 5 a), bd), c), d). In the triplet even state the tensor part is « attractive » 
for the one-pion exchange potential, «repulsive » for the two-pion exchange 
| potential, about one third as large as the first. The result is an «attractive » 
tensor interaction which vanishes at the origin. The central part exhibits a 
weak attraction mainly due to the two-pion exchange; it becomes repulsive 
for æ<0.6. 

In singlet even states there is a strong attraction (with a maximum of 
70 MeV at «= 0.6, with G?/4ahe = 0.08 and ©, = 6) almost entirely due to 
the two-pion exchange; the interaction is then strongly repulsive for 7<0.4. 

In triplet odd states the central part is everywhere attractive; the tensor 
part is «repulsive » for x>0.6, becomes «attractive » for 4<0.6 and tends to 
zero for æ — 0. 

In singlet odd states there is a repulsion for x>0.75 and an attraction for 
æ< 0.75 so strong as to give a bound state 1P,. 

It must be remarked however that the characteristic features of the Garten- 
haus potential occur inside the inner region whose phenomenological character 
was not recognized by Gartenhaus, claiming to extend the validity of his po- 
tential up to the origin. 


F.S.T. potential. The treatment of the probability factor makes the F.S.T. 
potential (2.32) a form of potential intermediate between the T.M.O. and the 
B.W. potentials, as is also apparent from the diagrams of Fig. 6 a), b), c), d). 


(56) A. KLEIN: Progr. Theor. Phys., 20, 257 (1958). 
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While it does not differ very much from the T.M.O. and B.W. potential for 
the tensor part of triplet even and odd states and the central part of triplet 
odd states, remarkable differences both from the T.M.O. potential and from 
the B.W. potential arise in the other cases. In triplet even states the strong 
central repulsion of T.M.O. is greatly weakened, although it remains a re- 
pulsion; in singlet even states the attraction of T.M.O. and B.W. is notably 
weakened for «> 0.6 and actually becomes a repulsion for # < 0.6; in singlet 
odd states the repulsion of T.M.O. and B.W. changes into a strong attraction 
for #< 0.8. 


K.M.O. potential. The K.M.O. potential (Fig. 7 a), 6), c), d)) is in almost 
complete agreement with that of T.M.O., F.S.T. and B.W. for the tensor part 
in triplet even states; the difference from that of T.M.O., F.S.T. and B.W. 
is a little more marked for the tensor part in triplet odd states; in singlet even 
states and in the central part of triplet odd state the attraction, of T.M.O. 
and B.W. is further enhanced. In singlet odd states its behaviour approaches 
that of F.S.T. although somewhat damped, and therefore differs sharply from 
that of T.M.O. and B.W. On the other hand the central part of the triplet 
even state comes near to that of B.W. and is therefore very different from that 
OL NO: 

In the K.M.O. potential, where no pair suppression is assumed, there is 
also a certain contribution from the s-wave, which however turns out to be 
rather small, except for the charge singlet states. 


Klein and McCormick potential. The term corresponding to the contri- 
bution from the s-wave, like in the K.M.O. potential, from which this poten- 
tial (2.95) does not essentially differ, gives a negligibly small contribution and 
is therefore omitted by KLEIN and McCormick. Its behaviour once this term 
is neglected, is plotted in the diagrams for the singlet even and triplet even 
states (Fig. 8 a), b)). 

For the singlet even states the strong attraction of B.W. is sharpened by 
the corrective terms occurring in the potential of Klein and McCormick. In 
triplet even. states the corrective terms on one hand sharpen the central at- 
traction, on the other they weaken the tensor « attraction » of the B.W. po- 
tential. The central part thus notably differs both from that of T.M.O. and 
that of K.M.O. (This is obvious since KLEIN and McCormick adopt B.W.’s 
prescription regarding the controversial terms.) In simple terms it may be 
said that the this potential is to the B.W. potential what the K.M.O. potential 
is to the T.M.O. potential. 


32. Non-adiabatic corrections and (L:S) term. The greatest efforts have 
been directed to the derivation of the static nucleon-nucleon potential. Ex- 
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tensive ambiguities prevent at the present stage of the theory a systematic 
study of non-adiabatic corrections. 

Experiment has recently pointed to the convenience (*) of considering ve- 
locity dependent terms: (L-S) (5?) coupling and also, more recently, of the type 


(6+ p)(o": p) and 4[(o™- L)(o®- L) + (o®- L)(o™- L)] 


A certain number of studies have been done in order to derive from meson 
theory terms of this kind. Since the Thomas spin-orbit coupling term (9), 
which is of order (u/M)?, is too small to be satisfactory, it is clear that a spin- 
orbit coupling term of adequate magnitude may be derived only taking into 
account the retardation effect in the interaction between the two nucleons. 
The usual way to evaluate this effect is to expand the interaction kernel in a 
power series of p/Me & u/M. 

In this way ARAKI (5) was able to show that from the interaction kernel 
up to the second order in the coupling constant no spin-orbit coupling term 
of order u/M relative to the static term can be derived from pseudoscalar 
meson theory. Later (7°), carrying out the calculations on the interaction 
kernels derived in the framework of pv coupling theory by I. SATO (47), he 
obtained a spin-orbit coupling term 


; oo \exp[= 27] 3 
dii n fea (sar) a rt 


MI 


LCL S)E 


which is isotopic spin independent, since the terms depending on isotopic spin 
were evaluated to be negligible. 

High energy p-p scattering experiments recently have suggested that spin- 
orbit coupling depends on isotopic spin (71). More accurate evaluations of the 
(L-S) coupling term have been performed by SHINDO and NISHIJIMA (7), by 
S. SATO (?) (using the F.S.T. method and the T.M.O. prescription for the 
probability operator), by I. Saro, K. ITABASHI and S. Sato (74) (using the 


(*) We do not go here into the question whether velocity-dependent terms are neces- 
sary to account for nucleon-nucleon experimental data. 


G. Bruit: Phys. Rev., 51, 248 (1937); 53, 153 (1938). 
x. ARAKI: Prog. Theor. Phys., 6, 379 (1951). 

G. ARAKI: Progr. Theor. Phys., 18, 13 (1955). 

D. S. SIGNELL, R. Zinn and R. E. MARSHAK: Phys. Rev. Lett., 1, 416 (1958); 
. Buzz, K. D. Pyart, C. R. FisHER and G. Brett: Phys. Rev. Lett., 3, 264 (1959). 
*) M. SHInpo and K. NIsHIJIMA: Progr. Theor. Phys., 13, 103 (1955). 

*) S. Sato: Progr. Theor. Phys., 18, 457 (1955). 

74) I. Sato, K. IrABASHI and S. Sato: Progr. Theor. Phys., 14, 303 (1955). 
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F.5.T. method and the F.S.T. prescription for the probability operator) on 
the basis of pv coupling theory. As an example we record the results of 
S. SATO: 


G2 \2 / u 
eh eee APTE 
ee | Là (aaa) (el 
1 + +e? 
ja opps Di “exp [— 2a][3 + 2(10-7©)] yi} (Les 


As can be seen from (3.2), it turns out to be positive in charge singlet states 
and negative in charge triplet states (the spin-orbit coupling term derived by 
I. SATO, K. ITABASHI and S. SATO on the contrary is negative for both charge 
states). In the shell model of nuclear structure the terms depending on 7, 
seem to be negligible; the V,, of S. SATo would thus have the «right sign ». 

More refined calculations on the basis of the Klein and McCormick method 
have been done by TZoAR, RAPHAEL and KLEIN (7). OKUBO and MARSHAK (8) 
with a preliminary evaluation were able to show that quadratic terms of 
the type 


4[(o?- L)(o G L) + (o'?: L)(o'?:L)] n (o- a) L? 


may be derived from meson theory; they follow from the expansion up to the 
second order in w/2M of the interaction kernel up to the second order in G. 

The most complete evaluation of the spin-orbit coupling terms up to the 
order (@2/4afic)(u/2M) and (G°/A4xhc)(u/2M)° is, up to now, that of SUGAWARA 
and OKUBO (7). From their results there appears a relevant difference between 
the ps-ps and ps-pv theories. In the former the (L-S) coupling term, of order 
ulM arises from the one nucleon-antinucleon pair processes alone. In the 
latter the main term follows from the expansion of the spin matrix elements 
in the one-pion exchange interaction kernel. On account of its origin in the 
two theories, the (L-S) coupling term seems to be more firmly grounded in the 
ps-pv theory than in the ps-ps theory. Moreover it turns out to be generally 
of greater magnitude in the ps-pv theory than in the ps-ps theory 

On account of the various reasons mentioned, we must conclude that the 
treatment of the spin-orbit coupling terms is still inedaquate. In the first place 
in ps-ps theory the terms in question arise, at the lowest order in w/M, from 


(75) N. TzoAR, R. RAPHAEL and A. KLEIN: Phys. Rev. Lett., 2, 433 (1959); 3, 145 


(78) M. SuaawaARA and S. OKUBO: Phys. Rev., 117, 605, 614 (1960). 
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the nucleon-antinucleon pair processes and are therefore sensitive to pair sup- 
pression. IwADARE (77) showed that the survival of these terms is tied to a 
partial pair suppression; besides the obvious case of complete suppression, 
they also vanish in the case of no suppression because of the cancellation be- 
tween the two-pair term, previously obtained by KLEIN (7%), and the one-pair 
term. 

Secondly the series expansion in powers of p/Me might be badly con- 
vergent, if at all. Even if the convergence of the expansion is admitted, can- 
cellations might arise when developing up to higher orders. 

Recent very suggestive ways to explain the origin of the spin-orbit term 
have been devised by Gupta (7), who assumes the existence of a neutral 
scalar meson 0°, with a mass considerably larger than the pion mass (about 
2 u) and by SAKURAI (8°), who assumes a resonance correlation of three pions 
(or a vector meson). 


Conclusions. 


As already pointed out the one-pion exchange potential is identical for all 
authors, irrespective of the method employed for its construction and of the 
coupling (ps or pv) used. One is therefore tempted to attach a certain degree 
of certainty to the one-pion exchange potential and to submit it to the test 
of experiment. That it should be in principle possible is made plausible by the 
fact that there is a peripheral region in which the one-pion exchange potential 
dominates the interaction; effects «localized » in this region should therefore 
be determined by the one-pion exchange potential. The experimental veri- 
fications of the one-pion exchange potential will be the subject of the second 
part of our article. 

Unfortunately analogous statements cannot be made for the two-pion ex- 
change potential which is still plagued by countless ambiguities. This is not 
surprising in view also of the fact that the evaluation of the various corrections, 
such as multiple scattering effects and radiative corrections, is still in a very 
preliminary stage and therefore hardly significant. Notwithstanding this var- 
ious authors have attempted experimental tests of certain one-pion plus 
two-pion exchange potentials; we do not think that the situation as to the 


) IWADARE: Progr. Theor. Phys., 18, 189 (1955); 14, 16 (1955). 
) A. KLEIN: Phys. Rev., 90, 1101 (1953). 

9) 5: N. Gupra: Phys. Rev Lett, 2, 12491959). 
dead ye 


SUKURAI: The spin-orbit force and a neutral vector meson, EFINS-60-18 
(April 1960). 
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two-pion exchange potential is ripe enough to make these tests very meaningful ; 
we refer therefore to the original literature for these tests (81). 


We thank Prof. P. CALDIROLA for his interest in our work and encour- 
agement to us; Prof. A. LOINGER and Prof. M. VERDE for helpful comments; 
Prof. M. CINI for extensive suggestions and advice; Prof. G. M. PROSPERI and 
Prof. P. GULMANELLI for useful discussions. 


(81) For instance: S. Fust, J. IWADARE, S. Orsuxi, M. TAKETANI, S. TANI and 
W. WATARI: Progr. Theor. Phys., 11, 11 (1954) for the T.M.O. potential; K. A. BRUECK- 
NER and K. M. Watson: Phys. Rev., 92, 1023 (1953); L. HULTHÉN and M. SUGAWARA: 
Handbuch der Physik, vol, Xxx1x (Berlin, 1957) for the B.W. potential; S. GARTENHAUS: 
Phys. Rev., 100, 9(0 (1955); J. L. GAMMEL and R. M. THALER: Phys. Rev., 103, 1874 
(1957) for the Gartenhaus potential; M. Konuma, H. Miyazawa, and S. OTSUKI: Progr. 
Theor. Phys., 19, 17 (1958) for the K.M.0. potential. 

Overall comparisons for the various potentials have been made by Japanese re- 
searchers and may be found in several articles of the Progress of Theoretical Physics. 
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3. Errors in the machine. — 4. Analogy with the mathematics of quantum 
mechanies. — APPENDIX. 


1. — Introduction. 


A learning machine, called PAPA (italian abbreviation for Automatic Pro- 
grammer and Analyser of Probabilities) has been built in our Institute according 
to the suggestions of one of us (A.G.) (1). Very briefly it works as follows. 

A set of photocells (A-units) receive the image of the pattern to be shown 
as filtered by a random mask on top of each photocell. According to whether 
the total amount of light is greater or smaller than the amount of light falling 
on a reference cell with an attenuator, the photocell will fire a « yes » or «no» 
answer into the «brain» part of the PAPA. The latter is simply a memory 
storing the «yes » and «no » frequencies of excitation of each A-unit for each 
class of patterns shown, together with a computing part that « multiplies » 
(or «adds logarithms ») in order to evaluate the probability that an unknown 
pattern belongs to a given class. For details, the reader is referred to the papers 
quoted in footnote (1). 

The purpose of the present note is to discuss some of the mathematics 


involved in the machine. 


(1) A. GAMBA: Proc. IRE., 49, 349 (1961); G. PALMIERI and R. SANNA: Methodos, 
to be published; A. GAMBA, L. GAMBERINI, G. PaLMIERI and R. Sanna: Suppl. Nuovo 
Cimento, 20, 112 (1961); A. GAMBA, G. PALMIERI and R. SANNA: Suppl. Nuovo Cimento, 


20, 146 (1961). 
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2. — Algebra of patterns and masks. 


In discussing the PAPA machine we have to consider both patterns and 
masks. 

Patterns are the objects shown to the machine. We may assume for sim- 
plicity that patterns are quantized like the following (20-quantized) example 


CALE 
227 
| |G 


p 7 ZA 
Z/MEE 


where the black and white zones are those that in the actual pattern are, re- 
spectively, dark and bright. Masks are those that appear on top of each photo- 
cell as filter and can be represented by symbols of the same type as the pat- 
terns, with the only difference that the black and white zones will now be 
interpreted as opaque and transparent zones, respectively. We will consider 
only quantized masks, too. Continuous patterns and masks will be considered 
as limits of quantized ones. Patterns and masks will therefore be identified 
by the same kind of symbols. 

Using the convention of reading the successive black and white zones like 
the lines of a page and using the correspondence 


bright (transparent) — 1 


(2) 
dark (opaque) — 0 


we can use à shorthand notation. Thus in the case of example (1) 
(3) (10110100011110100111). 


This convention will also allow us to introduce a Boolean algebra which holds 
for both patterns and masks. 

N-quantized symbols like (3)—which may represent either a pattern or a 
mask—will be denoted in the following by small greek letters «, B, y,.... The 
« threshold of a symbol » is defined as the number of 1’s appearing in the symbol 
independently of the order and will be indicated as an index to the greek 
letter used for the symbol; thus the symbol in (3) will be indicated as O15 
(or Biz, Yo, …). We may further define the product of two (N-quantized) 
symbols as the ordered Boolean product of the corresponding quantized « points ». 
Thus for N = 6, we have for example 


(4) (100110) - (101101) = (100100) 
or shortly 
(5) 3B, = Va 
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The relevance of this definition of product is apparent: when in PAPA a 
pattern (x) is shown through a mask (f,) to an A-unit, the total amount of 
light falling into the photocell is equal to the threshold of the product (72). 
Hence the usefulness of defining such a product and the reason for calling 
threshold the quantity which is the relevant one to decide whether the phys- 
ical threshold of the A-unit will be exceeded or not. 

The following manipulation rules are easily proved: 


(6) do ali 


i.e. the symbols are «idempotent » or, as a physicist would prefer to say, 
« projection operators », 


(7) (a Bs)y = CAVE) , 
(8) ors = bsx, , 


i.e. the product is both associative and commutative. For an arbitrary f,: 
(9) Oy Pr = Pas 

(10) oP, = Lo 5 

i.e. there exists the «unit» symbol 

(11) gees (hited) == 

and the «zero » symbol 

(12) ot = (000...0) = 0. 


From the definition follows that the product of two N-quantized symbols 
with threshold r and s, respectively, is a symbol with a threshold ¢, such that 


(13) t, = min (7, s)>t>max (r+s—N,0)=t. 


A more relevant formula for products will now be obtained. Let S, be the 

set of all the (3 = n different symbols with threshold r. Let NAGISA 
r (N —r):r: 

be the product of two such sets to be understood in the sense that every symbol 


N\ (N 
of the set S, is to be multiplied by every symbol of the set S,. Then the | S | . | è | 


products may be collected into sets S, and the following theorem holds 


ts (N —?)! 7 
y . Y — me — — — — A 9 
(14) IN = 2 {r—t)l(s —t) UN —7r—s + t)! t 
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where the sum is to be understood as «the entire set S, appears as many times as 
indicated by its coefficient ». The proof of formula (14) is given in the Appendix. 


aN 
Taking into account the number (7) of symbols of threshold # and the nor- 
malization coefficient 


A N! ora N 
nai nes ee dh | 


we may interpret formula (14) by saying that by multiplying a symbol of thresh- 
old r by a symbol of threshold s, one has a probability 


16 ris! (N = r)'(N —s)! 
do; XE ei 


of getting a symbol of threshold t. 

We already mentioned that in PAPA a pattern—that without loss of gener- 
ality we may assume to have a threshold r< N/2 —is shown to a photocell 
through a mask of threshold s. When the threshold ¢ of the product is above a 
certain value # this A-unit is excited. The state of excitation of a given A-unit 
constitutes the information which is carried into the machine. If nothing is 
known about the nature of the patterns to be shown to PAPA, 7.e. for a general 
purpose machine, it is clear that maximum information—in the sense of Shan- 
non’s theory of communication—is obtained from an A-unit when the prob- 
ability of its excitation is 3 for a random pattern. Then the amount of infor- 
mation obtained is exactly one bit. This situation can most easily be obtained 
by putting the threshold s of the mask equal to N/2 and the triggering thresh- 
old # of the A-unit equal to 7/2. This follows from formula (16) which is 
symmetric in ¢ around the point t=7/2 for s= N/2. 

From now on the symbol corresponding to the mask of the k-th A-unit 
will be denoted as follows 


(17) Oya K) Lol) = Val D , 


so that the threshold N/2 is clearly shown. 

The above condition refers to an individual A-unit. However the PAPA 
is a set of n such A-units, A,, A,,..., An. To get maximum information not 
only from each individual A-unit, but also from the whole set of them, we 
must require that the different A-units be uncorrelated, i.e. that the proba- 
bility of excitation of the k-th A-unit, A,, be independent of whether any 
other arbitrary A-unit, say A,, is excited or not (2). Obviously this condition 


(°) This independence is also necessary eventually to evaluate total probabilities 
as products of elementary (uncorrelated) probabilities. 
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is Satisfied if the transparent part of the unit A,, covers exactly half of the 
transparent part of the unit A,, i.e. if the product of the masks of the cor- 
responding A-units gives a symbol y,,, of threshold N/4 


(18) ff) wall) = Val À) - 


A set of n A-units, such that any pair of their masks satisfies eq. (18), will 
be called an « orthogonal » set of A-units. The fixing of the threshold at N]2 
(expression (17)) will be referred to as «normalization » of the A-units. From 
the above discussion it is clear that an «orthonormal » set of n A-units is a 
«best» set of n A-units for unspecified patterns, i.e. for a general purpose 
PAPA, with fixed masks and thresholds. One may now ask what is the max- 
imum number » of orthonormal A-units for A-quantized symbols. Since an 
N-quantized symbol contains N bits of information and since every orthogonal 
A-unit picks up one bit, it is clear that n cannot be greater than N (3). Let 
assume we take two random N-quantized symbols of threshold N/2. What 
is the probability that they are orthogonal? Using expression (16) for 
r=s=N/2 and t= N/4, we can easily see—for example using Stirling’s for- 
mula—that as N increases this probability tends to one. A random set is a set 
of masks almost orthogonal. This justifies our choice of random masks in PAPA, 


(*) The exact number, writing N= 2%-b, where 6 is odd, is n= 24_ 1 (A. LEVI 
and G. Manuzio: private communication). b occurs with quantized points, since 
one cannot split an odd number of points exactly into two. Leaving aside this trivial 
«unphysical » case, what is significant is the fact that » is « practically » of the same 
order of N. A «complete » set of 15 orthonormal symbols for N= 16 is given below 
as an example: 


UA, 
[TEE 


226 A. BORSELLINO and A. GAMBA 


instead of an orthogonal set, apart from consideration of practical convenience 
one has to keep in mind if PAPA devices will ever be built on a (chemical?) 
microscopic scale. 


8. — Errors in the machine. 


For the sake of simplicity let us consider only discrimination between two 
classes a and b. Let us indicate with A,, A,, Ag, ..., Axy ..., A, the associative 
units. After the learning process has been completed, PAPA will have mem- 
orized the following Table I, where a,, b, are the (relative) frequencies of 
excitation of the unit A,, for the examples shown to the machine during the 
learning process in the classes a and b respectively (0<a,<1; 0 <b;<1) 
The corresponding frequencies of unexcitation will be denoted with bars; thus 


(19) Gy = 1— 4; DU De 
TABLE I: 
: } Kone he 
A, Fe Salt A pi 
Classes 
a ay As | As dx Ay, 
iad | 

b bi ben AD) by db, 


An unknown symbol X, belonging either to class a or to class b, and with 
an excitation function £, like the following example 


| A-units | 
(20) ia 
(eure LL SCURI 


| | | 
4, | A | 4, A | se 
dl 0 1 0 


a ae 


| ’ 
| | | 
will be recognized as belonging to class a when 

(21) CORTEI 

and to class b when the other inequality holds: 


(22) SOIA 


What is the probability that PAPA will make a mistake, when it has learned 
in the above sense? 
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Introduce the following abbreviation 


De ls ln Uni 


b = bal babe: da) 


the index £ indicating a particular number and ordering of quantities with 
and without bars, 2e. a particular excitation function; thus to consider all 
possible cases £ could be an index running from 1 to 2”. 

a, is the probability of occurrence of the excitation function Æ for patterns 
of class a; similarly for b,. 

(a, +b,)/2 is the probability of occurrence of the excitation function £ for 
patterns of class a and 6 when the occurrence of each class is equally likely. 
Let 


M, = max (a,, b 


x DE 


(24) 
Mz = min (a,, b,), 


E 
Then from the inequalities (21) and (22) it follows that M,/2 is the prob- 
ability of a correct identification—and m,/2 the probability of a wrong iden- 
tification—of a pattern with excitation function Æ when both classes a and b 
are 4 priori equally likely. The overall probability of a correct identification 


is then 
sta <P Al 
(25) R= ae 


1 My 
(26) ‘dr > 3 

E dd 
where of course 
(27) eee 


D] 


It would be interesting to have an explicit expression for (25) in the form 
(28) = Er à 0)" 


where n is the number of A-units and « and f are the standard deviation of 
the frequencies a, and b, respectively around the point } (assuming a normal 
distribution). Then by measuring experimentally « and f one could derive 
the number n of A-units necessary for a PAPA with a preassigned reliability À. 
This mathematical problem is dealt with by E. GAGLIARDO in another paper. 
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4. — Analogy with the mathematics of quantum mechanics. 


There is an analogy between the process of computing probabilities used 
in PAPA and the method of «second quantization » in physics. 
Let the excitation function Æ, of pattern X in example (20) be written 


in the form 
(29) MAO ME corel) 
| E.) can be considered as a «ket» of a system of fermions (the «occupation 


numbers » being only 0’s and 1’s) with n independent states. Introducing the 
«creation operators » A}, Ai, …, Al, we can write 


(30) (Ee AA ee 
where 
(31) LO> F005 On 50> 


is the «vacuum state». Similarly any other excitation function can be ex- 
pressed by means of a ket of the form (29) or (30). 

We can then consider the whole set of patterns of class @ or class b (see 
Table I) as the kets 


k 


| COCO 


(32) and 


[b> = [[ [b,Az + A — B,)]|0 , 
k 


respectively. Expressions (23), giving the probabilities that a pattern Y with 
excitation function Z, belongs to class a and b are then given by the quantum 
mechanical « transition amplitudes » 


(33) <E, |a) and QUANDIE 


respectively. The class kets (32) are obtained from the kets corresponding 
to the excitation function of each individual pattern of the class, by a process 
different from the quantum mechanical process of the « Superposition prin- 
ciple », due to the lack of coherence of our kind of superposition. A better 
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formal analogy is obtained if the correspondence is established with « density 
matrices », which can take into account also uncoherent superposition. The 
analogy then runs as follows. 

Introduce for each A-unit, say A,, the operator 


: ) 
(34) re i i | 
0 0 


the eigenvalue 1 corresponding to the excited state and the eigenvalue 0 to 
the unexcited state of the unit. One can then define a one-to-one correspon- 
dence between excitation function Æ, and density matrix o in the product 
spaces of the operators (34). Thus in the case of the example (20) the excitation 
function would be represented as follows 


ae Lt ON 20 OF fl ON JE 0 0 
(35) D = ale 2 { ne 
= OO OMe 00) AL ONO hy 


the density matrices 


pe On 0 de 0 | Rg. 0 Un One 
ea e OC Oe bn) E Ud of 


(36) and 
> È 0 b, 0 | bs 0 fe 0 
Dr = o s Inv È 
si OM pi DIO, i Nb. 0 te 


would then represent the density matrices for the classes a and 6 respectively. 
In this case, however, they are obtained as superposition of the density matrices 
corresponding to the individual patterns belonging to the class, where super- 
position is now defined exactly like in quantum mechanics. The probabilities 


(33) are now given by the traces 


(37) Treo) and Tr(0,0); 


respectively. 

We wanted to point out this analogy with quantum mechanics, not only 
for the sake of mathematics, but also for its physical interest. These consi- 
derations in fact suggest that quantum mechanical PAPA devices of sub- 


microscopic size may eventually be discovered. 
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APPENDIX 


Proof of formula (14) of the text. 


Consider the set Sy_, of all the (N-quantized) symbols of threshold N — 1. 
They are N in number and they form a « generating set» for the totality of 
the symbols, since every symbol of threshold 7 can be obtained as a product 
of N — r such symbols. 

The following special case of formula (14) is easily proved by inspection 


(Ge) Sy Sr = (N r)S; (N “i 1)8:= ’ 


where the convention S_,,—0 is introduced here and in the following to 
deal also with the special case r= 0. 

In fact there are N — r symbols of the generating set having their « zero » 
in those N — y points where every symbol of threshold r has a zero. These 
products leave the symbol unchanged and this accounts for the first term in 
the right-hand side of eq. (1'). The other + symbols of the generating set will 
put a «Zero » in place of a « one » and so transform the symbols from threshold r 
to threshold r—1. The totality of symbols of threshold r — 1 (and therefore 
the set S,-,) must be obtained by symmetry: the number of times this set 
appears is then determined by the requirement that the number of products 
on both sides of eq. (1’) be the same. 

Multiplying formula (14) on both sides by Sy and using eq. (1’), one 
obtains after collecting terms 


2! Se 5 = id ee, ere alt, et MILE Pe ee oh 
O A ORNE 20 
hi (N a a Pr 
LIA CLT ERIE: r—s 1)(N E Lasa 
N ne 


= > —— = 54 + 


ing, (F —i—1)K{s—1T'(N—r—-s Er 7)! 
(M = t,) (7 — 4) 


ie, ee 
| = Si, 


(r —t,) (8 —t,) (NW — rs +4) 1(N — 7 +1) 
| (Ve toy! 


rh) —t,)!(N¥ = r— s +4,)'(N —r 41) 


Sa 
The last term in eq. (2') is equal to 


0 itp" 0), 


(2N —r—s+1)! 8 
N° ISS E 


ift=r+s—N, 
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the latter expression being nothing but the value under the sum in eq. (2') 

for t=r+s—N—1. We can therefore suppress this term and compensate 

for it by shifting the lower limit in the sum from ¢, tot; — max (r — 1+ s— N,0). 
Similarly the term in S,, is equal to 


0 a = iP 


Ce) SE, 
eu ua LS 
(r-s-1)N—r +1)!" i 


eu = 


This term can also be put into the sum by shifting the upper limit from #, 
I 

to tf, = max (r — 1,8). 
Therefore we can write 


54 AS tite S 
di de ‘CZ Enpa 


4 


? 


which is nothing but eq. (14) with r— 1 in place of 7. 
Since the formula is true for r= N —1 and s arbitrary (see eq. (LI) 
the theorem is now proved by induction. 
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SUPPLEMENTO AL VOLUME XX, SERIE X None SLI GI 
DEL NUOVO CIMENTO 20 Trimestre 


On the Evaluation of a Formula for the Errors 
of a Learning Machine. 


E. GAGLIARDO 


Istituto di Matematica dell’ Università - Genova 


(ricevuto il 7 Luglio 1961) 


1. — In a recent paper by A. BORSELLINO and A. GAMBA (1) the compu- 
tation of error in a learning machine gives rise to a statistical problem which 
can be stated in the following way: Let n be a given integer number (very 
large) and let 


CIRO ON, 


be 2n real numbers with 0<a,<1, 0<b,<1. 
Let 


y = A exp [— (a@—})?/a?], y = B exp [— (x — 3)?/B?] 


be the statistical distributions of the numbers a,, 6; in the interval (0, 1) 
where x, 6 are given constants, while A, B are the factors of normalization: 


L 


ue | | exp [— (@ — pfea}ac) ul | | exp [— (è — 3)?/p2] de) 


We denote by 
{PD} = {Psa (2), sory Ps,n(@)} ’ (s =], 2, 3, 0°) 2”) 


the totality of the sequences of n functions each of which can assume one or 


(*) A. BORSELLINO and A. GAMBA: Suppl. Nuovo Cimento, 20, 221 (1961). 
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‘ 


the other of the two following determinations: 


i 
> ‘ K 


| æ, | 
à Ps.i(X) — 
; | l_-a 
pai For every s we consider the numbers: à 2 
’ : 
| cm i 
M, = max [9,1(41)°...° @en(Gn), P2,1(0i)"2 Dx (07)}, sea 
ms = MIN | @s,1(1) wees "gold Psa(0:) Peso (Dr) A ENS 3 5 
A 
“ogy RS 
Our problem is to calculate for given values of a, b, n, the two following ar 
numbers: "4 
= + | 
at 27 fi on Pia Di 
RES My e=4)>m, î 
s s ry 
It can be noted that R--e—1, in fact for each sequence of n numbers " 
i", C}....-, C; one. has on 
| 192 tf, 
: Y Parle) a Pan(On) = 1 i 
as can be easily proven by induction with respect to n. "LE 
It can be also proven (see Section 4) that for n ++ co one has e 0. The 
. . . . > Ah, 4 
‘most important aspect of our problem is to determine the size of ¢ for given si 
Bo 
-; UD x, B. ny 
Consider, for every fixed 0<y<n, the following 2" vectors of the (x, y) È 
fe plane: e 
Aa n 5 7 , M È . un CP ui 
5, 0 = (1+@,,(a;)*-..°9,,(4,), 1-9, ,(0,)---<°@,,(8,)) « LATE 
J 


When speaking of vectors, it is always understood that they are applied 
i «at the origin. | ae 
For y=0 all these vectors coincide with the vector (1,1). However, for»; 
v=1 one has two positions generally distinct one from the other: (a, di), 
- (1—a,, 1—b,), which together, by addiction, give the vector (1,1); and as » 
_ inereases every position breaks down into two others. The problem now is 
to study the distribution of the vectors v”, in fact the value of 2e is obviously 
“given by the sum of the abscisses of the vectors vo which are found above the: 


8 


_ diagonal y —x and of the ordinates of those vectors below. (By analogy we 


can calculate 2/ by inverting the roles.) We can also take into consideration 
the symmetry relative to the diagonal. 
| For » 1, the statistical distribution of the extreme points of vectors vi, 
Ù 
fr 2 | 16 - Supplemento al Nuovo Cimento. 
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is given by the function 


AB exp [— (w— 9)*/a? — (y—3)*/P°] for O<w<l, 0<y<Tt 
f(æ, y) = 


0 elsewhere . 


The statistical distribution of the vectors v' which derive from the de- 


s 
(0) 


composition of a vector e” = (£, 7) is given by 


i Dae 
EE 
RS) 


Our question can now be reduced to a problem containing a single variable 
by means of the following observation: The evaluation of e (and of À) remains 
obviously unchanged if during the decomposition, for © vectors with the 
extreme points coinciding at the point (x,y), (© real positive), we substitute 
x vectors with the extreme points coinciding in the point projection of (a, y), 
from the origin, on the line x=1, that is in the point (1, y/x). 

By systematically using this observation, it is enough to consider only 
vectors with extreme points on the half line x=1, y>0. For »—0 there 
is still coincidence with the vector (1, 1). 

Going from » to y-+-1 the statistical distribution of the vectors (1, y,) coming 
from © vectors (1,7) can be obtained by integrating along the line y = ya, 
with do = V1+y2da, the function 


multiplied by wr (because of the preceeding observation), by the jacobian 
o=xv1+y?, and by 

dd Rr 

dYo sal + Un 


(thus justifying a successive integration with respect to y, instead of #), and 
is then given by 


h(n, Vo) 


Il È 1\2 / Y x L\2 3 
DAB a? ex æ 12 UAL eM 
©) > 7 re exp 7 (a 5) i a | A 3) JR da x 


0 


where 


h(y, Yo) = min [I i 4 
Yo 
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Consider the sequence of functions # (y) constructed in the following way: 


(Fo(y) = o(y — 1), (Dirac)) , 


RO, y) 
À > ; IRE 1\2 
By) = AB | Vex pa È = ji tige (ye = 5 |p. da, 
û ee _ 
asa ho, 4) 
; di} Les: VM 
F,11(9) =A B| F (1) ni x? exp (e zi je È i 3) JB: da. 
0 0 


The solution of the problem considered at the beginning is obviously given by 


+o 1° 
e= 2" [Faw dy = 2" |y Fy) dy. 


0 0 


2. —~ We come now to several considerations useful in the numerical treat- 


ment of our problem. 
The points of the positive y-axis can be sketched by the following values: 


pales tia hae Oa PT pe ORR ah 


where g is sufficiently near to 1 and k is an integer sufficiently large. More 
precisely we shall identify with y=g° the points of the segment: 


BG 1) gn), (i=—k, ..., +h), 


For every 0<v<n the function 2#”(y) can then be sketched by 2%+1 


numbers: 


gv) gr 90) 
CER OI SOON 


that for » —0 have the initial values 
SA AM Og Orem lin AVE PAU 


and for » —0 can be constructed according to the following rule: 
Me. +k 
Q(v+1) Gv) 
Ù = DI 120, 
j 


1® 
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where 


min [1,9 LI 


I ek 
i; = AB (a = a q' | exp [— (æ — 4)?/a? — (gia — $)?/B?] da. 


0 


The limits of integration, if a, b are sufficiently small, can be substituted 
By, 005). CO. 


1 : The solution of our problem is given by 

hoy 

ti gn) 9,02) 

las Vo i ¢ " da 

es == SO SIN = » go Joe 

7 i>0 “i i<0 # 

eo 

Di It is easily seen that the values 9 (putting bY — 0 for 7<—k and for 
cat i> k) furnish the coefficients I, of the twice repeated transformation. Con- 


tinuing in this way one arrives at the following rule: 


ne a i 


—k...k 
“af (1) gi) 
pe = > 9 9° 
a 


moi which is very convenient to use, after having calculated 9, if n is very large. 

CR 3. — To diminish the error inevitably introduced by numerical calculations 
ollowin servations are useful in progra ing: 

the following observations are ful programmin 


a) At every step it is convenient to normalize the numbers 9°, dividing 
by $ > 9"(1+4'). 


Cela Y 


b) The value of q can be changed (more than once) in the course of the 


à procedure: initially very near to 1, then when 8%, 9% begin to become sig- 
nificant with respect to 9”, it is useful to substitute q? for q; then the values 

} dy? must be placed in the boxes reserved for 9”, while the 0” with i odd can 

: be forgotten or, even better can be divided into two equal parts and placed 

i in the two adiacent boxes before making the aforementioned displacement. 

c) Because of small numerical errors the quantities 
Dotata PT Lire 
da è ‘ i>0 D i<0 = 


LS 


which theoretically should remain equal, progressively differ one from another, 
slowly at the beginning and then, by consequence, more quickly. 


556 


A, 
+ 


A: 
x 


Gi 
= 


ee PEN aN [in net, Se déc N LI US ee } x s or 
si 299 xt * RE pe più Py w a 4 ‘ + 1 x PI } 1 7 
4 + he y 2 ti ; Cr ; ha ÿ ; 
7 Per ve ey Ea nr” ; y : 
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Fortunatenà the rule for the construction of 92" is invariant by translation 
(y) 


with respect to i: when displacing the 9” in the boxes reserved for the #5, 
_ the 9?” also undergo a similar displacement. 


It is possible to profit from this fact by making at each step a displace: 
ment in such a way to make the two quantities under consideration as equal 


as possible. 


4. — The method explained in the Section 2 proves also that (as remarked 


in Section 1) for n —> + co one has ¢-oo. (The treatment by means of in- : 


tegrals, as in Section 1, is left to the reader.) 
We can put k=-+ co and ¢ arbitrarily near 1. 
The numbers J” satisfy the condition: 


> (1+9)=1, 
thus 


3 ay Wy Hie 


z 


Let 


(VP); _ (7) g(r) 
a = 0; DE 2 
j 


it is enough to prove that for y ++ co 


The rule which gives the 77” is 


(ysl) —_ (v) 
Ni => K;_,n; » 


where 


Since there are at least two K,+0, it is easily seen that the center of the 


system of the weights 7” at the points i of the æ-axis, when progressing from 


y to v-+1, undergoes a constant tr anslation (which could also be a vanishing 


one) while the center of the system of the same weights placed at the points 


|i| diverges. 
Hence the above statement. 


#2 
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5. — A brief calculation made by means of the IBM 650 at the Centro di 
Caleolo Numerico of the University of Genova has given the following results: 


MABTE 


| 
| e(n) 


ML 
| a= B==1/10 | x=f=1/100 x— B—1/1000 
d — 92 4.44-10-1 4.94-10-1 4.994: 10-1 
92 = 95 | 4.23-10 | 4.92-10-1 4.992: 10 1 
23 : 94 3.94-10-2 4.88-10-1 4.988-10-1 | 
94 95 3.58-10-1 4.84-10-1 4.984: 10-1 | 
| 25 + 96 311-101 477.102 4.97 -10-2 
26 : 97 2.25-10-1 | 4.68-10-1 4.96 -10 1 
27 +98 | 1.94-1072 | 4.56-10-1 4.95 -10-1 
28 99 | 131-104 4.39-10-1 4.93 -10-1 | 
29 = 910 7.57-10-2 416-101 4.91 -10-2 
210 : 911 3.46 +102 3.86: 10-1 4.87 -10-1 
gi 92 | 1.14-10-2 3147-10 4.82 101 
212 918 3.40 - 107 | 2.99 - 1071 ANSE 
213 : 914 2.63-10-4 2.41-1073 4.65 -10-2 
gli : 915 | 1.15-10-5 1.78-10-! 4.51 -10-1 
915 + 916 1.39-10 7 1.16: 10-1 4.32 +101 
916 + 917 2.68: 10-10 6.39: 10-2 4.07 -10-2 
917 : 918 3.87: 10-14 2,72-10-2 8.74. "10: 
D18 = 919 1.43- 10-19 8.18- 10-3 3.31 -10-1 | 
919 = 920 | 3.00 -10-27 1.48-10-3 2.80 -10-1 | 
920 - 921 4.14-10-38 1.33-10-4 2.20 -10-1 
9212 922 : 4.43-10-6 1.56 :10-1 
922 923 3.58-10-8 9.70 -10-2 
223 - 924 2 3.94- 10-11 4.92 -10-2 
224: 925 = 2.57-10-35 1.88 -10-2 
225 : 926 | = 3.10: 10-22 4.84 -10-3 | 
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